Splign: algorithms for computing spliced alignments with identification of paralogs

Structural genomic variation and behavioral interactions underpin a balanced sexual mimicry polymorphism

How phenotypic diversity originates and persists within populations are classic puzzles in evolutionary biology. While balanced polymorphisms segregate within many species, it remains rare for both the genetic basis and the selective forces to be known, leading to an incomplete understanding of many classes of traits under balancing selection. Here, we uncover the genetic architecture of a balanced sexual mimicry polymorphism and identify behavioral mechanisms that may be involved in its maintenance in the swordtail fish Xiphophorus birchmanni. We find that ∼40% of X. birchmanni males develop a "false gravid spot," a melanic pigmentation pattern that mimics the "pregnancy spot" associated with sexual maturity in female live-bearing fish. Using genome-wide association mapping, we detect a single intergenic region associated with variation in the false gravid spot phenotype, which is upstream of kitlga, a melanophore patterning gene. By performing long-read sequencing within and across populations, we identify complex structural rearrangements between alternate alleles at this locus. The false gravid spot haplotype drives increased allele-specific expression of kitlga, which provides a mechanistic explanation for the increased melanophore abundance that causes the spot. By studying social interactions in the laboratory and in nature, we find that males with the false gravid spot experience less aggression; however, they also receive increased attention from other males and are disdained by females. These behavioral interactions may contribute to the maintenance of this phenotypic polymorphism in natural populations. We speculate that structural variants affecting gene regulation may be an underappreciated driver of balanced polymorphisms across diverse species.

Phylogenomic analyses of all species of swordtail fishes (genus Xiphophorus) show that hybridization preceded speciation

Hybridization has been recognized to play important roles in evolution, however studies of the genetic consequence are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the genus Xiphophorus are proposed to have evolved with multiple ancient and ongoing hybridization events. They have served as an informative research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genomic resource including annotations for all described 26 Xiphophorus species and three undescribed taxa and resolve all uncertain phylogenetic relationships. We investigate the molecular evolution of genes related to cancers such as melanoma and for the genetic control of puberty timing, focusing on genes that are predicted to be involved in pre-and postzygotic isolation and thus affect hybridization. We discovered dramatic size-variation of some gene families. These persisted despite reticulate evolution, rapid speciation and short divergence time. Finally, we clarify the hybridization history in the entire genus settling disputed hybridization history of two Southern swordtails. Our comparative genomic analyses revealed hybridization ancestries that are manifested in the mosaic fused genomes and show that hybridization often preceded speciation.

Swordtail fish hybrids reveal that genome evolution is surprisingly predictable after initial hybridization

Over the past 2 decades, biologists have come to appreciate that hybridization, or genetic exchange between distinct lineages, is remarkably common-not just in particular lineages but in taxonomic groups across the tree of life. As a result, the genomes of many modern species harbor regions inherited from related species. This observation has raised fundamental questions about the degree to which the genomic outcomes of hybridization are repeatable and the degree to which natural selection drives such repeatability. However, a lack of appropriate systems to answer these questions has limited empirical progress in this area. Here, we leverage independently formed hybrid populations between the swordtail fish Xiphophorus birchmanni and X. cortezi to address this fundamental question. We find that local ancestry in one hybrid population is remarkably predictive of local ancestry in another, demographically independent hybrid population. Applying newly developed methods, we can attribute much of this repeatability to strong selection in the earliest generations after initial hybridization. We complement these analyses with time-series data that demonstrates that ancestry at regions under selection has remained stable over the past approximately 40 generations of evolution. Finally, we compare our results to the well-studied X. birchmanni × X. malinche hybrid populations and conclude that deeper evolutionary divergence has resulted in stronger selection and higher repeatability in patterns of local ancestry in hybrids between X. birchmanni and X. cortezi.

Conserved Signaling Pathways in the Ciona robusta Gut

The urochordate Ciona robusta exhibits numerous functional and morphogenetic traits that are shared with vertebrate models. While prior investigations have identified several analogies between the gastrointestinal tract (i.e., gut) of Ciona and mice, the molecular mechanisms responsible for these similarities remain poorly understood. This study seeks to address this knowledge gap by investigating the transcriptional landscape of the adult stage gut. Through comparative genomics analyses, we identified several evolutionarily conserved components of signaling pathways of pivotal importance for gut development (such as WNT, Notch, and TGFβ-BMP) and further evaluated their expression in three distinct sections of the gastrointestinal tract by RNA-seq. Despite the presence of lineage-specific gene gains, losses, and often unclear orthology relationships, the investigated pathways were characterized by well-conserved molecular machinery, with most components being expressed at significant levels throughout the entire intestinal tract of C. robusta. We also showed significant differences in the transcriptional landscape of the stomach and intestinal tract, which were much less pronounced between the proximal and distal portions of the intestine. This study confirms that C. robusta is a reliable model system for comparative studies, supporting the use of ascidians as a model to study gut physiology.

SimSpliceEvol2: alternative splicing-aware simulation of biological sequence evolution and transcript phylogenies

BACKGROUND

SimSpliceEvol is a tool for simulating the evolution of eukaryotic gene sequences that integrates exon-intron structure evolution as well as the evolution of the sets of transcripts produced from genes. It takes a guide gene tree as input and generates a gene sequence with its transcripts for each node of the tree, from the root to the leaves. However, the sets of transcripts simulated at different nodes of the guide gene tree lack evolutionary connections. Consequently, SimSpliceEvol is not suitable for evaluating methods for transcript phylogeny inference or gene phylogeny inference that rely on transcript conservation.

RESULTS

Here, we introduce SimSpliceEvol2, which, compared to the first version, incorporates an explicit model of transcript evolution for simulating alternative transcripts along the branches of a guide gene tree, as well as the transcript phylogenies inferred. We offer a comprehensive software with a graphical user interface and an updated version of the web server, ensuring easy and user-friendly access to the tool.

CONCLUSION

SimSpliceEvol2 generates synthetic datasets that are useful for evaluating methods and tools for spliced RNA sequence analysis, such as spliced alignment methods, methods for identifying conserved transcripts, and transcript phylogeny reconstruction methods. The web server is accessible at https://simspliceevol.cobius.usherbrooke.ca , where you can also download the standalone software. Comprehensive documentation for the software is available at the same address. For developers interested in the source code, which requires the installation of all prerequisites to run, it is provided at  https://github.com/UdeS-CoBIUS/SimSpliceEvol .

The complete sequence and comparative analysis of ape sex chromosomes

Apes possess two sex chromosomes-the male-specific Y chromosome and the X chromosome, which is present in both males and females. The Y chromosome is crucial for male reproduction, with deletions being linked to infertility1. The X chromosome is vital for reproduction and cognition2. Variation in mating patterns and brain function among apes suggests corresponding differences in their sex chromosomes. However, owing to their repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the methodology developed for the telomere-to-telomere (T2T) human genome, we produced gapless assemblies of the X and Y chromosomes for five great apes (bonobo (Pan paniscus), chimpanzee (Pan troglodytes), western lowland gorilla (Gorilla gorilla gorilla), Bornean orangutan (Pongo pygmaeus) and Sumatran orangutan (Pongo abelii)) and a lesser ape (the siamang gibbon (Symphalangus syndactylus)), and untangled the intricacies of their evolution. Compared with the X chromosomes, the ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements-owing to the accumulation of lineage-specific ampliconic regions, palindromes, transposable elements and satellites. Many Y chromosome genes expand in multi-copy families and some evolve under purifying selection. Thus, the Y chromosome exhibits dynamic evolution, whereas the X chromosome is more stable. Mapping short-read sequencing data to these assemblies revealed diversity and selection patterns on sex chromosomes of more than 100 individual great apes. These reference assemblies are expected to inform human evolution and conservation genetics of non-human apes, all of which are endangered species.

SEDA 2024 update: enhancing the SEquence DAtaset builder for seamless integration into automated data analysis pipelines

BACKGROUND

The initial version of SEDA assists life science researchers without programming skills with the preparation of DNA and protein sequence FASTA files for multiple bioinformatics applications. However, the initial version of SEDA lacks a command-line interface for more advanced users and does not allow the creation of automated analysis pipelines.

RESULTS

The present paper discusses the updates of the new SEDA release, including the addition of a complete command-line interface, new functionalities like gene annotation, a framework for automated pipelines, and improved integration in Linux environments.

CONCLUSION

SEDA is an open-source Java application and can be installed using the different distributions available ( https://www.sing-group.org/seda/download.html ) as well as through a Docker image ( https://hub.docker.com/r/pegi3s/seda ). It is released under a GPL-3.0 license, and its source code is publicly accessible on GitHub ( https://github.com/sing-group/seda ). The software version at the time of submission is archived at Zenodo (version v1.6.0, http://doi.org/10.5281/zenodo.10201605 ).

Annotations of four high-quality indigenous chicken genomes identify more than one thousand missing genes in subtelomeric regions and micro-chromosomes with high G/C contents

BACKGROUND

Although multiple chicken genomes have been assembled and annotated, the numbers of protein-coding genes in chicken genomes and their variation among breeds are still uncertain due to the low quality of these genome assemblies and limited resources used in their gene annotations. To fill these gaps, we recently assembled genomes of four indigenous chicken breeds with distinct traits at chromosome-level. In this study, we annotated genes in each of these assembled genomes using a combination of RNA-seq- and homology-based approaches.

RESULTS

We identified varying numbers (17,497-17,718) of protein-coding genes in the four indigenous chicken genomes, while recovering 51 of the 274 "missing" genes in birds in general, and 36 of the 174 "missing" genes in chickens in particular. Intriguingly, based on deeply sequenced RNA-seq data collected in multiple tissues in the four breeds, we found 571 ~ 627 protein-coding genes in each genome, which were missing in the annotations of the reference chicken genomes (GRCg6a and GRCg7b/w). After removing redundancy, we ended up with a total of 1,420 newly annotated genes (NAGs). The NAGs tend to be found in subtelomeric regions of macro-chromosomes (chr1 to chr5, plus chrZ) and middle chromosomes (chr6 to chr13, plus chrW), as well as in micro-chromosomes (chr14 to chr39) and unplaced contigs, where G/C contents are high. Moreover, the NAGs have elevated quadruplexes G frequencies, while both G/C contents and quadruplexes G frequencies in their surrounding regions are also high. The NAGs showed tissue-specific expression, and we were able to verify 39 (92.9%) of 42 randomly selected ones in various tissues of the four chicken breeds using RT-qPCR experiments. Most of the NAGs were also encoded in the reference chicken genomes, thus, these genomes might harbor more genes than previously thought.

CONCLUSION

The NAGs are widely distributed in wild, indigenous and commercial chickens, and they might play critical roles in chicken physiology. Counting these new genes, chicken genomes harbor more genes than originally thought.

Color vision evolution in egg-laying mammals: insights from visual photoreceptors and daily activities of Australian echidnas

Egg-laying mammals (monotremes) are considered "primitive" due to traits such as oviparity, cloaca, and incomplete homeothermy, all of which they share with reptiles. Two groups of monotremes, the terrestrial echidna (Tachyglossidae) and semiaquatic platypus (Ornithorhynchidae), have evolved highly divergent characters since their emergence in the Cenozoic era. These evolutionary differences, notably including distinct electrosensory and chemosensory systems, result from adaptations to species-specific habitat conditions. To date, very few studies have examined the visual adaptation of echidna and platypus. In the present study, we show that echidna and platypus have different light absorption spectra in their dichromatic visual sensory systems at the molecular level. We analyzed absorption spectra of monotreme color opsins, long-wavelength sensitive opsin (LWS) and short-wavelength sensitive opsin 2 (SWS2). The wavelength of maximum absorbance (λmax) in LWS was 570.2 in short-beaked echidna (Tachyglossus aculeatus) and 560.6 nm in platypus (Ornithorhynchus anatinus); in SWS2, λmax was 451.7 and 442.6 nm, respectively. Thus, the spectral range in echidna color vision is ~ 10 nm longer overall than in platypus. Natural selection analysis showed that the molecular evolution of monotreme color opsins is generally functionally conserved, suggesting that these taxa rely on species-specific color vision. In order to understand the usage of color vision in monotremes, we made 24-h behavioral observations of captive echidnas at warm temperatures and analyzed the resultant ethograms. Echidnas showed cathemeral activity and various behavioral repertoires such as feeding, traveling, digging, and self-grooming without light/dark environment selectivity. Halting (careful) behavior is more frequent in dark conditions, which suggests that echidnas may be more dependent on vision during the day and olfaction at night. Color vision functions have contributed to dynamic adaptations and dramatic ecological changes during the ~ 60 million years of divergent monotreme evolution. The ethogram of various day and night behaviors in captive echidnas also contributes information relevant to habitat conservation and animal welfare in this iconic species, which is locally endangered.

A premature stop codon in the CYP2C19 gene may explain the unexpected sensitivity of vultures to diclofenac toxicity

The unintended environmental exposure of vultures to diclofenac has resulted in the deaths of millions of old-world vultures on the Asian subcontinent. While toxicity has been since associated with a long half-life of elimination and zero order metabolism, the actual constraint in biotransformation is yet to be clarified. For this study we evaluated if the evident zero order metabolism could be due to defects in the CYP2C9/2C19 enzyme system. For this, using whole genome sequencing and de-novo transcriptome alignment, the vulture CYP2C19 open reading frame was identified through Splign analysis. The result sequence analysis revealed the presence of a premature stop codon on intron 7 of the identified open reading frame. Even if the stop codon was not present, amino acid residue analysis tended to suggest that the enzyme would be lower in activity than the equivalent human enzyme, with differences present at sites 105, 286 and 289. The defect was also conserved across the eight non-related vultures tested. From these results, we conclude that the sensitivity of the old-world vultures to diclofenac is due to the non-expression of a viable CYP2C19 enzyme system. This is not too dissimilar to the effects seen in certain people with a similar defective enzyme.

Insight into the aroma quality of 'Callista' cultivar of hop (Humulus lupulus L.): Impact of harvest timing, year, and location

Hops (Humulus lupulus L.) are essential ingredients in brewing, contributing to beer's flavor, aroma, and stability. This study pioneers an in-depth analysis of the 'Callista' cultivar, aiming to unravel how harvest timing, annual variations, and cultivation location synergistically influence its molecular profile, sensory perception, and biochemistry. Leveraging high-performance liquid chromatography and gas chromatography-mass spectrometry-olfactometry, we identified significant year-to-year and location-based fluctuations in bitter acids-the quintessential aroma constituents in hops. Our comprehensive aroma profiling discerned 55 volatile compounds, marking the first-ever sensory detection of 2-butanone in hops, with its presence showing remarkable interannual variability. This study showed significant differences among the three years tested, whereas hops were perceived "fruitier" and more "citrusy" in 2021, even though the bitter acid and aroma analysis showed that 2022 sticks out due to extremely high lupulone values up to 10% dry cone weight and 78% β-myrcene in the oil fraction compared to 60% and 45% in 2020 and 2021, respectively. Molecular analysis of key enzymes involved in hop aroma biosynthesis revealed no significant associations with location, but a strong diurnal pattern for all genes. The results indicated that especially the hot temperatures of 2022 may have induced significant changes of cone quality, while 2021 was more interesting from the sensory evaluations, which may justify the usage of viticultural terms such as "vintage" for hop marketing. These findings contribute to a better understanding of the factors influencing hop aroma and quality.

Genome evolution is surprisingly predictable after initial hybridization

Over the past two decades, evolutionary biologists have come to appreciate that hybridization, or genetic exchange between distinct lineages, is remarkably common - not just in particular lineages but in taxonomic groups across the tree of life. As a result, the genomes of many modern species harbor regions inherited from related species. This observation has raised fundamental questions about the degree to which the genomic outcomes of hybridization are repeatable and the degree to which natural selection drives such repeatability. However, a lack of appropriate systems to answer these questions has limited empirical progress in this area. Here, we leverage independently formed hybrid populations between the swordtail fish Xiphophorus birchmanni and X. cortezi to address this fundamental question. We find that local ancestry in one hybrid population is remarkably predictive of local ancestry in another, demographically independent hybrid population. Applying newly developed methods, we can attribute much of this repeatability to strong selection in the earliest generations after initial hybridization. We complement these analyses with time-series data that demonstrates that ancestry at regions under selection has remained stable over the past ~40 generations of evolution. Finally, we compare our results to the well-studied X. birchmanni×X. malinche hybrid populations and conclude that deeper evolutionary divergence has resulted in stronger selection and higher repeatability in patterns of local ancestry in hybrids between X. birchmanni and X. cortezi.

Chromosome-level genome assembly of chub mackerel (Scomber japonicus) from the Indo-Pacific Ocean

Chub mackerels (Scomber japonicus) are a migratory marine fish widely distributed in the Indo-Pacific Ocean. They are globally consumed for their high Omega-3 content, but their population is declining due to global warming. Here, we generated the first chromosome-level genome assembly of chub mackerel (fScoJap1) using the Vertebrate Genomes Project assembly pipeline with PacBio HiFi genomic sequencing and Arima Hi-C chromosome contact data. The final assembly is 828.68 Mb with 24 chromosomes, nearly all containing telomeric repeats at their ends. We annotated 31,656 genes and discovered that approximately 2.19% of the genome contained DNA transposon elements repressed within duplicated genes. Analyzing 5-methylcytosine (5mC) modifications using HiFi reads, we observed open/close chromatin patterns at gene promoters, including the FADS2 gene involved in Omega-3 production. This chromosome-level reference genome provides unprecedented opportunities for advancing our knowledge of chub mackerels in biology, industry, and conservation.

The Complete Sequence and Comparative Analysis of Ape Sex Chromosomes

Apes possess two sex chromosomes-the male-specific Y and the X shared by males and females. The Y chromosome is crucial for male reproduction, with deletions linked to infertility. The X chromosome carries genes vital for reproduction and cognition. Variation in mating patterns and brain function among great apes suggests corresponding differences in their sex chromosome structure and evolution. However, due to their highly repetitive nature and incomplete reference assemblies, ape sex chromosomes have been challenging to study. Here, using the state-of-the-art experimental and computational methods developed for the telomere-to-telomere (T2T) human genome, we produced gapless, complete assemblies of the X and Y chromosomes for five great apes (chimpanzee, bonobo, gorilla, Bornean and Sumatran orangutans) and a lesser ape, the siamang gibbon. These assemblies completely resolved ampliconic, palindromic, and satellite sequences, including the entire centromeres, allowing us to untangle the intricacies of ape sex chromosome evolution. We found that, compared to the X, ape Y chromosomes vary greatly in size and have low alignability and high levels of structural rearrangements. This divergence on the Y arises from the accumulation of lineage-specific ampliconic regions and palindromes (which are shared more broadly among species on the X) and from the abundance of transposable elements and satellites (which have a lower representation on the X). Our analysis of Y chromosome genes revealed lineage-specific expansions of multi-copy gene families and signatures of purifying selection. In summary, the Y exhibits dynamic evolution, while the X is more stable. Finally, mapping short-read sequencing data from >100 great ape individuals revealed the patterns of diversity and selection on their sex chromosomes, demonstrating the utility of these reference assemblies for studies of great ape evolution. These complete sex chromosome assemblies are expected to further inform conservation genetics of nonhuman apes, all of which are endangered species.

Phylloxera and Aphids Show Distinct Features of Genome Evolution Despite Similar Reproductive Modes

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared with the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with 1 sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared with their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

The late-evolving salmon and trout join the GnRH1 club

Although it is known that the whitefish, an ancient salmonid, expresses three distinct gonadotropin-releasing hormone (GnRH) forms in the brain, it has been thought that the later-evolving salmonids (salmon and trout) had only two types of GnRH: GnRH2 and GnRH3. We now provide evidence for the expression of GnRH1 in the gonads of Atlantic salmon by rapid amplification of cDNA ends, real-time quantitative PCR and immunohistochemistry. We examined six different salmonid genomes and found that each assembly has one gene that likely encodes a viable GnRH1 prepropeptide. In contrast to both functional GnRH2 and GnRH3 paralogs, the GnRH1 homeolog can no longer express the hormone. Furthermore, the viable salmonid GnRH1 mRNA is composed of only three exons, rather than the four exons that build the GnRH2 and GnRH3 mRNAs. Transcribed gnrh1 is broadly expressed (in 17/18 tissues examined), with relative abundance highest in the ovaries. Expression of the gnrh2 and gnrh3 mRNAs is more restricted, primarily to the brain, and not in the gonads. The GnRH1 proximal promoter presents composite binding elements that predict interactions with complexes that contain diverse cell fate and differentiation transcription factors. We provide immunological evidence for GnRH1 peptide in the nucleus of 1-year-old type A spermatogonia and cortical alveoli oocytes. GnRH1 peptide was not detected during other germ cell or reproductive stages. GnRH1 activity in the salmonid gonad may occur only during early stages of development and play a key role in a regulatory network that controls mitotic and/or meiotic processes within the germ cell.

G, Bashyal BM, Bhowmick PK, Vinod KK, Bollinedi H, Singh NK, Singh AK. Fine mapping of qBK1.2, a major QTL governing resistance to bakanae disease in rice

Bakanae disease caused by Fusarium fujikuroi is an emerging disease of rice causing losses in all rice-growing regions around the world. A BC2F2 population was developed by backcrossing the recurrent parent Pusa Basmati 1121 (PB1121) with the recombinant inbred line RIL28, which harbors a major quantitative trait locus (QTL) governing resistance to bakanae, qBK1.2. MassARRAY-based single-nucleotide polymorphism (SNP) assays targeting the genomic region of qBK1.2 helped in fine mapping the QTL to a region of 130 kb between the SNP markers rs3164311 and rs3295562 using 24 recombinants. In-silico mining of the fine-mapped region identified 11 putative candidate genes with functions related to defense. The expression analysis identified two significantly differentially expressed genes, that is, LOC_Os01g06750 and LOC_Os01g06870, between the susceptible genotype PB1121 and the resistant genotypes Pusa1342 and R-NIL4. Furthermore, the SNPs identified in LOC_Os01g06750 produced minor substitutions of amino acids with no major effect on the resistance-related functional motifs. However, LOC_Os01g06870 had 21 amino acid substitutions, which led to the creation of the leucine-rich repeat (LRR) domain in the resistant genotype Pusa1342, thereby making it a potential candidate underlying the major bakanae-resistant QTL qBK1.2. The markers used in the fine mapping program are of immense utility in marker-assisted breeding for bakanae resistance in rice.

Low mutation rate in epaulette sharks is consistent with a slow rate of evolution in sharks

Sharks occupy diverse ecological niches and play critical roles in marine ecosystems, often acting as apex predators. They are considered a slow-evolving lineage and have been suggested to exhibit exceptionally low cancer rates. These two features could be explained by a low nuclear mutation rate. Here, we provide a direct estimate of the nuclear mutation rate in the epaulette shark (Hemiscyllium ocellatum). We generate a high-quality reference genome, and resequence the whole genomes of parents and nine offspring to detect de novo mutations. Using stringent criteria, we estimate a mutation rate of 7×10-10 per base pair, per generation. This represents one of the lowest directly estimated mutation rates for any vertebrate clade, indicating that this basal vertebrate group is indeed a slowly evolving lineage whose ability to restore genetic diversity following a sustained population bottleneck may be hampered by a low mutation rate.

Genome sequencing provides insights into the evolution of gene families encoding plant cell wall-degrading enzymes in longhorned beetles

With more than 36,000 species, the longhorned beetles (family Cerambycidae) are a mega-diverse lineage of mostly xylophagous insects, all of which are represented by the sole sequenced genome of the Asian longhorned beetle (Anoplophora glabripennis; Lamiinae). Their successful radiation has been linked to their ability to degrade plant cell wall components using a range of so-called plant cell wall-degrading enzymes (PCWDEs). Our previous analysis of larval gut transcriptomes demonstrated that cerambycid beetles horizontally acquired genes encoding PCWDEs from various microbial donors; these genes evolved through multiple duplication events to form gene families. To gain further insights into the evolution of these gene families during the Cerambycidae radiation, we assembled draft genomes for four beetle species belonging to three subfamilies using long-read nanopore sequencing. All the PCWDE-encoding genes we annotated from the corresponding larval gut transcriptomes were present in these draft genomes. We confirmed that the newly discovered horizontally acquired glycoside hydrolase family 7 (GH7), subfamily 26 of GH43 (GH43_26), and GH53 (all of which are absent from the A. glabripennis genome) were indeed encoded by these beetles' genome. Most of the PCWDE-encoding genes of bacterial origin gained introns after their transfer into the beetle genome. Altogether, we show that draft genome assemblies generated from nanopore long-reads offer meaningful information to the study of the evolution of gene families in insects. We anticipate that our data will support studies aiming to better understand the biology of the Cerambycidae and other beetles in general.

Distinct monkeypox virus lineages co-circulating in humans before 2022

The 2022 global mpox outbreak raises questions about how this zoonotic disease established effective human-to-human transmission and its potential for further adaptation. The 2022 outbreak virus is related to an ongoing outbreak in Nigeria originally reported in 2017, but the evolutionary path linking the two remains unclear due to a lack of genomic data between 2018, when virus exportations from Nigeria were first recorded, and 2022, when the global mpox outbreak began. Here, 18 viral genomes obtained from patients across southern Nigeria in 2019-2020 reveal multiple lineages of monkeypox virus (MPXV) co-circulated in humans for several years before 2022, with progressive accumulation of mutations consistent with APOBEC3 activity over time. We identify Nigerian A.2 lineage isolates, confirming the lineage that has been multiply exported to North America independently of the 2022 outbreak originated in Nigeria, and that it has persisted by human-to-human transmission in Nigeria for more than 2 years before its latest exportation. Finally, we identify a lineage-defining APOBEC3-style mutation in all A.2 isolates that disrupts gene A46R, encoding a viral innate immune modulator. Collectively, our data demonstrate MPXV capacity for sustained diversification within humans, including mutations that may be consistent with established mechanisms of poxvirus adaptation.

The complete sequence of a human Y chromosome

The human Y chromosome has been notoriously difficult to sequence and assemble because of its complex repeat structure that includes long palindromes, tandem repeats and segmental duplications1-3. As a result, more than half of the Y chromosome is missing from the GRCh38 reference sequence and it remains the last human chromosome to be finished4,5. Here, the Telomere-to-Telomere (T2T) consortium presents the complete 62,460,029-base-pair sequence of a human Y chromosome from the HG002 genome (T2T-Y) that corrects multiple errors in GRCh38-Y and adds over 30 million base pairs of sequence to the reference, showing the complete ampliconic structures of gene families TSPY, DAZ and RBMY; 41 additional protein-coding genes, mostly from the TSPY family; and an alternating pattern of human satellite 1 and 3 blocks in the heterochromatic Yq12 region. We have combined T2T-Y with a previous assembly of the CHM13 genome4 and mapped available population variation, clinical variants and functional genomics data to produce a complete and comprehensive reference sequence for all 24 human chromosomes.

Phylloxera and aphids show distinct features of genome evolution despite similar reproductive modes

Genomes of aphids (family Aphididae) show several unusual evolutionary patterns. In particular, within the XO sex determination system of aphids, the X chromosome exhibits a lower rate of interchromosomal rearrangements, fewer highly expressed genes, and faster evolution at nonsynonymous sites compared to the autosomes. In contrast, other hemipteran lineages have similar rates of interchromosomal rearrangement for autosomes and X chromosomes. One possible explanation for these differences is the aphid's life cycle of cyclical parthenogenesis, where multiple asexual generations alternate with one sexual generation. If true, we should see similar features in the genomes of Phylloxeridae, an outgroup of aphids which also undergoes cyclical parthenogenesis. To investigate this, we generated a chromosome-level assembly for the grape phylloxera, an agriculturally important species of Phylloxeridae, and identified its single X chromosome. We then performed synteny analysis using the phylloxerid genome and 30 high-quality genomes of aphids and other hemipteran species. Unexpectedly, we found that the phylloxera does not share aphids' patterns of chromosome evolution. By estimating interchromosomal rearrangement rates on an absolute time scale, we found that rates are elevated for aphid autosomes compared to their X chromosomes, but this pattern does not extend to the phylloxera branch. Potentially, the conservation of X chromosome gene content is due to selection on XO males that appear in the sexual generation. We also examined gene duplication patterns across Hemiptera and uncovered horizontal gene transfer events contributing to phylloxera evolution.

Identification and Characterization of the Very-Low-Density Lipoprotein Receptor Gene from Branchiostoma belcheri: Insights into the Origin and Evolution of the Low-Density Lipoprotein Receptor Gene Family

Low-density lipoprotein receptors (LDLRs) are a class of cell-surface endocytosis receptors that are mainly involved in cholesterol homeostasis and cellular signal transduction. Very-low-density lipoprotein receptors (VLDLRs), which are members of the LDLR family, have been regarded as multi-function receptors that fulfill diverse physiological functions. However, no VLDLR gene has been identified in protochordates to date. As a representative protochordate species, amphioxi are the best available example of vertebrate ancestors. Identifying and characterizing the VLDLR gene in amphioxi has high importance for exploring the evolutionary process of the LDLR family. With this study, a new amphioxus VLDLR gene (designated AmphiVLDLR) was cloned and characterized using RACE-PCR. The 3217 nt transcript of the AmphiVLDLR had a 2577 nt ORF, and the deduced 858 amino acids were highly conserved within vertebrate VLDLRs according to their primary structure and three-dimensional structure, both of which contained five characteristic domains. In contrast to other vertebrate VLDLRs, which had a conserved genomic structure organization with 19 exons and 18 introns, the AmphiVLDLR had 13 exons and 12 introns. The results of a selective pressure analysis showed that the AmphiVLDLR had numerous positive selection sites. Furthermore, the tissue expression of AmphiVLDLR using RT-qPCR showed that AmphiVLDLR RNA expression levels were highest in the gills and muscles, moderate in the hepatic cecum and gonads, and lowest in the intestines. The results of the evolutionary analysis demonstrated that the AmphiVLDLR gene is a new member of the VLDLR family whose family members have experienced duplications and deletions over the evolutionary process. These results imply that the functions of LDLR family members have also undergone differentiation. In summary, we found a new VLDLR gene homolog (AmphiVLDLR) in amphioxi. Our results provide insight into the function and evolution of the LDLR gene family.

Decoding allelic diversity, transcript variants and transcriptional complexity of CENH3 gene in Brassica oleracea var. botrytis

Histone proteins play a critical role in the primary organization of nucleosomes, which is the fundamental unit of chromatin. Among the five types of the histones, histone H3 has multiple variants, and the number differs among the species. Amongst histone H3 variants, centromeric histone H3 (CENH3) is crucial for centromere identification and proper chromosomal segregation during cell division. In the present study, we have identified 17 putative histone H3 genes of Brassica oleracea. Furthermore, we have done a detailed characterization of the CENH3 gene of B. oleracea. We showed that a single CENH3 gene exhibits allelic diversity with at least two alleles and alternative splicing pattern. Also, we have identified a CENH3 gene-specific co-dominant cleaved amplified polymorphic sequence marker SNP34(A/C) to distinguish CENH3 alleles and follow their expression in leaf and flower tissues. The gene structure analysis of the CENH3 gene revealed the conserved 5'-CAGCAG-3' sequence at the intron 3-exon 4 junction in B. oleracea, which serves as an alternative splicing site with one-codon (alanine) addition/deletion. However, this one-codon alternative splicing feature is not conserved in the CENH3 genes of wild allied Brassica species. Our finding suggests that transcriptional complexity and alternative splicing might play a key role in the transcriptional regulation and function of the CENH3 gene in B. oleracea. Altogether, data generated from the present study can serve as a primary information resource and can be used to engineer CENH3 gene towards developing haploid inducer lines in B. oleracea.

Expression profiling and antibacterial analysis of cd36 in mandarin fish, Siniperca chuatsi

Cd36 is classified as a class B scavenger receptor and has also been identified as a pattern recognition receptor. In this study, we investigated the genomic structure and molecular characteristics of cd36 in mandarin fish (Siniperca chuatsi), examined its tissue distribution, and evaluated its antibacterial activity. Genomic structure analysis showed that Sccd36 consists of 12 exons and 11 introns. Sequencing analysis confirmed that the open reading frame of Sccd36 contains 1410 bp, encoding 469 amino acids. Sccd36 is deeply conserved with other vertebrates in terms of genomic structure, gene loci and molecular evolution, and the feature of two transmembrane was observed in ScCd36 through structural prediction. Sccd36 was constitutively expressed in all tissues tested, with the strongest expression in the intestine, followed by the heart and the kidney. Dramatic changes of Sccd36 mRNA were detected in mucosal tissues, including the intestine, gill and skin, when stimulated by the microbial ligands lipopolysaccharide and lipoteichoic acid. In addition, ScCd36 was identified as having strong binding ability to microbial ligands and antibacterial activity against the gram-negative bacteria Aeromonas hydrophila and the gram-positive bacteria Streptococcus lactis. Furthermore, we verified that the genetic ablation of cd36 impaired the resistance of fish to bacterial challenge by using zebrafish cd36 knockout line. In conclusion, our findings suggest that ScCd36 plays a crucial role in the innate immune response of mandarin fish against bacterial infections. This also sets the stage for further exploration into the antibacterial function of Cd36 in lower vertebrate species.

Genome-wide analysis of Catalase gene family reveal insights into abiotic stress response mechanism in Brassica juncea and B. rapa

Environmental stresses affect the yield and productivity of Brassica crops. Catalases are important antioxidant enzymes involved in reducing excess hydrogen peroxide produced by environmental stresses. In the present study, nine and seven CAT family members in two oilseed Brassica species (B. juncea and B. rapa) were identified with complete characterization based on gene and protein structure. Phylogenetic classification categorized CAT proteins into three classes and differentiated the monocot and dicot-specific CAT proteins. Further, the gene and protein characterizations revealed a high degree of conservation across the CAT family members. Differences were observed in the CAT-HEME binding affinity in CAT1, CAT2, and CAT3 isozymes, which could suggest their differential enzyme activities in different conditions. Furthermore, protein-protein interaction with other antioxidant proteins suggested their coordinated role in ROS scavenging mechanisms. Notably, the differential gene expression of BjuCATs and BraCATs and CAT enzyme activities suggested their crucial roles in major abiotic stresses faced by Brassica species. Promoter analysis in BjuCATs and BraCATs suggested the presence of abiotic-stress responsive cis-regulatory elements. Gene regulatory network analysis suggested miRNA and TF mediated stress response in BjuCATs and BraCATs. CAT family screening and characterization in Brassica sp. has established a basic ground for further functional validation in abiotic and heavy-metal stresses which can help in developing stress tolerant crops.

Transformer gene regulates feminization under two complementary sex determination loci in the ant, Vollenhovia emeryi

Organisms that reproduce sexually have evolved well-organized mechanisms to determine two sexes. Some hymenopterans (such as ants, bees, and wasps) have a complementary sex-determination system in which heterozygosity at one CSD locus induces female development, whereas hemi- or homozygosity at the locus induces male development. This system can generate a high cost of inbreeding, as individuals that are homozygous at the locus become sterile, diploid males. On the other hand, some hymenopterans have evolved a multi-locus, complementary, sex-determination system in which heterozygosity in at least one CSD locus induces female development. This system effectively reduces the proportion of sterile diploid males; however, how these multiple
primary signals based on CSD pass through a molecular cascade to regulate downstream genes has remained unclear. To clarify this matter, we used a backcross to investigate the molecular cascade in the ant, Vollenhovia emeryi, with two CSD loci. Here we show by gene disruption that transformer (tra) is necessary for proper feminization. Expression analysis of tra and doublesex (dsx) showed that heterozygosity in at least one of the two CSD loci is sufficient to promote female sex determination. Analysis of overexpression suggested that female-type Tra protein promotes splicing of tra pre-mRNA to female isoform by a positive-regulatory-feedback loop. Our data also showed that tra affects splicing of dsx. We conclude that two-loci sex determination system in V. emeryi evolved based on tra-dsx splicing cascade that is well conserved in other insect species. Finally, we suggest a cascade model to arrive at a binary determination of sex under multiple primary signals.

Nuclear import of IRF11 via the importin α/β pathway is essential for its antiviral activity

Interferon regulatory factor 11 (IRF11), an intriguing IRF member found only in fish species, has recently been shown to have antiviral properties that are dependent on its nuclear entry and DNA binding affinity. However, the mechanisms by which IRF11 enters the nucleus are unknown. In the present study, we found orthologs of IRF11 in lamprey and lancelet species by combining positional, phylogenetic and structural comparison data, showing that this gene has an ancient origin. The IRF11 gene (AjIRF11) from the Japanese eel, Anguilla japonica, was subsequently characterized, and it was found that AjIRF11 has antiviral activities against spring viremia of carp virus (SVCV), which are accomplished by regulating the production of type I IFN and IFN-stimulated genes. In addition to its known DNA binding residues in the α3 helix, two residues in Loop 1, His40 and Trp46, are also involved in DNA binding and activation of the IFN promoter. Using immunofluorescence microscopy and site-directed mutagenesis analysis, we confirmed that full nuclear localization of AjIRF11 requires the bipartite nuclear localization sequence (NLS) spanning residues 75 to 101, as well as the monopartite NLS situated between residues 119 and 122. Coimmunoprecipitation assays confirmed that AjIRF11 interacts with importin α via its NLSs and can also bind to importin β directly, implying that IRF11 can be imported to the nucleus by one or more transport receptors.

When does the female bias arise? Insights from the sex determination cascade of a flea beetle with a strongly skewed sex ratio

Reproduction-manipulating bacteria like Wolbachia can shift sex ratios in insects towards females, but skewed sex ratios may also arise from genetic conflicts. The flea beetle Altica lythri harbors three main mtDNA strains that are coupled to three different Wolbachia infections. Depending on the mtDNA types, the females produce either offspring with a balanced sex ratio or exclusively daughters. To obtain markers that can monitor when sex bias arises in the beetle's ontogeny, we elucidated the sex determination cascade of A. lythri. We established a RT-PCR method based on length variants of dsx (doublesex) transcripts to determine the sex of morphologically indistinguishable eggs and larvae. In females of one mtDNA type (HT1/HT1*) known to produce only daughters, male offspring were already missing at the egg stage while for females of another type (HT2), the dsx splice variants revealed a balanced sex ratio among eggs and larvae. Our data suggest that the sex determination cascade in A. lythri is initiated by maternally transmitted female-specific tra (transformer) mRNA as primary signal. This tra mRNA seems to be involved in a positive feedback loop that maintains the production of the female splice variant, as known for female offspring in Tribolium castaneum. The translation of the maternally transmitted female tra mRNA must be inhibited in male offspring, but the underlying primary genetic signal remains to be identified. We discuss which differences between the mtDNA types can influence sex determination and lead to the skewed sex ratio of HT1.

Pitfalls of using sequence databases for heterologous expression studies - a technical review

Synthesis of DNA fragments based on gene sequences that are available in public resources has become an efficient and affordable method that has gradually replaced traditional cloning efforts such as PCR cloning from cDNA. However, database entries based on genome sequencing results are prone to errors which can lead to false sequence information and, ultimately, errors in functional characterisation of proteins such as ion channels and transporters in heterologous expression systems. We have identified five common problems that repeatedly appear in public resources: (1) Not every gene has yet been annotated; (2) not all gene annotations are necessarily correct; (3) transcripts may contain automated corrections; (4) there are mismatches between gene, mRNA and protein sequences; and (5) splicing patterns often lack experimental validation. This technical review highlights and provides a strategy to bypass these issues in order to avoid critical mistakes that could impact future studies of any gene/protein of interest in heterologous expression systems.

In Silico Integrated Analysis of Genomic, Transcriptomic, and Proteomic Data Reveals QTL-Specific Genes for Bacterial Canker Resistance in Tomato (Solanum lycopersicum L.)

Bacterial canker of tomato, caused by Clavibacter michiganensis subsp. michiganensis (Cmm), is a devasting disease that leads to significant yield losses. Although QTLs originating from three wild species (Solanum arcanum, S. habrochaites, and S. pimpinellifolium) were identified, none of the QTLs was annotated for candidate gene identification. In the present study, a QTL-based physical map was constructed to reveal the meta-QTLs for Cmm resistance. As a result, seven major QTLs were mapped. Functional annotation of QTLs revealed 48 candidate genes. Additionally, experimentally validated Cmm resistance-related genes based on transcriptomic and proteomic studies were mapped in the genome and 25 genes were found to be located in the QTL regions. The present study is the first report to construct a physical map for Cmm resistance QTLs and identify QTL-specific candidate genes. The candidate genes identified in the present study are valuable targets for fine mapping and developing markers for marker-assisted selection in tomatoes for Cmm resistance breeding.

A highly diverse set of novel immunoglobulin-like transcript (NILT) genes in zebrafish indicates a wide range of functions with complex relationships to mammalian receptors

Multiple novel immunoglobulin-like transcripts (NILTs) have been identified from salmon, trout, and carp. NILTs typically encode activating or inhibitory transmembrane receptors with extracellular immunoglobulin (Ig) domains. Although predicted to provide immune recognition in ray-finned fish, we currently lack a definitive framework of NILT diversity, thereby limiting our predictions for their evolutionary origin and function. In order to better understand the diversity of NILTs and their possible roles in immune function, we identified five NILT loci in the Atlantic salmon (Salmo salar) genome, defined 86 NILT Ig domains within a 3-Mbp region of zebrafish (Danio rerio) chromosome 1, and described 41 NILT Ig domains as part of an alternative haplotype for this same genomic region. We then identified transcripts encoded by 43 different NILT genes which reflect an unprecedented diversity of Ig domain sequences and combinations for a family of non-recombining receptors within a single species. Zebrafish NILTs include a sole putative activating receptor but extensive inhibitory and secreted forms as well as membrane-bound forms with no known signaling motifs. These results reveal a higher level of genetic complexity, interindividual variation, and sequence diversity for NILTs than previously described, suggesting that this gene family likely plays multiple roles in host immunity.

Molecular characterization and expression analysis of immunoglobulins (IgM and IgT) heavy chains in black rockfish (Sebastes schlegelii) that response to bacterial challenge

Sebastes schlegelii is a kind of fish with great economic values. Recently, with the continuous expansion of aquaculture scale and the continuous improvement of aquaculture density, outbreak of various diseases has caused huge economic losses to its aquaculture industry. Study of fish immune system can help to understand the mechanism of immune response to external pathogens and can promote the development of immune prevention and control methods. Immunoglobulins (Igs) are complex glycoproteins that appear to be unique to the vertebrates that can recognize a wide variety of pathogens and recruit immune cells and molecules to destroy pathogens, which are generated by a series of rearrangement and somatic mutations. We therefore studied the immunoglobulin genes of S. schlegelii in view of their important roles in resisting to external pathogen infections. In this study, the immunoglobulin heavy chain genes (sIgM, mIgM, sIgT, and mIgT) of S. schlegelii were successfully identified and cloned. Phylogenetic analysis showed that the IgM and IgT genes of S. schlegelii were clustered together with homologous genes of other species, indicating that they were highly conserved during the evolutionary process. Collinearity analysis showed that the immunoglobulin genes and their adjacent genes were aligned with zebrafish, Atlantic salmon and tilapia, which further confirmed the conserved immunoglobulin gene of teleost. Expression analysis of healthy tissues showed that the expression levels of sIgM, sIgT and mIgT were the highest in the skin, while mIgM was the highest in spleen. After different bacterial infection, IgM and IgT were significantly expressed in skin and gill, which may be because skin and gill are the first line of defense against the infection pathogens. Subcellular localization showed that the mIgT protein was expressed in both the cell membrane and cytoplasm. Meanwhile, recombinant protein of mIgT was obtained in vitro, which laid a foundation for subsequent protein function studies. These results provide a theoretical basis for understanding the immunity role of immunoglobulin in S. schlegelii.

A truncated and catalytically inactive isoform of KDM5B histone demethylase accumulates in breast cancer cells and regulates H3K4 tri-methylation and gene expression

KDM5B histone demethylase is overexpressed in many cancers and plays an ambivalent role in oncogenesis, depending on the specific context. This ambivalence could be explained by the expression of KDM5B protein isoforms with diverse functional roles, which could be present at different levels in various cancer cell lines. We show here that one of these isoforms, namely KDM5B-NTT, accumulates in breast cancer cell lines due to remarkable protein stability relative to the canonical PLU-1 isoform, which shows a much faster turnover. This isoform is the truncated and catalytically inactive product of an mRNA with a transcription start site downstream of the PLU-1 isoform, and the consequent usage of an alternative ATG for translation initiation. It also differs from the PLU-1 transcript in the inclusion of an additional exon (exon-6), previously attributed to other putative isoforms. Overexpression of this isoform in MCF7 cells leads to an increase in bulk H3K4 methylation and induces derepression of a gene cluster, including the tumor suppressor Cav1 and several genes involved in the interferon-alpha and -gamma response. We discuss the relevance of this finding considering the hypothesis that KDM5B may possess regulatory roles independent of its catalytic activity.

Genomic Sequence of Canadian Chenopodium berlandieri: A North American Wild Relative of Quinoa

Chenopodium berlandieri (pitseed goosefoot) is a widespread native North American plant, which was cultivated and consumed by indigenous peoples prior to the arrival of European colonists. Chenopodium berlandieri is closely related to, and freely hybridizes with the domesticated South American food crop C. quinoa. As such it is a potential source of wild germplasm for breeding with C. quinoa, for improved quinoa production in North America. The C. berlandieri genome sequence could also be a useful source of information for improving quinoa adaptation. To this end, we first optimized barcode markers in two chloroplast genes, rbcL and matK. Together these markers can distinguish C. berlandieri from the morphologically similar Eurasian invasive C. album (lamb's quarters). Second, we performed whole genome sequencing and preliminary assembly of a C. berlandieri accession collected in Manitoba, Canada. Our assembly, while fragmented, is consistent with the expected allotetraploid structure containing diploid Chenopodium sub-genomes A and B. The genome of our accession is highly homozygous, with only one variant site per 3-4000 bases in non-repetitive sequences. This is consistent with predominant self-fertilization. As previously reported for the genome of a partly domesticated Mexican accession of C. berlandieri, our genome assembly is similar to that of C. quinoa. Somewhat unexpectedly, the genome of our accession had almost as many variant sites when compared to the Mexican C. berlandieri, as compared to C. quinoa. Despite the overall similarity of our genome sequence to that of C. quinoa, there are differences in genes known to be involved in the domestication or genetics of other food crops. In one example, our genome assembly appears to lack one functional copy of the SOS1 (salt overly sensitive 1) gene. SOS1 is involved in soil salinity tolerance, and by extension may be relevant to the adaptation of C. berlandieri to the wet climate of the Canadian region where it was collected. Our genome assembly will be a useful tool for the improved cultivation of quinoa in North America.

Protein-to-genome alignment with miniprot

MOTIVATION

Protein-to-genome alignment is critical to annotating genes in non-model organisms. While there are a few tools for this purpose, all of them were developed over 10 years ago and did not incorporate the latest advances in alignment algorithms. They are inefficient and could not keep up with the rapid production of new genomes and quickly growing protein databases.

RESULTS

Here, we describe miniprot, a new aligner for mapping protein sequences to a complete genome. Miniprot integrates recent techniques such as k-mer sketch and vectorized dynamic programming. It is tens of times faster than existing tools while achieving comparable accuracy on real data.

AVAILABILITY AND IMPLEMENTATION

https://github.com/lh3/miniport.

The profile of HSPA1A gene expression and its association with heat tolerance in crossbred cattle and the tropically adapted dwarf Vechur and Kasaragod

Certain livestock breeds are adapted to hot and humid environments, and these breeds have genetics that could be useful in a changing climate. The expression of several genes has been identified as a useful biomarker for heat stress. In this study, the responses to heat exposure of heat-tolerant Vechur and Kasaragod cattle found in Kerala state in India (also known as dwarf Bos taurus indicus) were compared to crossbred cattle (crosses of Bos t. taurus with Bos t. indicus). At various time points during heat exposure, rectal temperature and the expression of HSPA1A were determined, and the relationship between them was characterized. We characterized HSPA1A mRNA in Vechur cattle and performed molecular clock analysis. The expression of HSPA1A between the lineages and at different temperature humidity index (THI) was significant. There were significant differences between the expression profiles of HSPA1A in Kasaragod and crossbred (p < 0.01) and Vechur and crossbred (p < 0.01) cattle, but no significant difference in expression was observed between Vechur and Kasaragod cattle. The genetic distance between Vechur, B. grunniens, B. t. taurus, and B. t. indicus was 0.0233, 0.0059, and 0.007, respectively. The genetic distance between Vechur and the Indian dwarf breed Malnad Gidda was 0.0081. A molecular clock analysis revealed divergent adaptive evolution of Vechur cattle to B. t. taurus, with adaptations to the high temperatures and humidity that are prevalent in their breeding tract in Kerala, India. These results could also prove useful in selecting heat-tolerant animals using HSPA1A as a marker.

Discovery and functional characterization of neuropeptides in crinoid echinoderms

Neuropeptides are one of the largest and most diverse families of signaling molecules in animals and, accordingly, they regulate many physiological processes and behaviors. Genome and transcriptome sequencing has enabled the identification of genes encoding neuropeptide precursor proteins in species from a growing variety of taxa, including bilaterian and non-bilaterian animals. Of particular interest are deuterostome invertebrates such as the phylum Echinodermata, which occupies a phylogenetic position that has facilitated reconstruction of the evolution of neuropeptide signaling systems in Bilateria. However, our knowledge of neuropeptide signaling in echinoderms is largely based on bioinformatic and experimental analysis of eleutherozoans-Asterozoa (starfish and brittle stars) and Echinozoa (sea urchins and sea cucumbers). Little is known about neuropeptide signaling in crinoids (feather stars and sea lilies), which are a sister clade to the Eleutherozoa. Therefore, we have analyzed transcriptome/genome sequence data from three feather star species, Anneissia japonica, Antedon mediterranea, and Florometra serratissima, to produce the first comprehensive identification of neuropeptide precursors in crinoids. These include representatives of bilaterian neuropeptide precursor families and several predicted crinoid neuropeptide precursors. Using A. mediterranea as an experimental model, we have investigated the expression of selected neuropeptides in larvae (doliolaria), post-metamorphic pentacrinoids and adults, providing new insights into the cellular architecture of crinoid nervous systems. Thus, using mRNA in situ hybridization F-type SALMFamide precursor transcripts were revealed in a previously undescribed population of peptidergic cells located dorso-laterally in doliolaria. Furthermore, using immunohistochemistry a calcitonin-type neuropeptide was revealed in the aboral nerve center, circumoral nerve ring and oral tube feet in pentacrinoids and in the ectoneural and entoneural compartments of the nervous system in adults. Moreover, functional analysis of a vasopressin/oxytocin-type neuropeptide (crinotocin), which is expressed in the brachial nerve of the arms in A. mediterranea, revealed that this peptide causes a dose-dependent change in the mechanical behavior of arm preparations in vitro-the first reported biological action of a neuropeptide in a crinoid. In conclusion, our findings provide new perspectives on neuropeptide signaling in echinoderms and the foundations for further exploration of neuropeptide expression/function in crinoids as a sister clade to eleutherozoan echinoderms.

Semi-automated assembly of high-quality diploid human reference genomes

The current human reference genome, GRCh38, represents over 20 years of effort to generate a high-quality assembly, which has benefitted society1,2. However, it still has many gaps and errors, and does not represent a biological genome as it is a blend of multiple individuals3,4. Recently, a high-quality telomere-to-telomere reference, CHM13, was generated with the latest long-read technologies, but it was derived from a hydatidiform mole cell line with a nearly homozygous genome5. To address these limitations, the Human Pangenome Reference Consortium formed with the goal of creating high-quality, cost-effective, diploid genome assemblies for a pangenome reference that represents human genetic diversity6. Here, in our first scientific report, we determined which combination of current genome sequencing and assembly approaches yield the most complete and accurate diploid genome assembly with minimal manual curation. Approaches that used highly accurate long reads and parent-child data with graph-based haplotype phasing during assembly outperformed those that did not. Developing a combination of the top-performing methods, we generated our first high-quality diploid reference assembly, containing only approximately four gaps per chromosome on average, with most chromosomes within ±1% of the length of CHM13. Nearly 48% of protein-coding genes have non-synonymous amino acid changes between haplotypes, and centromeric regions showed the highest diversity. Our findings serve as a foundation for assembling near-complete diploid human genomes at scale for a pangenome reference to capture global genetic variation from single nucleotides to structural rearrangements.

Genome-wide identification and expression analysis of PIN gene family under phytohormone and abiotic stresses in Vitis Vinifera L

The auxin efflux transport proteins PIN-formed (PIN) has wide adaptability to hormone and abiotic stress, but the response mechanism of PINs in grape remains unclear. In this study, 12 members of VvPINs were identified and distributed on 8 chromosomes. The PIN genes of five species were divided into two subgroups, and the similarity of exons/introns and motifs of VvPIN genes were found in the same subgroup. Meanwhile, according to the examination of conserved motifs, the motif 3 included the conserved structure NPNTY. The promoter region of VvPIN gene family contained various cis-acting elements, which were related to light, abiotic stress, and hormones which are essential for growth and development. Additionally, VvPIN1, VvPIN9, and VvPIN11 proteins simultaneously interacted with the ARF, ABC, PINOID, GBF1, and VIT_08s0007g09010. The results of qRT-PCR revealed that the majority of the VvPINs were highly induced by NAA, GA3, ABA, MeJA, SA, NaCl, low-temperature (4 ℃), and PEG treatments, and the results were consistent with the prediction of the cis-acting elements. Moreover, the expression profile and quantitative real-time PCR (qRT-PCR) demonstrated that VvPIN genes were expressed in roots, stems, and leaves. The subcellular localization of VvPIN1 in Nicotiana benthamiana revealed that VvPIN1 was localized at the plasma membrane. Collectively, this study revealed that PIN genes could respond to various abiotic stresses, and provided a framework for regulating the expression of PIN genes to enhance the resistance of the grape.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-022-01239-8.

The genomic basis of evolutionary novelties in a leafhopper

Evolutionary innovations generate phenotypic and species diversity. Elucidating the genomic processes underlying such innovations is central to understanding biodiversity. In this study, we addressed the genomic basis of evolutionary novelties in the glassy-winged sharpshooter (Homalodisca vitripennis, GWSS), an agricultural pest. Prominent evolutionary innovations in leafhoppers include brochosomes, proteinaceous structures that are excreted and used to coat the body, and obligate symbiotic associations with two bacterial types that reside within cytoplasm of distinctive cell types. Using PacBio long-read sequencing and Dovetail Omni-C technology, we generated a chromosome-level genome assembly for the GWSS and then validated the assembly using flow cytometry and karyotyping. Additional transcriptomic and proteomic data were used to identify novel genes that underlie brochosome production. We found that brochosome-associated genes include novel gene families that have diversified through tandem duplications. We also identified the locations of genes involved in interactions with bacterial symbionts. Ancestors of the GWSS acquired bacterial genes through horizontal gene transfer (HGT), and these genes appear to contribute to symbiont support. Using a phylogenomics approach, we inferred HGT sources and timing. We found that some HGT events date to the common ancestor of the hemipteran suborder Auchenorrhyncha, representing some of the oldest known examples of HGT in animals. Overall, we show that evolutionary novelties in leafhoppers are generated by the combination of acquiring novel genes, produced both de novo and through tandem duplication, acquiring new symbiotic associations that enable use of novel diets and niches, and recruiting foreign genes to support symbionts and enhance herbivory.

Identification of the Gene Repertoire of the IMD Pathway and Expression of Antimicrobial Peptide Genes in Several Tissues and Hemolymph of the Cockroach Blattella germanica

Antimicrobial peptide (AMP) genes, triggered by Toll and IMD pathways, are essential components of the innate immune system in the German cockroach Blattella germanica. Besides their role in killing pathogenic bacteria, AMPs could be involved in controlling its symbiotic systems (endosymbiont and microbiota). We found that the IMD pathway was active in the adult female transcriptomes of six tissues (salivary glands, foregut, midgut, hindgut, Malpighian tubules and fat body) and hemolymph. Total expression of AMP genes was high in hemolymph and salivary glands and much lower in the other sample types. The expression of specific AMP genes was very heterogeneous among sample types. Two genes, defensin_g10 and drosomycin_g5, displayed relevant expression in the seven sample types, although higher in hemolymph. Other genes only displayed high expression in one tissue. Almost no expression of attacin-like and blattellicin genes was observed in any sample type, although some of them were among the genes with the highest expression in adult female whole bodies. The expression of AMP genes in salivary glands could help control pathogens ingested with food and even determine gut microbiota composition. The low expression levels in midgut and hindgut are probably related to the presence of beneficial microbiota. Furthermore, a reduction in the expression of AMP genes in fat body could be the way to prevent damage to the population of the endosymbiont Blattabacterium cuenoti within bacteriocytes.

Molecular Cloning and Expression Analysis of Thyrotropin-Releasing Hormone, and Its Possible Role in Gonadal Differentiation in Rice Field eel Monopterus albus

Rice field eel (Monopterus albus), a protogynous hermaphrodite fish, is a good model for the research of sex determination and gonadal differentiation in teleosts. In this study, we cloned the full-length cDNA sequence of trh, which encoded a predicted protein with 270 amino acids. Trh mainly expressed in the brain, followed by the ovary, testis, muscle and pituitary, and had low levels in other peripheral tissues. During natural sex reversal, trh mRNA expression levels exhibited a significant increase at the late intersexual stage in the hypothalamus. In the gonad, trh mRNA expression levels showed a trend of increase followed by decrease, and only increased significantly at the middle intersexual stage. No matter static incubation or intraperitoneal (IP) injection, TRH had no significant effect on trh and thyroid-stimulating hormone βsubunit (tshβ) mRNA expression levels, and serum T3, T4 and TRH release. After static incubation of ovarian fragments by TRH, the expression of gonadal soma derived factor (gsdf) was up-regulated significantly at both the doses of 10 and 100 nM. IP injection of TRH stimulated the expression of gsdf, and inhibited the expression of ovarian aromatase gene (cyp19a1a), accompanied by the increase of serum 11-KT levels. The results indicated that TRH may play a novel role in gonadal differentiation by the regulation of gonadal differentiation-related gene expression and sex steroid hormone secretion in rice field eel.

Merfin: improved variant filtering, assembly evaluation and polishing via k-mer validation

Variant calling has been widely used for genotyping and for improving the consensus accuracy of long-read assemblies. Variant calls are commonly hard-filtered with user-defined cutoffs. However, it is impossible to define a single set of optimal cutoffs, as the calls heavily depend on the quality of the reads, the variant caller of choice and the quality of the unpolished assembly. Here, we introduce Merfin, a k-mer based variant-filtering algorithm for improved accuracy in genotyping and genome assembly polishing. Merfin evaluates each variant based on the expected k-mer multiplicity in the reads, independently of the quality of the read alignment and variant caller's internal score. Merfin increased the precision of genotyped calls in several benchmarks, improved consensus accuracy and reduced frameshift errors when applied to human and nonhuman assemblies built from Pacific Biosciences HiFi and continuous long reads or Oxford Nanopore reads, including the first complete human genome. Moreover, we introduce assembly quality and completeness metrics that account for the expected genomic copy numbers.

The Auxin-Response Repressor IAA30 Is Down-Regulated in Reproductive Tissues of Apomictic Paspalum notatum

The capacity for apomixis in Paspalum notatum is controlled by a single-dominant genomic region, which shows strong synteny to a portion of rice chromosome 12 long arm. The locus LOC_Os12g40890, encoding the Auxin/Indole-3-Acetic Acid (Aux/IAA) family member OsIAA30, is located in this rice genomic segment. The objectives of this work were to identify transcripts coding for Aux/IAA proteins expressed in reproductive tissues of P. notatum, detect the OsIAA30 putative ortholog and analyze its temporal and spatial expression pattern in reproductive organs of sexual and apomictic plants. Thirty-three transcripts coding for AUX/IAA proteins were identified. Predicted protein alignment and phylogenetic analysis detected a highly similar sequence to OsIAA30 (named as PnIAA30) present in both sexual and apomictic samples. The expression assays of PnIAA30 showed a significant down-regulation in apomictic spikelets compared to sexual ones at the stages of anthesis and post-anthesis, representation levels negatively correlated with apospory expressivity and different localizations in sexual and apomictic ovules. Several PnIAA30 predicted interactors also appeared differentially regulated in the sexual and apomictic floral transcriptomes. Our results showed that an auxin-response repressor similar to OsIAA30 is down-regulated in apomictic spikelets of P. notatum and suggests a contrasting regulation of auxin signaling during sexual and asexual seed formation.

SEDA: A Desktop Tool Suite for FASTA Files Processing

SEDA (SEquence DAtaset builder) is a multiplatform desktop application for the manipulation of FASTA files containing DNA or protein sequences. The convenient graphical user interface gives access to a collection of simple (filtering, sorting, or file reformatting, among others) and advanced (BLAST searching, protein domain annotation, gene annotation, and sequence alignment) utilities not present in similar applications, which eases the work of life science researchers working with DNA and/or protein sequences, especially those who have no programming skills. This paper presents general guidelines on how to build efficient data handling protocols using SEDA, as well as practical examples on how to prepare high-quality datasets for single gene phylogenetic studies, the characterization of protein families, or phylogenomic studies. The user-friendliness of SEDA also relies on two important features: (i) the availability of easy-to-install distributable versions and installers of SEDA, including a Docker image for Linux, and (ii) the facility with which users can manage large datasets. SEDA is open-source, with GNU General Public License v3.0 license, and publicly available at GitHub (https://github.com/sing-group/seda). SEDA installers and documentation are available at https://www.sing-group.org/seda/.

Reconstructing the Origins of the Somatostatin and Allatostatin-C Signaling Systems Using the Accelerated Evolution of Biodiverse Cone Snail Toxins

Somatostatin and its related peptides (SSRPs) form an important family of hormones with diverse physiological roles. The ubiquitous presence of SSRPs in vertebrates and several invertebrate deuterostomes suggests an ancient origin of the SSRP signaling system. However, the existence of SSRP genes outside of deuterostomes has not been established, and the evolutionary history of this signaling system remains poorly understood. Our recent discovery of SSRP-like toxins (consomatins) in venomous marine cone snails (Conus) suggested the presence of a related signaling system in mollusks and potentially other protostomes. Here, we identify the molluscan SSRP-like signaling gene that gave rise to the consomatin family. Following recruitment into venom, consomatin genes experienced strong positive selection and repeated gene duplications resulting in the formation of a hyperdiverse family of venom peptides. Intriguingly, the largest number of consomatins was found in worm-hunting species (>400 sequences), indicating a homologous system in annelids, another large protostome phylum. Consistent with this, comprehensive sequence mining enabled the identification of SSRP-like sequences (and their corresponding orphan receptor) in annelids and several other protostome phyla. These results established the existence of SSRP-like peptides in many major branches of bilaterians and challenge the prevailing hypothesis that deuterostome SSRPs and protostome allatostatin-C are orthologous peptide families. Finally, having a large set of predator-prey SSRP sequences available, we show that although the cone snail's signaling SSRP-like genes are under purifying selection, the venom consomatin genes experience rapid directional selection to target receptors in a changing mix of prey.

Sex differences in the early life stages of the salmon louse Lepeophtheirus salmonis (Copepoda: Caligidae)

Salmon lice are ectoparasites on salmonids and feed on blood, mucus, and skin from their hosts. This causes high annual costs for treatment and control for the aquaculture industry. Salmon lice have a life cycle consisting of eight life stages. Sex determination by eye is only possible from the sixth stage onwards. A molecular sex determination has not been carried out so far, even though few individual sex-linked SNPs have been reported. In the present study, we used known sex-specific SNPs as a basis to sequence the complete sex-specific gene variants and used the sequence information to develop a sex determination assay. This assay could be used to determine the developmental speed of the two sexes already in the earliest life stages. Additionally, we sampled salmon lice in the nauplius II stage, determined the sex of each individual, pooled their RNA according to their sex, and used RNA sequencing to search for differences in gene expression and further sex-specific SNPs. We succeeded in developing a sex-determination assay that works on DNA or RNA from even the earliest larval stages of the salmon louse after hatching. At these early developmental stages, male salmon lice develop slightly quicker than females. We detected several previously unknown, sex-specific SNPs in our RNA-data seq, but only very few genes showed a differential expression between the sexes. Potential connections between SNPs, gene expression, and development are discussed.

Leptin Receptor Deficiency Results in Hyperphagia and Increased Fatty Acid Mobilization during Fasting in Rainbow Trout (Oncorhynchus mykiss)

Leptin is a pleiotropic hormone known for regulating appetite and metabolism. To characterize the role of leptin signaling in rainbow trout, we used CRISPR/Cas9 genome editing to disrupt the leptin receptor (LepR) genes, lepra1 and lepra2. We compared wildtype (WT) and mutant fish that were either fed to satiation or feed deprived for six weeks. The LepR mutants exhibited a hyperphagic phenotype, which led to heavier body weight, faster specific growth rate, increased viscero- and hepatosomatic indices, and greater condition factor. Muscle glycogen, plasma leptin, and leptin transcripts (lepa1) were also elevated in fed LepR mutant fish. Expression levels of several hypothalamic genes involved in feed regulation were analyzed (agrp, npy, orexin, cart-1, cart-2, pomc-a1, pomc-b). No differences were detected between fed WT and mutants except for pomc-b (proopiomelanocortin-b), where levels were 7.5-fold higher in LepR fed mutants, suggesting that pomc-b expression is regulated by leptin signaling. Fatty acid (FA) content did not statistically differ in muscle of fed mutant fish compared to WT. However, fasted mutants exhibited significantly lower muscle FA concentrations, suggesting that LepR mutants exhibit increased FA mobilization during fasting. These data demonstrate a key role for leptin signaling in lipid and energy mobilization in a teleost fish.