PipMaker--a web server for aligning two genomic DNA sequences

Molecular basis of positional memory in limb regeneration

The amputation of a salamander limb triggers anterior and posterior connective tissue cells to form distinct signalling centres that together fuel regeneration1. Anterior and posterior identities are established during development and are thought to persist for the whole life in the form of positional memory2. However, the molecular basis of positional memory and whether positional memory can be altered remain unknown. Here, we identify a positive-feedback loop that is responsible for posterior identity in the limb of an axolotl (Ambystoma mexicanum). Posterior cells express residual Hand2 transcription factor from development, and this primes them to form a Shh signalling centre after limb amputation. During regeneration, Shh signalling is also upstream of Hand2 expression. After regeneration, Shh is shut down but Hand2 is sustained, safeguarding posterior memory. We used this regeneration circuitry to convert anterior cells to a posterior-cell memory state. Transient exposure of anterior cells to Shh during regeneration kick-started an ectopic Hand2-Shh loop, leading to stable Hand2 expression and lasting competence to express Shh. Our results implicate positive-feedback in the stability of positional memory and reveal that positional memory is reprogrammed more easily in one direction (anterior to posterior) than in the other. Modifying positional memory in regenerative cells changes their signalling outputs, which has implications for tissue engineering.

Evolution of Pineal Nonvisual Opsins in Lizards and the Tuatara and Identification of Lepidopsin: A New Opsin Gene

Many lizards (Squamata), as well as the tuatara (Rhynchocephalia), are distinguished among vertebrate groups for the presence of the parietal eye, or "third eye", a structure derived from the pineal complex containing a simplified retina with photoreceptor cells. The parietal eye expresses nonvisual opsins that differ from the visual opsin repertoire of the lateral eyes. These are pinopsin (OPNP), parapinopsin (OPNPP), and parietopsin (OPNPT), all being evolutionary close to visual opsins. Here, we searched over 60 lepidosaurian genomes for pineal nonvisual opsins to check for the evolutionary trajectory of these genes in reptiles. Unexpectedly, we identified a novel opsin gene, which we termed "lepidopsin" (OPNLEP), that is present solely in the genomes of the tuatara and most lizard groups but absent from other vertebrates. Remnants of the gene are found in the coelacanth and some ray-finned fishes, implying that OPNLEP is an ancient opsin that has been repeatedly lost during vertebrate evolution. We found that the tuatara and most lizards of the Iguania, Anguimorpha, Scincoidea, and Lacertidae clades, which possess a parietal eye, harbor all pineal opsin genes. Lizards missing the parietal eye, like geckos, teiids, and a fossorial amphisbaenian, lack most or all pineal nonvisual opsins. In summary, our survey of pineal nonvisual opsins reveals (i) the persistence of a previously unknown ancient opsin gene-OPNLEP-in lepidosaurians; (ii) losses of nonvisual opsins in specific lizard clades; and (iii) a correlation between the presence of a parietal eye and the genomic repertoire of pineal nonvisual opsins.

Codon usage patterns and genomic variation analysis of chloroplast genomes provides new insights into the evolution of Aroideae

Aroideae is an important subfamily of the Araceae family and contains many plants with medicinal and edible value. It is difficult to identify and classify Aroideae species accurately on the basis of morphology alone because of their polymorphic phenotypic traits. The chloroplast genome (CPG) is useful for studying on plant taxonomy and phylogeny, and the analysis of codon usage bias (CUB) in CPGs provides further insights into the intricate phylogenetic relationships among Aroideae. The results showed that the codon third position of the chloroplast genome coding sequence in Aroideae was rich in A and T, with a GC content of 37.91%. The ENC-plot and PR2-plot revealed that the codon usage bias of Aroideae was influenced by multiple factors, with natural selection as the dominant factor. Thirteen to twenty optimal codons ending in A/T were identified in 61 Aroideae species. Additionally, the comparative analysis of CPGs revealed that two single copy regions and non-coding regions were variable in Aroideae. Eight highly divergent regions (Pi > 0.064) were identified (ndhF, rpl32, ccsA, ndhE, ndhG, ndhF-rpl32, ccsA-ndhD, and ndhE-ndhG) , in which ndhE have the potential to serve as a reliable DNA marker to discriminate chloroplasts in Aroideae subfamily. Furthermore, the maximum likelihood-based phylogenetic trees constructed from complete chloroplast genomes and protein-coding sequences presented similar topologies. Principal component clustering analysis based on relative synonymous codon usage values (RSCUs) revealed that Calla was clearly deviated from Montrichardia and Anubias, and that Alocasia was closer to Colocasieae than to Arisaemateae. These findings suggest that the use of RSCU for clustering analysis could offer new theoretical support for species classification and evolution. Our research could provide a theoretical foundation for the chloroplast genetic engineering, taxonomy, and phylogenetic relationships of Aroideae chloroplasts.

A comprehensive review of peroxiredoxin 4, a redox protein evolved in oxidative protein folding coupled with hydrogen peroxide detoxification

Peroxiredoxin (PRDX) primarily employs electrons from thioredoxin in order to reduce peroxides. PRDX4 mainly resides either in the endoplasmic reticulum (ER) lumen or in extracellular spaces. Due to the usage of alternative promoters, a first exon is transcribed from different regions of the Prdx4 gene, which results in two types of mRNAs. The first type is designated as Prdx4. It is translated with a cleavable, hydrophobic signal sequence and is expressed in most cells throughout the body. The second type is designated as Prdx4t. The peroxidase activity of PRDX4 is involved in both the reduction of hydrogen peroxides and in the oxidative folding of nascent proteins in the ER. Prdx4 appears to have evolved from an ancestral gene in Eutherians simultaneously with the evolution of sperm protamine to cysteine-rich peptides, and, therefore, the testis-specific PRDX4t is likely involved in spermatogenesis through the oxidative folding of protamine. The dysfunction of PRDX4 leads to oxidative damage and ER stress, and is related to various diseases including diabetes and cancer. In this review article we refer to the results of biological and medical research in order to unveil the functional consequences of this unique member of the PRDX family.

Evolution of CEACAM pathogen decoy receptors in primates

BACKGROUND

CEACAM1 in leukocytes controls cell activation during inflammation. This and its expression in epithelial cells led to frequent independent appropriation of CEACAM1 as receptor by pathogens in humans and other species to gain host access and to downregulate its immune response. As a countermeasure, decoy receptors with CEACAM1-like pathogen-binding domains evolved. The granulocyte-specific human CEACAM3 endocytic receptor diverts CEACAM1-binding pathogens to neutrophils for internalization and destruction. The role of the glycosylphosphatidylinositol-anchored CEACAM5 and CEACAM6 which can also bind CEACAM1-targeting pathogens in humans is less clear.

METHODS

We analyzed the selection of CEACAMs to avoid pathogen binding and to maintain similarity between pathogen receptors and decoy receptors in 148 primate species.

RESULTS

Notably, functional CEACAM3 genes were not found in gibbons and New World monkeys. Interestingly, CEACAM6 in these primates exhibits similar high ratios of rates of nonsynonymous and synonymous substitution (dN/dS) in their pathogen-binding N domain exons as found for CEACAM1. High dN/dS ratios are indicative of selection for diversification typically seen in pathogen receptors. Human CEACAM6 is expressed on granulocytes and epithelial cells. Therefore, CEACAM6 could substitute for the missing endocytic receptor CEACAM3. In nearly all investigated primate groups also N exons of the epithelially expressed CEACAM5 exhibit selection for diversification. In African populations, five high-frequency polymorphisms are observed in the pathogen-binding region of CEACAM5 (I80V, V83A, I100T, I112V, I113T) with 3-4 polymorphisms combined in the same individual. These polymorphisms correspond to CEACAM1 pathogen-binding domain sequences.

CONCLUSION

The glycosylphosphatidylinositol-anchored CEACAM5 and CEACAM6 are under selection to maintain similarity to the pathogen receptor CEACAM1 in most primate species, indicating a function as decoy receptors.

Translation of genome-wide association study: from genomic signals to biological insights

Since the turn of the 21st century, genome-wide association study (GWAS) have successfully identified genetic signals associated with a myriad of common complex traits and diseases. As we transition from establishing robust genetic associations with diverse phenotypes, the central challenge is now focused on characterizing the underlying functional mechanisms driving these signals. Previous GWAS efforts have revealed multiple variants, each conferring relatively subtle susceptibility, collectively contributing to the pathogenesis of various common diseases. Such variants can further exhibit associations with multiple other traits and differ across ancestries, plus disentangling causal variants from non-causal due to linkage disequilibrium complexities can lead to challenges in drawing direct biological conclusions. Combined with cellular context considerations, such challenges can reduce the capacity to definitively elucidate the biological significance of GWAS signals, limiting the potential to define mechanistic insights. This review will detail current and anticipated approaches for functional interpretation of GWAS signals, both in terms of characterizing the underlying causal variants and the corresponding effector genes.

A Comprehensive Analysis of the Genomic and Expressed Repertoire of the T-Cell Receptor Beta Chain in Equus caballus

In this paper, we report a comprehensive and consistent annotation of the locus encoding the β-chain of the equine T-cell receptor (TRB), as inferred from recent genome assembly using bioinformatics tools. The horse TRB locus spans approximately 1 Mb, making it the largest locus among the mammalian species studied to date, with a significantly higher number of genes related to extensive duplicative events. In the region, 136 TRBV (belonging to 29 subgroups), 2 TRBD, 13 TRBJ, and 2 TRBC genes, were identified. The general genomic organization resembles that of other mammals, with a V cluster of 135 TRBV genes located upstream of two in-tandem aligned TRBD-J-C clusters and an inverted TRBV gene at the 3' end of the last TRBC gene. However, the horse b-chain repertoire would be affected by a high number of non-functional TRBV genes. Thus, we queried a transcriptomic dataset derived from splenic tissue of a healthy adult horse, using each TRBJ gene as a probe to analyze clonotypes encompassing the V(D)J junction. This analysis provided insights into the usage of the TRBV, TRBD, and TRBJ genes and the variability of the non-germline-encoded CDR3. Our results clearly demonstrated that the horse β-chain constitutes a complex level of variability, broadly like that described in other mammalian species.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state regulatory potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbor distinctive transcription factor binding motifs that are similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we show that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

Interspecies regulatory landscapes and elements revealed by novel joint systematic integration of human and mouse blood cell epigenomes

Knowledge of locations and activities of cis-regulatory elements (CREs) is needed to decipher basic mechanisms of gene regulation and to understand the impact of genetic variants on complex traits. Previous studies identified candidate CREs (cCREs) using epigenetic features in one species, making comparisons difficult between species. In contrast, we conducted an interspecies study defining epigenetic states and identifying cCREs in blood cell types to generate regulatory maps that are comparable between species, using integrative modeling of eight epigenetic features jointly in human and mouse in our Validated Systematic Integration (VISION) Project. The resulting catalogs of cCREs are useful resources for further studies of gene regulation in blood cells, indicated by high overlap with known functional elements and strong enrichment for human genetic variants associated with blood cell phenotypes. The contribution of each epigenetic state in cCREs to gene regulation, inferred from a multivariate regression, was used to estimate epigenetic state Regulatory Potential (esRP) scores for each cCRE in each cell type, which were used to categorize dynamic changes in cCREs. Groups of cCREs displaying similar patterns of regulatory activity in human and mouse cell types, obtained by joint clustering on esRP scores, harbored distinctive transcription factor binding motifs that were similar between species. An interspecies comparison of cCREs revealed both conserved and species-specific patterns of epigenetic evolution. Finally, we showed that comparisons of the epigenetic landscape between species can reveal elements with similar roles in regulation, even in the absence of genomic sequence alignment.

PtoMYB031, the R2R3 MYB transcription factor involved in secondary cell wall biosynthesis in poplar

Introduction

The biosynthesis of the secondary cell wall (SCW) is orchestrated by an intricate hierarchical transcriptional regulatory network. This network is initiated by first-layer master switches, SCW-NAC transcription factors, which in turn activate the second-layer master switches MYBs. These switches play a crucial role in regulating xylem specification and differentiation during SCW formation. However, the roles of most MYBs in woody plants are yet to be fully understood.

Methods

In this study, we identified and isolated the R2R3-MYB transcription factor, PtoMYB031, from Populus tomentosa. We explored its expression, mainly in xylem tissues, and its role as a transcriptional repressor in the nucleus. We used overexpression and RNA interference techniques in poplar, along with Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays, to analyze the regulatory effects of PtoMYB031.

Results

Overexpression of PtoMYB031 in poplar significantly reduced lignin, cellulose, and hemicellulose content, and inhibited vascular development in stems, resulting in decreased SCW thickness in xylem tissues. Gene expression analysis showed that structural genes involved in SCW biosynthesis were downregulated in PtoMYB031-OE lines. Conversely, RNA interference of PtoMYB031 increased these compounds. Additionally, PtoMYB031 was found to recruit the repressor PtoZAT11, forming a transcriptional inhibition complex.

Discussion

Our findings provide new insights into how PtoMYB031, through its interaction with PtoZAT11, forms a complex that can suppress the expression of key regulatory genes, PtoWND1A and PtoWND2B, in SCW biosynthesis. This study enhances our understanding of the transcriptional regulation involved in SCW formation in poplar, highlighting the significant role of PtoMYB031.

Evolution of lysine-specific demethylase 1 and REST corepressor gene families and their molecular interaction

Lysine-specific demethylase 1A (LSD1) binds to the REST corepressor (RCOR) protein family of corepressors to erase transcriptionally active marks on histones. Functional diversity in these complexes depends on the type of RCOR included, which modulates the catalytic activity of the complex. Here, we studied the duplicative history of the RCOR and LSD gene families and analyzed the evolution of their interaction. We found that RCOR genes are the product of the two rounds of whole-genome duplications that occurred early in vertebrate evolution. In contrast, the origin of the LSD genes traces back before to the divergence of animals and plants. Using bioinformatics tools, we show that the RCOR and LSD1 interaction precedes the RCOR repertoire expansion that occurred in the last common ancestor of jawed vertebrates. Overall, we trace LSD1-RCOR complex evolution and propose that animal non-model species offer advantages in addressing questions about the molecular biology of this epigenetic complex.

Microbiome and plant cell transformation trigger insect gall induction in cassava

Several specialised insects can manipulate normal plant development to induce a highly organised structure known as a gall, which represents one of the most complex interactions between insects and plants. Thus far, the mechanism for insect-induced plant galls has remained elusive. To study the induction mechanism of insect galls, we selected the gall induced by Iatrophobia brasiliensis (Diptera: Cecidomyiidae) in cassava (Euphorbiaceae: Manihot esculenta Crantz) as our model. PCR-based molecular markers and deep metagenomic sequencing data were employed to analyse the gall microbiome and to test the hypothesis that gall cells are genetically transformed by insect vectored bacteria. A shotgun sequencing discrimination approach was implemented to selectively discriminate between foreign DNA and the reference host plant genome. Several known candidate insertion sequences were identified, the most significant being DNA sequences found in bacterial genes related to the transcription regulatory factor CadR, cadmium-transporting ATPase encoded by the cadA gene, nitrate transport permease protein (nrtB gene), and arsenical pump ATPase (arsA gene). In addition, a DNA fragment associated with ubiquitin-like gene E2 was identified as a potential accessory genetic element involved in gall induction mechanism. Furthermore, our results suggest that the increased quality and rapid development of gall tissue are mostly driven by microbiome enrichment and the acquisition of critical endophytes. An initial gall-like structure was experimentally obtained in M. esculenta cultured tissues through inoculation assays using a Rhodococcus bacterial strain that originated from the inducing insect, which we related to the gall induction process. We provide evidence that the modification of the endophytic microbiome and the genetic transformation of plant cells in M. esculenta are two essential requirements for insect-induced gall formation. Based on these findings and having observed the same potential DNA marker in galls from other plant species (ubiquitin-like gene E2), we speculate that bacterially mediated genetic transformation of plant cells may represent a more widespread gall induction mechanism found in nature.

Comparative analysis of the myoglobin gene in whales and humans reveals evolutionary changes in regulatory elements and expression levels

Cetacea and other diving mammals have undergone numerous adaptations to their aquatic environment, among them high levels of the oxygen-carrying intracellular hemoprotein myoglobin in skeletal muscles. Hypotheses regarding the mechanisms leading to these high myoglobin levels often invoke the induction of gene expression by exercise, hypoxia, and other physiological gene regulatory pathways. Here we explore an alternative hypothesis: that cetacean myoglobin genes have evolved high levels of transcription driven by the intrinsic developmental mechanisms that drive muscle cell differentiation. We have used luciferase assays in differentiated C2C12 cells to test this hypothesis. Contrary to our hypothesis, we find that the myoglobin gene from the minke whale, Balaenoptera acutorostrata, shows a low level of expression, only about 8% that of humans. This low expression level is broadly shared among cetaceans and artiodactylans. Previous work on regulation of the human gene has identified a core muscle-specific enhancer comprised of two regions, the "AT element" and a C-rich sequence 5' of the AT element termed the "CCAC-box". Analysis of the minke whale gene supports the importance of the AT element, but the minke whale CCAC-box ortholog has little effect. Instead, critical positive input has been identified in a G-rich region 3' of the AT element. Also, a conserved E-box in exon 1 positively affects expression, despite having been assigned a repressive role in the human gene. Last, a novel region 5' of the core enhancer has been identified, which we hypothesize may function as a boundary element. These results illustrate regulatory flexibility during evolution. We discuss the possibility that low transcription levels are actually beneficial, and that evolution of the myoglobin protein toward enhanced stability is a critical factor in the accumulation of high myoglobin levels in adult cetacean muscle tissue.

How Many Sirtuin Genes Are Out There? Evolution of Sirtuin Genes in Vertebrates With a Description of a New Family Member

Studying the evolutionary history of gene families is a challenging and exciting task with a wide range of implications. In addition to exploring fundamental questions about the origin and evolution of genes, disentangling their evolution is also critical to those who do functional/structural studies to allow a deeper and more precise interpretation of their results in an evolutionary context. The sirtuin gene family is a group of genes that are involved in a variety of biological functions mostly related to aging. Their duplicative history is an open question, as well as the definition of the repertoire of sirtuin genes among vertebrates. Our results show a well-resolved phylogeny that represents an improvement in our understanding of the duplicative history of the sirtuin gene family. We identified a new sirtuin gene family member (SIRT3.2) that was apparently lost in the last common ancestor of amniotes but retained in all other groups of jawed vertebrates. According to our experimental analyses, elephant shark SIRT3.2 protein is located in mitochondria, the overexpression of which leads to an increase in cellular levels of ATP. Moreover, in vitro analysis demonstrated that it has deacetylase activity being modulated in a similar way to mammalian SIRT3. Our results indicate that there are at least eight sirtuin paralogs among vertebrates and that all of them can be traced back to the last common ancestor of the group that existed between 676 and 615 millions of years ago.

A survey of current methods to detect and genotype inversions

Polymorphic inversions are ubiquitous in humans, and they have been linked to both adaptation and disease. Following their discovery in Drosophila more than a century ago, inversions have proved to be more elusive than other structural variants. A wide variety of methods for the detection and genotyping of inversions have recently been developed: multiple techniques based on selective amplification by PCR, short- and long-read sequencing approaches, principal component analysis of small variant haplotypes, template strand sequencing, optical mapping, and various genome assembly methods. Many methods apply complex wet lab protocols or increasingly refined bioinformatic analyses. This review is an attempt to provide a practical summary and comparison of the methods that are in current use, with a focus on metrics such as the maximum size of segmental duplications at inversion breakpoints that each method can tolerate, the size range of inversions that they recover, their throughput, and whether the locations of putative inversions must be known beforehand. This article is protected by copyright. All rights reserved.

Multiple genome alignment in the telomere-to-telomere assembly era

With the arrival of telomere-to-telomere (T2T) assemblies of the human genome comes the computational challenge of efficiently and accurately constructing multiple genome alignments at an unprecedented scale. By identifying nucleotides across genomes which share a common ancestor, multiple genome alignments commonly serve as the bedrock for comparative genomics studies. In this review, we provide an overview of the algorithmic template that most multiple genome alignment methods follow. We also discuss prospective areas of improvement of multiple genome alignment for keeping up with continuously arriving high-quality T2T assembled genomes and for unlocking clinically-relevant insights.

Mutated KLF4(K409Q) in meningioma binds STRs and activates FGF3 gene expression

Krüppel-like factor 4 (KLF4) is a transcription factor that has been proven necessary for both induction and maintenance of pluripotency and self-renewal. Whole-genome sequencing defined a unique mutation in KLF4 (KLF4K409Q) in human meningiomas. However, the molecular mechanism of this tumor-specific KLF4 mutation is unknown. Using genome-wide high-throughput and focused quantitative transcriptional approaches in human cell lines, primary meningeal cells, and meningioma tumor tissue, we found that a change in the evolutionarily conserved DNA-binding domain of KLF4 alters its DNA recognition preference, resulting in a shift in downstream transcriptional activity. In the KLF4K409Q-specific targets, the normally silent fibroblast growth factor 3 (FGF3) is activated. We demonstrated a neomorphic function of KLF4K409Q in stimulating FGF3 transcription through binding to its promoter and in using short tandem repeats (STRs) located within the locus as enhancers.

A Journey Through Genetics to Biology

Although my engagement with human genetics emerged gradually, and sometimes serendipitously, it has held me spellbound for decades. Without my teachers, students, postdocs, colleagues, and collaborators, I would not be writing this review of my scientific adventures. Early gene and disease mapping was a satisfying puzzle-solving exercise, but building biological insight was my main goal. The project trajectory was hugely influenced by the evolutionarily conserved nature of the implicated genes and by the pace of progress in genetic technologies. The rich detail of clinical observations, particularly in eye disease, makes humans an excellent model, especially when complemented by the use of multiple other animal species for experimental validation. The contributions of collaborators and rivals also influenced our approach. We are very fortunate to work in this era of unprecedented progress in genetics and genomics. Expected final online publication date for the Annual Review of Genomics and Human Genetics, Volume 23 is October 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Epigenetic Regulation of Endothelial Cell Lineages During Zebrafish Development-New Insights From Technical Advances

Epigenetic regulation is integral in orchestrating the spatiotemporal regulation of gene expression which underlies tissue development. The emergence of new tools to assess genome-wide epigenetic modifications has enabled significant advances in the field of vascular biology in zebrafish. Zebrafish represents a powerful model to investigate the activity of cis-regulatory elements in vivo by combining technologies such as ATAC-seq, ChIP-seq and CUT&Tag with the generation of transgenic lines and live imaging to validate the activity of these regulatory elements. Recently, this approach led to the identification and characterization of key enhancers of important vascular genes, such as gata2a, notch1b and dll4. In this review we will discuss how the latest technologies in epigenetics are being used in the zebrafish to determine chromatin states and assess the function of the cis-regulatory sequences that shape the zebrafish vascular network.

Plant cytoplasmic ribosomal proteins: an update on classification, nomenclature, evolution and resources

Standardized naming systems are essential to integrate and unify distinct research fields, and to link multi-species data within and across kingdoms. We conducted a comprehensive survey of cytoplasmic ribosomal proteins (CRPs) in the dicot model Arabidopsis thaliana and the monocot model rice, noting that the standardized naming system has not been widely adopted in the plant community. We generated a database linking the old classical names to their updated and compliant names. We also explored the sequences, molecular evolution, and structural and functional characteristics of all plant CRP families, emphasizing evolutionarily conserved and plant-specific features through cross-kingdom comparisons. Unlike fungal CRP paralogs that were mainly created by whole-genome duplication (WGD) or retroposition under a concerted evolution mode, plant CRP genes evolved primarily through both WGD and tandem duplications in a rapid birth-and-death process. We also provide a web-based resource (http://www.plantcrp.cn/) with the aim of sharing the latest knowledge on plant CRPs and facilitating the continued development of a standardized framework across the entire community.

Complete chloroplast genome sequences of three aroideae species (Araceae): lights into selective pressure, marker development and phylogenetic relationships

Background

Colocasia gigantea, Caladium bicolor and Xanthosoma sagittifolium are three worldwide famous ornamental and/or vegetable plants in the Araceae family, these species in the subfamily Aroideae are phylogenetically perplexing due to shared interspecific morphological traits and variation.

Result

This study, for the first time ever, assembled and analyzed complete chloroplast genomes of C. gigantea, C. bicolor and X. sagittifolium with genome sizes of 165,906 bp, 153,149 bp and 165,169 bp in length, respectively. The genomes were composed of conserved quadripartite circular structures with a total of 131 annotated genes, including 8 rRNA, 37 tRNA and 86 protein-coding genes. A comparison within Aroideae showed seven protein-coding genes (accD, ndhF, ndhK, rbcL, rpoC1, rpoC2 and matK) linked to environmental adaptation. Phylogenetic analysis confirmed a close relationship of C. gigantea with C. esculenta and S. colocasiifolia, and the C. bicolor with X. sagittifolium. Furthermore, three DNA barcodes (atpH-atpI + psaC-ndhE, atpH-atpI + trnS-trnG, atpH-atpI + psaC-ndhE + trnS-trnG) harbored highly variable regions to distinguish species in Aroideae subfamily.

Conclusion

These results would be beneficial for species identification, phylogenetic relationship, genetic diversity, and potential of germplasm resources in Aroideae.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08400-3.

Perfect and imperfect views of ultraconserved sequences

Across the human genome, there are nearly 500 'ultraconserved' elements: regions of at least 200 contiguous nucleotides that are perfectly conserved in both the mouse and rat genomes. Remarkably, the majority of these sequences are non-coding, and many can function as enhancers that activate tissue-specific gene expression during embryonic development. From their first description more than 15 years ago, their extreme conservation has both fascinated and perplexed researchers in genomics and evolutionary biology. The intrigue around ultraconserved elements only grew with the observation that they are dispensable for viability. Here, we review recent progress towards understanding the general importance and the specific functions of ultraconserved sequences in mammalian development and human disease and discuss possible explanations for their extreme conservation.

Evolutionary rewiring of regulatory networks contributes to phenotypic differences between human and mouse orthologous genes

Mouse models have been engineered to reveal the biological mechanisms of human diseases based on an assumption. The assumption is that orthologous genes underlie conserved phenotypes across species. However, genetically modified mouse orthologs of human genes do not often recapitulate human disease phenotypes which might be due to the molecular evolution of phenotypic differences across species from the time of the last common ancestor. Here, we systematically investigated the evolutionary divergence of regulatory relationships between transcription factors (TFs) and target genes in functional modules, and found that the rewiring of gene regulatory networks (GRNs) contributes to the phenotypic discrepancies that occur between humans and mice. We confirmed that the rewired regulatory networks of orthologous genes contain a higher proportion of species-specific regulatory elements. Additionally, we verified that the divergence of target gene expression levels, which was triggered by network rewiring, could lead to phenotypic differences. Taken together, a careful consideration of evolutionary divergence in regulatory networks could be a novel strategy to understand the failure or success of mouse models to mimic human diseases. To help interpret mouse phenotypes in human disease studies, we provide quantitative comparisons of gene expression profiles on our website (http://sbi.postech.ac.kr/w/RN).

Isolation and molecular identification of native As-resistant bacteria: As(III) and As(V) removal capacity and possible mechanism of detoxification

The study of arsenic (As)-resistant microorganisms with high As removal capacity is fundamental for the development of economically sustainable technologies used for the treatment of water contaminated with metalloid. In the current study, four bacterial strains were isolated from As-contaminated water samples of the Xichu region, Mexico. Based on 16S rRNA gene sequencing and phylogenetic analysis of the isolated strains, Rhodococcus gordoniae, Microbacterium hydrocarbonoxydans, Exiguobacterium indicum, and Pseudomonas kribbensis were identified as potential As removal strains. R. gordoniae shows the highest growth capacity in both As(III) and As(V). R. gordoniae, M. hydrocarbonoxydans, and E. indicum removed approximately 81.6, 79.9, and 61.7% of As(III), as well as 77.2, 68.9, and 74.8% of As(V), respectively. P. kribbensis removed only about 80.2% of As(V). This study contributes to the possible detoxification mechanisms employed by these bacteria. Such insight could be crucial in the successful implementation of in situ bioremediation programs using these little-known bacteria.

Multiscale analysis of the randomization limits of the chromosomal gene organization between Lepidoptera and Diptera

How chromosome gene organization and gene content evolve among distantly related and structurally malleable genomes remains unresolved. This is particularly the case when considering different insect orders. We have compared the highly contiguous genome assemblies of the lepidopteran Danaus plexippus and the dipteran Drosophila melanogaster, which shared a common ancestor around 290 Ma. The gene content of 23 out of 30 D. plexippus chromosomes was significantly associated with one or two of the six chromosomal elements of the Drosophila genome, denoting common ancestry. Despite the phylogenetic distance, 9.6% of the 1-to-1 orthologues still reside within the same ancestral genome neighbourhood. Furthermore, the comparison D. plexippus-Bombyx mori indicated that the rates of chromosome repatterning are lower in Lepidoptera than in Diptera, although still within the same order of magnitude. Concordantly, 14 developmental gene clusters showed a higher tendency to retain full or partial clustering in D. plexippus, further supporting that the physical association between the SuperHox and NK clusters existed in the ancestral bilaterian. Our results illuminate the scope and limits of the evolution of the gene organization and content of the ancestral chromosomes to the Lepidoptera and Diptera while helping reconstruct portions of the genome in their most recent common ancestor.

The Organization of the Pig T-Cell Receptor γ (TRG) Locus Provides Insights into the Evolutionary Patterns of the TRG Genes across Cetartiodactyla

The domestic pig (Sus scrofa) is a species representative of the Suina, one of the four suborders within Cetartiodactyla. In this paper, we reported our analysis of the pig TRG locus in comparison with the loci of species representative of the Ruminantia, Tylopoda, and Cetacea suborders. The pig TRG genomic structure reiterates the peculiarity of the organization of Cetartiodactyla loci in TRGC “cassettes”, each containing the basic V-J-J-C unit. Eighteen genes arranged in four TRGC cassettes, form the pig TRG locus. All the functional TRG genes were expressed, and the TRGV genes preferentially rearrange with the TRGJ genes within their own cassette, which correlates the diversity of the γ-chain repertoire with the number of cassettes. Among them, the TRGC5, located at the 5′ end of the locus, is the only cassette that retains a marked homology with the corresponding TRGC cassettes of all the analyzed species. The preservation of the TRGC5 cassette for such a long evolutionary time presumes a highly specialized function of its genes, which could be essential for the survival of species. Therefore, the maintenance of this cassette in pigs confirms that it is the most evolutionarily ancient within Cetartiodactyla, and it has undergone a process of duplication to give rise to the other TRGC cassettes in the different artiodactyl species in a lineage-specific manner.

Genomic analysis revealed conserved acid tolerance mechanisms from native micro-organisms in fermented feed

AIMS

Fermented feed is an agricultural practice used in many regions of the world to improve the growth performance of farm animals. This study aimed to identify and evaluate the lactic acid bacteria and yeast involved in the production of fermented feed.

METHODS AND RESULTS

We isolated and described two micro-organisms from autochthonous microbiota origin present in a regional feed product, Lactobacillus paracasei IBR07 (Lacticaseibacillus paracasei) and Kazachstania unispora IBR014 (Saccharomyces unisporum). Genome sequence analyses were performed to characterize both micro-organisms. Potential pathways involved in the acid response, tolerance and persistence were predicted in both genomes. Although L. paracasei and K. unispora are considered safe for animal feed, we analysed the presence of virulence factors, antibiotic resistance and pathogenicity islands. Furthermore, the Galleria mellonella model was used to support the safety of both isolates.

CONCLUSIONS

We conclude that IBR07 and IBR014 strains are good candidates to be used as starter cultures for feed fermentation.

SIGNIFICANCE AND IMPACT OF THE STUDY

The data presented here will be helpful to explore other biotechnological aspects and constitute a starting point for further studies to establish the consumption benefit of fermented feed in farm animal production.

Olfactory expression of trace amine-associated receptors requires cooperative cis-acting enhancers

Olfactory sensory neurons express a large family of odorant receptors (ORs) and a small family of trace amine-associated receptors (TAARs). While both families are subject to so-called singular expression (expression of one allele of one gene), the mechanisms underlying TAAR gene choice remain obscure. Here, we report the identification of two conserved sequence elements in the mouse TAAR cluster (T-elements) that are required for TAAR gene expression. We observed that cell-type-specific expression of a TAAR-derived transgene required either T-element. Moreover, deleting either element reduced or abolished expression of a subset of TAAR genes, while deleting both elements abolished olfactory expression of all TAARs in cis with the mutation. The T-elements exhibit several features of known OR enhancers but also contain highly conserved, unique sequence motifs. Our data demonstrate that TAAR gene expression requires two cooperative cis-acting enhancers and suggest that ORs and TAARs share similar mechanisms of singular expression.

Characterization of the Mitochondrial Genome of a Wheat AL-Type Male Sterility Line and the Candidate CMS Gene

Heterosis utilization is very important in hybrid seed production. An AL-type cytoplasmic male sterile (CMS) line has been used in wheat-hybrid seed production, but its sterility mechanism has not been explored. In the present study, we sequenced and verified the candidate CMS gene in the AL-type sterile line (AL18A) and its maintainer line (AL18B). In the late uni-nucleate stage, the tapetum cells of AL18A showed delayed programmed cell death (PCD) and termination of microspore at the bi-nucleate stage. As compared to AL18B, the AL18A line produced 100% aborted pollens. The mitochondrial genomes of AL18A and AL18B were sequenced using the next generation sequencing such as Hiseq and PacBio. It was found that the mitochondrial genome of AL18A had 99% similarity with that of Triticum timopheevii, AL18B was identical to that of Triticum aestivum cv. Chinese Yumai. Based on transmembrane structure prediction, 12 orfs were selected as candidate CMS genes, including a previously suggested orf256. Only the lines harboring orf279 showed sterility in the transgenic Arabidopsis system, indicating that orf279 is the CMS gene in the AL-type wheat CMS lines. These results provide a theoretical basis and data support to further analyze the mechanism of AL-type cytoplasmic male sterility in wheat.

Haplotype Shuffling and Dimorphic Transposable Elements in the Human Extended Major Histocompatibility Complex Class II Region

The major histocompatibility complex (MHC) on chromosome 6p21 is one of the most single-nucleotide polymorphism (SNP)-dense regions of the human genome and a prime model for the study and understanding of conserved sequence polymorphisms and structural diversity of ancestral haplotypes/conserved extended haplotypes. This study aimed to follow up on a previous analysis of the MHC class I region by using the same set of 95 MHC haplotype sequences downloaded from a publicly available BioProject database at the National Center for Biotechnology Information to identify and characterize the polymorphic human leukocyte antigen (HLA)-class II genes, the MTCO3P1 pseudogene alleles, the indels of transposable elements as haplotypic lineage markers, and SNP-density crossover (XO) loci at haplotype junctions in DNA sequence alignments of different haplotypes across the extended class II region (∼1 Mb) from the telomeric PRRT1 gene in class III to the COL11A2 gene at the centromeric end of class II. We identified 42 haplotypic indels (20 Alu, 7 SVA, 13 LTR or MERs, and 2 indels composed of a mosaic of different transposable elements) linked to particular HLA-class II alleles. Comparative sequence analyses of 136 haplotype pairs revealed 98 unique XO sites between SNP-poor and SNP-rich genomic segments with considerable haplotype shuffling located in the proximity of putative recombination hotspots. The majority of XO sites occurred across various regions including in the vicinity of MTCO3P1 between HLA-DQB1 and HLA-DQB3, between HLA-DQB2 and HLA-DOB, between DOB and TAP2, and between HLA-DOA and HLA-DPA1, where most XOs were within a HERVK22 sequence. We also determined the genomic positions of the PRDM9-recombination suppression sequence motif ATCCATG/CATGGAT and the PRDM9 recombination activation partial binding motif CCTCCCCT/AGGGGAG in the class II region of the human reference genome (NC_ 000006) relative to published meiotic recombination positions. Both the recombination and anti-recombination PRDM9 binding motifs were widely distributed throughout the class II genomic regions with 50% or more found within repeat elements; the anti-recombination motifs were found mostly in L1 fragmented repeats. This study shows substantial haplotype shuffling between different polymorphic blocks and confirms the presence of numerous putative ancestral recombination sites across the class II region between various HLA class II genes.

Brassinosteroid-BZR1/2-WAT1 module determines the high level of auxin signalling in vascular cambium during wood formation

The tight regulation of local auxin homeostasis and signalling maxima in xylem precursor cells specifies the organising activity of the vascular cambium and consequently promotes xylem differentiation and wood formation. However, the molecular mechanisms underlying the local auxin signalling maxima in the vascular cambium are largely unknown. Here, we reveal that brassinosteroid (BR)-activated WALLS ARE THIN1 (WAT1) facilitates wood formation by enhancing local auxin signalling in the vascular cambium in Solanum lycopersicum. Growth defects and low auxin signalling readouts in the BR-deficient tomato cultivar, Micro-Tom, were associated with a novel recessive allele, Slwat1-copi, created by the insertion of a retrotransposon in the last exon of the SlWAT1 locus. Molecular and genetic studies by generating the gain-of-function and loss-of-function tomato mutants revealed that SlWAT1 is a critical regulator for fine tuning local auxin homeostasis and signalling outputs in vascular cambium to facilitate secondary growth. Finally, we discovered that BR-regulated SlBZR1/2 directly activated downstream auxin responses by SlWAT1 upregulation in xylem precursor cells to facilitate xylem differentiation and subsequent wood formation. Our data suggest that the BR-SlBZR1/2-WAT1 signalling network contributes to the high level of auxin signalling in the vascular cambium for secondary growth.

The T Cell Receptor (TRB) Locus in Tursiops truncatus: From Sequence to Structure of the Alpha/Beta Heterodimer in the Human/Dolphin Comparison

The bottlenose dolphin (Tursiops truncatus) belongs to the Cetartiodactyla and, similarly to other cetaceans, represents the most successful mammalian colonization of the aquatic environment. Here we report a genomic, evolutionary, and expression study of T. truncatus T cell receptor beta (TRB) genes. Although the organization of the dolphin TRB locus is similar to that of the other artiodactyl species, with three in tandem D-J-C clusters located at its 3' end, its uniqueness is given by the reduction of the total length due essentially to the absence of duplications and to the deletions that have drastically reduced the number of the germline TRBV genes. We have analyzed the relevant mature transcripts from two subjects. The simultaneous availability of rearranged T cell receptor α (TRA) and TRB cDNA from the peripheral blood of one of the two specimens, and the human/dolphin amino acids multi-sequence alignments, allowed us to calculate the most likely interactions at the protein interface between the alpha/beta heterodimer in complex with major histocompatibility class I (MH1) protein. Interacting amino acids located in the complementarity-determining region according to IMGT numbering (CDR-IMGT) of the dolphin variable V-alpha and beta domains were identified. According to comparative modelization, the atom pair contact sites analysis between the human MH1 grove (G) domains and the T cell receptor (TR) V domains confirms conservation of the structure of the dolphin TR/pMH.

The Genomic Organisation of the TRA/TRD Locus Validates the Peculiar Characteristics of Dromedary δ-Chain Expression

The role of γδ T cells in vertebrate immunity is still an unsolved puzzle. Species such as humans and mice display a low percentage of these T lymphocytes (i.e., “γδ low species”) with a restricted diversity of γδ T cell receptors (TR). Conversely, artiodactyl species (i.e., “γδ high species”) account for a high proportion of γδ T cells with large γ and δ chain repertoires. The genomic organisation of the TR γ (TRG) and δ (TRD) loci has been determined in sheep and cattle, noting that a wide number of germline genes that encode for γ and δ chains characterise their genomes. Taking advantage of the current improved version of the genome assembly, we have investigated the genomic structure and gene content of the dromedary TRD locus, which, as in the other mammalian species, is nested within the TR α (TRA) genes. The most remarkable finding was the identification of a very limited number of variable germline genes (TRDV) compared to sheep and cattle, which supports our previous expression analyses for which the somatic hypermutation mechanism is able to enlarge and diversify the primary repertoire of dromedary δ chains. Furthermore, the comparison between genomic and expressed sequences reveals that D genes, up to four incorporated in a transcript, greatly contribute to the increased diversity of the dromedary δ chain antigen binding-site.

EvoLSTM: context-dependent models of sequence evolution using a sequence-to-sequence LSTM

Motivation

Accurate probabilistic models of sequence evolution are essential for a wide variety of bioinformatics tasks, including sequence alignment and phylogenetic inference. The ability to realistically simulate sequence evolution is also at the core of many benchmarking strategies. Yet, mutational processes have complex context dependencies that remain poorly modeled and understood.

Results

We introduce EvoLSTM, a recurrent neural network-based evolution simulator that captures mutational context dependencies. EvoLSTM uses a sequence-to-sequence long short-term memory model trained to predict mutation probabilities at each position of a given sequence, taking into consideration the 14 flanking nucleotides. EvoLSTM can realistically simulate mammalian and plant DNA sequence evolution and reveals unexpectedly strong long-range context dependencies in mutation probabilities. EvoLSTM brings modern machine-learning approaches to bear on sequence evolution. It will serve as a useful tool to study and simulate complex mutational processes.

Availability and implementation

Code and dataset are available at https://github.com/DongjoonLim/EvoLSTM.

Supplementary information

Supplementary data are available at Bioinformatics online.

Genetic Mechanisms Underlying Cortical Evolution in Mammals

The remarkable sensory, motor, and cognitive abilities of mammals mainly depend on the neocortex. Thus, the emergence of the six-layered neocortex in reptilian ancestors of mammals constitutes a fundamental evolutionary landmark. The mammalian cortex is a columnar epithelium of densely packed cells organized in layers where neurons are generated mainly in the subventricular zone in successive waves throughout development. Newborn cells move away from their site of neurogenesis through radial or tangential migration to reach their specific destination closer to the pial surface of the same or different cortical area. Interestingly, the genetic programs underlying neocortical development diversified in different mammalian lineages. In this work, I will review several recent studies that characterized how distinct transcriptional programs relate to the development and functional organization of the neocortex across diverse mammalian lineages. In some primates such as the anthropoids, the neocortex became extremely large, especially in humans where it comprises around 80% of the brain. It has been hypothesized that the massive expansion of the cortical surface and elaboration of its connections in the human lineage, has enabled our unique cognitive capacities including abstract thinking, long-term planning, verbal language and elaborated tool making capabilities. I will also analyze the lineage-specific genetic changes that could have led to the modification of key neurodevelopmental events, including regulation of cell number, neuronal migration, and differentiation into specific phenotypes, in order to shed light on the evolutionary mechanisms underlying the diversity of mammalian brains including the human brain.

SNP-Density Crossover Maps of Polymorphic Transposable Elements and HLA Genes Within MHC Class I Haplotype Blocks and Junction

The genomic region (~4 Mb) of the human major histocompatibility complex (MHC) on chromosome 6p21 is a prime model for the study and understanding of conserved polymorphic sequences (CPSs) and structural diversity of ancestral haplotypes (AHs)/conserved extended haplotypes (CEHs). The aim of this study was to use a set of 95 MHC genomic sequences downloaded from a publicly available BioProject database at NCBI to identify and characterise polymorphic human leukocyte antigen (HLA) class I genes and pseudogenes, MICA and MICB, and retroelement indels as haplotypic lineage markers, and single-nucleotide polymorphism (SNP) crossover loci in DNA sequence alignments of different haplotypes across the Olfactory Receptor (OR) gene region (~1.2 Mb) and the MHC class I region (~1.8 Mb) from the GPX5 to the MICB gene. Our comparative sequence analyses confirmed the identity of 12 haplotypic retroelement markers and revealed that they partitioned the HLA-A/B/C haplotypes into distinct evolutionary lineages. Crossovers between SNP-poor and SNP-rich regions defined the sequence range of haplotype blocks, and many of these crossover junctions occurred within particular transposable elements, lncRNA, OR12D2, MUC21, MUC22, PSORS1A3, HLA-C, HLA-B, and MICA. In a comparison of more than 250 paired sequence alignments, at least 38 SNP-density crossover sites were mapped across various regions from GPX5 to MICB. In a homology comparison of 16 different haplotypes, seven CEH/AH (7.1, 8.1, 18.2, 51.x, 57.1, 62.x, and 62.1) had no detectable SNP-density crossover junctions and were SNP poor across the entire ~2.8 Mb of sequence alignments. Of the analyses between different recombinant haplotypes, more than half of them had SNP crossovers within 10 kb of LTR16B/ERV3-16A3_I, MLT1, Charlie, and/or THE1 sequences and were in close vicinity to structurally polymorphic Alu and SVA insertion sites. These studies demonstrate that (1) SNP-density crossovers are associated with putative ancestral recombination sites that are widely spread across the MHC class I genomic region from at least the telomeric OR12D2 gene to the centromeric MICB gene and (2) the genomic sequences of MHC homozygous cell lines are useful for analysing haplotype blocks, ancestral haplotypic landscapes and markers, CPSs, and SNP-density crossover junctions.

The barley pan-genome reveals the hidden legacy of mutation breeding

Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the 'pan-genome'1). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley (Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions2. Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley-comprising landraces, cultivars and a wild barley-that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

Characterization of the Complete Mitochondrial Genomes of Two Sibling Species of Parasitic Roundworms, Haemonchus contortus and Teladorsagia circumcincta

Haemonchus contortus and Teladorsagia circumcincta are among the two most pathogenic internal parasitic nematodes infecting small ruminants, such as sheep and goats, and are a global animal health issue. Accurate identification and delineation of Haemonchidae species is essential for development of diagnostic and control strategies with high resolution for Trichostrongyloidea infection in ruminants. Here, we describe in detail and compare the complete mitochondrial (mt) genomes of the New Zealand H. contortus and T. circumcincta field strains to improve our understanding of species- and strain-level evolution in these closely related roundworms. In the present study, we performed extensive comparative bioinformatics analyses on the recently sequenced complete mt genomes of the New Zealand H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP field strains. Amino acid sequences inferred from individual genes of each of the two mt genomes were compared, concatenated and subjected to phylogenetic analysis using Bayesian inference (BI), Maximum Likelihood (ML), and Maximum Parsimony (MP). The AT-rich mt genomes of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP are 14,001 bp (A+T content of 77.4%) and 14,081 bp (A+T content of 77.3%) in size, respectively. All 36 of the typical nematode mt genes are transcribed in the forward direction in both species and comprise of 12 protein-encoding genes (PCGs), 2 ribosomal RNA (rrn) genes, and 22 transfer RNA (trn) genes. The secondary structures for the 22 trn genes and two rrn genes differ between H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP, however the gene arrangements of both are consistent with other Trichostrongylidea sequenced to date. Comparative analyses of the complete mitochondrial nucleotide sequences, PCGs, A+T rich and non-coding repeat regions of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP further reinforces the high levels of diversity and gene flow observed among Trichostrongylidea, and supports their potential as ideal markers for strain-level identification from different hosts and geographical regions with high resolution for future studies. The complete mt genomes of H. contortus NZ_Hco_NP and T. circumcincta NZ_Teci_NP presented here provide useful novel markers for further studies of the meta-population connectivity and the genetic mechanisms driving evolution in nematode species.

Plastome comparison and evolution within the tribes of Plantaginaceae: Insights from an Asian gypsyweed

In spite of availability of several plastomes representing different tribes of Plantaginaceae, sparse attempts have been made to understand the plastome structure, evolution, and phylogenomics. In the present study, we have made an effort to understand the gene content and plastome evolution in the family Plantaginaceae using the newly generated plastome sequence of Veronica ovata subsp. kiusiana, a taxon native to SE Asia. In the first-ever attempt, plastomes of seven out of 10 tribes of Plantaginaceae have been compared to understand the evolution across the tribes of Plantaginaceae. The size of the plastome of V. ovata subsp. kiusiana is 152,249 bp, showing a typical quadripartite structure containing LSC, SSC, and two IRs with the sizes of 83,187, 17,704, and 25,679 respectively. The plastome comparison revealed the unique deletions in ycf2 and ndhF genes of members of different tribes, and also revealed high nucleotide variable hotspots. The study also revealed six highly variable genes and intergenic spacer viz. rps16, rps15-ycf1, ccsA-ndhD, ndhC-trnV, petN-psbM, and ycf1-trnN as potential DNA barcodes for the genus Veronica. The phylogenomic study revealed the sister relationship between V. ovata subsp. kiusiana and V. persica and also suggested the tentative placement of seven tribes in the family Plantaginaceae.

The expansion of the TRB and TRG genes in domestic goats (Capra hircus) is characteristic of the ruminant species

BACKGROUND

Goats (Capra hircus), one of the first domesticated species, are economically important for milk and meat production, and their broad geographical distribution reflects their successful adaptation to diverse environmental conditions. Despite the relevance of this species, the genetic research on the goat traits is limited compared to other domestic species. Thanks to the latest goat reference genomic sequence (ARS1), which is considered to be one of the most continuous assemblies in livestock, we deduced the genomic structure of the T cell receptor beta (TRB) and gamma (TRG) loci in this ruminant species.

RESULTS

Our analyses revealed that although the organization of the goat TRB locus is broadly similar to that of the other artiodactyl species, with three in-tandem D-J-C clusters located at the 3' end, a complex and extensive series of duplications have occurred in the V genes at the 5' end, leading to a marked expansion in the number of the TRBV genes. This phenomenon appears to be a feature of the ruminant lineage since similar gene expansions have also occurred in sheep and cattle. Likewise, the general organization of the goat TRG genes is typical of ruminant species studied so far, with two paralogous TRG loci, TRG1 and TRG2, located in two distinct and distant positions on the same chromosome as result of a split in the ancestral locus. Each TRG locus consists of reiterated V-J-J-C cassettes, with the goat TRG2 containing an additional cassette relative to the corresponding sheep and cattle loci.

CONCLUSIONS

Taken together, these findings demonstrate that strong evolutionary pressures in the ruminant lineage have selected for the development of enlarged sets of TRB and TRG genes that contribute to a diverse T cell receptor repertoire. However, differences observed among the goat, sheep and cattle TRB and TRG genes indicate that distinct evolutionary histories, with independent expansions and/or contractions, have also affected each ruminant species.

Tbx2a Modulates Switching of RH2 and LWS Opsin Gene Expression

Sensory systems are tuned by selection to maximize organismal fitness in particular environments. This tuning has implications for intraspecies communication, the maintenance of species boundaries, and speciation. Tuning of color vision largely depends on the sequence of the expressed opsin proteins. To improve tuning of visual sensitivities to shifts in habitat or foraging ecology over the course of development, many organisms change which opsins are expressed. Changes in this developmental sequence (heterochronic shifts) can create differences in visual sensitivity among closely related species. The genetic mechanisms by which these developmental shifts occur are poorly understood. Here, we use quantitative trait locus analyses, genome sequencing, and gene expression studies in African cichlid fishes to identify a role for the transcription factor Tbx2a in driving a switch between long wavelength sensitive (LWS) and Rhodopsin-like (RH2) opsin expression. We identify binding sites for Tbx2a in the LWS promoter and the highly conserved locus control region of RH2 which concurrently promote LWS expression while repressing RH2 expression. We also present evidence that a single change in Tbx2a regulatory sequence has led to a species difference in visual tuning, providing the first mechanistic model for the evolution of rapid switches in sensory tuning. This difference in visual tuning likely has important roles in evolution as it corresponds to differences in diet, microhabitat choice, and male nuptial coloration.

The novel coronavirus SARS-CoV-2: From a zoonotic infection to coronavirus disease 2019

The novel coronavirus (CoV), severe acute respiratory syndrome (SARS)‐CoV‐2 is an international public health emergency. Until now, the intermediate host and mechanisms of the interspecies jump of this virus are unknown. Phylogenetic analysis of all available bat CoV complete genomes was performed to analyze the relationships between bat CoV and SARS‐CoV‐2. To suggest a possible intermediate host, another phylogenetic reconstruction of CoV genomes obtained from animals that were hypothetically commercialized in the Chinese markets was also carried out. Moreover, mutation analysis was executed to suggest genomic regions that may have permitted the adaptation of SARS‐CoV‐2 to the human host. The phylogenetic analysis demonstrated that SARS‐CoV‐2 formed a cluster with the bat CoV isolate RaTG13. Possible CoV interspecies jumps among bat isolates were also observed. The phylogenetic tree reconstructed from CoV strains belonging to different animals demonstrated that SARS‐CoV‐2, bat RaTG13, and pangolin CoV genomes formed a monophyletic cluster, demonstrating that pangolins may be suggested as SARS‐CoV‐2 intermediate hosts. Three AA substitutions localized in the S1 portion of the S gene were observed, some of which have been correlated to structural modifications of the S protein which may facilitate SARS‐CoV‐2 tropism to human cells. Our analysis shows the tight relationship between SARS‐CoV‐2 and bat SARS‐like strains. It also hypothesizes that pangolins might have been possible intermediate hosts of the infection. Some of the observed AA substitutions in the S‐binding protein may serve as possible adaptation mutations in humans but more studies are needed to elucidate their function.

The intermediate host of SARS‐CoV‐2 is unknown.

Zoonotic genomic surveillance is important to understant SARS‐CoV‐2 evolutionary history.

Bat coronavirus internspecies jumps are important for the evolution of coronaviruses.

Specific aminoacid substitutions might have helped for the adaptation of SARS‐CoV‐2 in the human population.

Comprehensive genomic analysis of the dromedary T cell receptor gamma (TRG) locus and identification of a functional TRGC5 cassette

The emergent availability in public databases of more complete genome assemblies allows us to improve genomic data obtained by classical molecular cloning. The main goal of this study was to refine the genomic map of the dromedary TRG locus by integrating our previous genomic data with the analysis of recent genomic assemblies. We identified an additional TRGC cassette, defined as a V-J-C recombination unit, located at the 5' of the locus and made up of five TRGV genes followed by three TRGJ genes and one TRGC gene. Hence, the complete dromedary TRG locus spans about 105 Kb and consists of three in tandem TRGC cassettes delimited by AMPH and STARD3NL genes at the 5' and 3' end, respectively. An expression assay carried out on peripheral blood showed the functional competency for the dromedary TRGC5 cassette and confirmed the presence of the somatic hypermutation mechanism able to enlarge the repertoire diversity of the dromedary γδ T cells.

Organization of the Addax Major Histocompatibility Complex Provides Insights Into Ruminant Evolution

Ruminants are critical as prey in transferring solar energy fixed by plants into carnivorous species, yet the genetic signature of the driving forces leading to the evolutionary success of the huge number of ruminant species remains largely unknown. Here we report a complete DNA map of the major histocompatibility complex (MHC) of the addax (Addax nasomaculatus) genome by sequencing a total of 47 overlapping BAC clones previously mapped to cover the MHC region. The addax MHC is composed of 3,224,151 nucleotides, harboring a total of 150 coding genes, 50 tRNA genes, and 14 non-coding RNA genes. The organization of addax MHC was found to be highly conserved to those of sheep and cattle, highlighted by a large piece of chromosome inversion that divided the MHC class II into IIa and IIb subregions. It is now highly possible that all of the ruminant species in the family of Bovidae carry the same chromosome inversion in the MHC region, inherited from a common ancestor of ruminants. Phylogenetic analysis indicated that DY, a ruminant-specific gene located at the boundary of the inversion and highly expressed in dendritic cells, was possibly evolved from DQ, with an estimated divergence time ~140 million years ago. Homology modeling showed that the overall predicted structure of addax DY was similar to that of HLA-DQ2. However, the pocket properties of P1, P4, P6, and P9, which were critical for antigen binding in the addax DY, showed certain distinctive features. Structural analysis suggested that the populations of peptide antigens presented by addax DY and HLA-DQ2 were quite diverse, which in theory could serve to promote microbial regulation in the rumen by ruminant species, contributing to enhanced grass utilization ability. In summary, the results of our study helped to enhance our understanding of the MHC evolution and provided additional supportive evidence to our previous hypothesis that an ancient chromosome inversion in the MHC region of the last common ancestor of ruminants may have contributed to the evolutionary success of current ruminants on our planet.

Comparative Analysis of the TRB Locus in the Camelus Genus

T cells can be separated into two major subsets based on the heterodimer that forms their T cell receptors. αβ T cells have receptors consisting of α and β chains, while γδ T cells are composed of γ and δ chains. αβ T cells play an essential role within the adaptive immune responses against pathogens. The recent genomic characterization of the Camelus dromedarius T cell receptor β (TRB) locus has allowed us to infer the structure of this locus from the draft genome sequences of its wild and domestic Bactrian congeners, Camelus ferus and Camelus bactrianus. The general structural organization of the wild and domestic Bactrian TRB locus is similar to that of the dromedary, with a pool of TRBV genes positioned at the 5′ end of D-J-C clusters, followed by a single TRBV gene located at the 3′ end with an inverted transcriptional orientation. Despite the fragmented nature of the assemblies, comparative genomics reveals the existence of a perfect co-linearity between the three Old World camel TRB genomic sequences, which enables the transfer of information from one sequence to another and the filling of gaps in the genomic sequences. A virtual camelid TRB locus is hypothesized with the presence of 33 TRBV genes distributed in 26 subgroups. Likewise, in the artiodactyl species, three in-tandem D-J-C clusters, each composed of one TRBD gene, six or seven TRBJ genes, and one TRBC gene, are placed at the 3′ end of the locus. As reported in the ruminant species, a group of four functional TRY genes at the 5′ end and only one gene at the 3′ end, complete the camelid TRB locus. Although the gene content is similar, differences are observed in the TRBV functional repertoire, and genes that are functional in one species are pseudogenes in the other species. Hence, variations in the functional repertoire between dromedary, wild and domestic Bactrian camels, rather than differences in the gene content, may represent the molecular basis explaining the disparity in the TRB repertoire between the Camelus species. Finally, our data contribute to the knowledge about the evolutionary history of Old World camelids.

Evolution of Deeper Rooting 1-like homoeologs in wheat entails the C-terminus mutations as well as gain and loss of auxin response elements

Root growth angle (RGA) in response to gravity controlled by auxin is a pertinent target trait for obtainment of higher yield in cereals. But molecular basis of this root architecture trait remain obscure in wheat and barley. We selected four cultivars two each for wheat and barley to unveil the molecular genetic mechanism of Deeper Rooting 1-like gene which controls RGA in rice leading to higher yield under drought imposition. Morphological analyses revealed a deeper and vertically oriented root growth in “NARC 2009” variety of wheat than “Galaxy” and two other barley cultivars “Scarlet” and “ISR42-8”. Three new homoeologs designated as TaANDRO1-like, TaBNDRO1-like and TaDNDRO1-like corresponding to A, B and D genomes of wheat could be isolated from “NARC 2009”. Due to frameshift and intronization/exonization events the gene structures of these paralogs exhibit variations in size. DRO1-like genes with five distinct domains prevail in diverse plant phyla from mosses to angiosperms but in lower plants their differentiation from LAZY, NGR and TAC1 (root and shoot angle genes) is enigmatic. Instead of IGT as denominator motif of this family, a new C-terminus motif WxxTD in the V-domain is proposed as family specific motif. The EAR-like motif IVLEM at the C-terminus of the TaADRO1-like and TaDDRO1-like that diverged to KLHTLIPNK in TaBDRO1-like and HvDRO1-like is the hallmark of these proteins. Split-YFP and yeast two hybrid assays complemented the interaction of TaDRO1-like with TOPLESS—a repressor of auxin regulated root promoting genes in plants—through IVLEM/KLHTLIPNK motif. Quantitative RT-PCR revealed abundance of DRO1-like RNA in root tips and spikelets while transcript signals were barely detectable in shoot and leaf tissues. Interestingly, wheat exhibited stronger expression of TaBDRO1-like than barley (HvDRO1-like), but TaBDRO1-like was the least expressing among three paralogs. The underlying cause of this expression divergence seems to be the presence of AuxRE motif TGTCTC and core TGTC with a coupling AuxRE-like motif ATTTTCTT proximal to the transcriptional start site in TaBDRO1-like and HvDRO1-like promoters. This is evident from binding of ARF1 to TGTCTC and TGTC motifs of TaBDRO1-like as revealed by yeast one-hybrid assay. Thus, evolution of DRO1-like wheat homoeologs might incorporate the C-terminus mutations as well as gain and loss of AuxREs and other cis-regulatory elements during expression divergence. Since root architecture is an important target trait for wheat crop improvement, therefore DRO1-like genes have potential applications in plant breeding for enhancement of plant productivity by the use of modern genome editing approaches.

The Marine Mammal Class II Major Histocompatibility Complex Organization

Sirenians share with cetaceans and pinnipeds several convergent traits selected for the aquatic lifestyle. Living in water poses new challenges not only for locomotion and feeding but also for combating new pathogens, which may render the immune system one of the best tools aquatic mammals have for dealing with aquatic microbial threats. So far, only cetaceans have had their class II Major Histocompatibility Complex (MHC) organization characterized, despite the importance of MHC genes for adaptive immune responses. This study aims to characterize the organization of the marine mammal class II MHC using publicly available genomes. We located class II sequences in the genomes of one sirenian, four pinnipeds and eight cetaceans using NCBI-BLAST and reannotated the sequences using local BLAST search with exon and intron libraries. Scaffolds containing class II sequences were compared using dotplot analysis and introns were used for phylogenetic analysis. The manatee class II region shares overall synteny with other mammals, however most DR loci were translocated from the canonical location, past the extended class II region. Detailed analysis of the genomes of closely related taxa revealed that this presumed translocation is shared with all other living afrotherians. Other presumptive chromosome rearrangements in Afrotheria are the deletion of DQ loci in Afrosoricida and deletion of DP in E. telfairi. Pinnipeds share the main features of dog MHC: lack of a functional pair of DPA/DPB genes and inverted DRB locus between DQ and DO subregions. All cetaceans share the Cetartiodactyla inversion separating class II genes into two subregions: class IIa, with DR and DQ genes, and class IIb, with non-classic genes and a DRB pseudogene. These results point to three distinct and unheralded class II MHC structures in marine mammals: one canonical organization but lacking DP genes in pinnipeds; one bearing an inversion separating IIa and IIb subregions lacking DP genes found in cetaceans; and one with a translocation separating the most diverse class II gene from the MHC found in afrotherians and presumptive functional DR, DQ, and DP genes. Future functional research will reveal how these aquatic mammals cope with pathogen pressures with these divergent MHC organizations.