Genome-Wide Identification and Expression Analysis of ADK Gene Family Members in Cotton under Abiotic Stress.

Adenosine kinase (ADK) is a key enzyme widely distributed in plants, playing an important role in maintaining cellular energy homeostasis and regulating plant growth, development, and responses to environmental stresses. However, research on ADK genes in cotton (Gossypium hirsutum), an economically significant crop, has been limited. This study identified 92 ADK genes from four cotton species (G. arboreum, G. raimondii, G. hirsutum, and G. barbadense) using HMMER and Local BLASTP methods and classified them into six groups. Chromosomal localization revealed a random distribution of ADK genes in G. hirsutum, with 13 genes located on the At subgenome and 14 genes on the Dt subgenome. Gene structure analysis showed consistency in exon-intron organization within subgroups, while conserved motif analysis identified subgroup-specific motifs, indicating functional diversity. Synteny and collinearity mapping analysis revealed that the primary expansion mechanisms of the ADK gene family in cotton are polyploidy and segmental duplication. Cis-regulatory elements in GhADK promoters were classified into light response, hormone response, developmental regulation, and stress response. We also analyzed the expression patterns of GhADK genes under a low temperature (4 °C) and drought conditions. Most GhADK genes responded to cold stress with different expression patterns, indicating their roles in rapid response and long-term cold adaptation. Under drought stress, expression patterns varied, with some genes showing sustained high expression levels. The qRT-PCR validation of transcriptomic data confirmed the stress-induced expression patterns of selected GhADK genes. Functional analysis through the VIGS silencing of GhADK25 demonstrated its importance in cold and drought stress responses, with silencing resulting in poor growth under stress, highlighting its significance in stress tolerance. This study provides a basis for further understanding the evolutionary relationships and functions of the cotton ADK gene family.

The genomic basis of the plant island syndrome in Darwin's giant daisies.

The repeated, rapid and often pronounced patterns of evolutionary divergence observed in insular plants, or the 'plant island syndrome', include changes in leaf phenotypes, growth, as well as the acquisition of a perennial lifestyle. Here, we sequence and describe the genome of the critically endangered, Galápagos-endemic species Scalesia atractyloides Arnot., obtaining a chromosome-resolved, 3.2-Gbp assembly containing 43,093 candidate gene models. Using a combination of fossil transposable elements, k-mer spectra analyses and orthologue assignment, we identify the two ancestral genomes, and date their divergence and the polyploidization event, concluding that the ancestor of all extant Scalesia species was an allotetraploid. There are a comparable number of genes and transposable elements across the two subgenomes, and while their synteny has been mostly conserved, we find multiple inversions that may have facilitated adaptation. We identify clear signatures of selection across genes associated with vascular development, growth, adaptation to salinity and flowering time, thus finding compelling evidence for a genomic basis of the island syndrome in one of Darwin's giant daisies.

Genome-Wide Identification and Characterization of the RCI2 Gene Family in Allotetraploid Brassica napus Compared with Its Diploid Progenitors.

Brassica napus and its diploid progenitors (B. rapa and B. oleracea) are suitable for studying the problems associated with polyploidization. As an important anti-stress protein, RCI2 proteins widely exist in various tissues of plants, and are crucial to plant growth, development, and stress response. In this study, the RCI2 gene family was comprehensively identified and analyzed, and 9, 9, and 24 RCI2 genes were identified in B. rapa, B. oleracea, and B. napus, respectively. Phylogenetic analysis showed that all of the identified RCI2 genes were divided into two groups, and further divided into three subgroups. Ka/Ks analysis showed that most of the identified RCI2 genes underwent a purifying selection after the duplication events. Moreover, gene structure analysis showed that the structure of RCI2 genes is largely conserved during polyploidization. The promoters of the RCI2 genes in B. napus contained more cis-acting elements, which were mainly involved in plant development and growth, plant hormone response, and stress responses. Thus, B. napus might have potential advantages in some biological aspects. In addition, the changes of RCI2 genes during polyploidization were also discussed from the aspects of gene number, gene structure, gene relative location, and gene expression, which can provide reference for future polyploidization analysis.

Genome-Wide Identification, Characterization and Expression Analysis of the JAZ Gene Family in Resistance to Gray Leaf Spots in Tomato.

The plant disease resistance system involves a very complex regulatory network in which jasmonates play a key role in response to external biotic or abiotic stresses. As inhibitors of the jasmonic acid (JA) signaling pathway, JASMONATE ZIM domain (JAZ) proteins have been identified in many plant species, and their functions are gradually being clarified. In this study, 26 JAZ genes were identified in tomato. The physical and chemical properties, predicted subcellular localization, gene structure, cis-acting elements, and interspecies collinearity of 26 SlJAZ genes were subsequently analyzed. RNA-seq data combined with qRT-PCR analysis data showed that the expression of most SlJAZ genes were induced in response to Stemphylium lycopersici, methyl jasmonate (MeJA) and salicylic acid (SA). Tobacco rattle virus RNA2-based VIGS vector (TRV2)-SlJAZ25 plants were more resistant to tomato gray leaf spots than TRV2-00 plants. Therefore, we speculated that SlJAZ25 played a negative regulatory role in tomato resistance to gray leaf spots. Based on combining the results of previous studies and those of our experiments, we speculated that SlJAZ25 might be closely related to JA and SA hormone regulation. SlJAZ25 interacted with SlJAR1, SlCOI1, SlMYC2, and other resistance-related genes to form a regulatory network, and these genes played an important role in the regulation of tomato gray leaf spots. The subcellular localization results showed that the SlJAZ25 gene was located in the nucleus. Overall, this study is the first to identify and analyze JAZ family genes in tomato via bioinformatics approaches, clarifying the regulatory role of SlJAZ25 genes in tomato resistance to gray leaf spots and providing new ideas for improving plant disease resistance.

Genome-Wide Identification of the Nramp Gene Family in Spirodela polyrhiza and Expression Analysis under Cadmium Stress.

Natural resistance-associated macrophage proteins (Nramps) are specific metal transporters in plants with different functions among various species. The evolutionary and functional information of the Nramp gene family in Spirodela polyrhiza has not been previously reported in detail. To identify the Nramp genes in S. polyrhiza, we performed genome-wide identification, characterization, classification, and cis-elements analysis among 22 species with 138 amino acid sequences. We also conducted chromosomal localization and analyzed the synteny relationship, promoter, subcellular localization, and expression patterns in S. polyrhiza. ß-Glucuronidase staining indicated that SpNramp1 and SpNramp3 mainly accumulated in the root and joint between mother and daughter frond. Moreover, SpNramp1 was also widely displayed in the frond. SpNramp2 was intensively distributed in the root and frond. Quantitative real-time PCR results proved that the SpNramp gene expression level was influenced by Cd stress, especially in response to Fe or Mn deficiency. The study provides detailed information on the SpNramp gene family and their distribution and expression, laying a beneficial foundation for functional research.

Genome-Wide Identification and Expression Profiling of the WOX Gene Family in Citrus sinensis and Functional Analysis of a CsWUS Member.

WUSCHEL-related homeobox (WOX) transcription factors (TFs) are well known for their role in plant development but are rarely studied in citrus. In this study, we identified 11 putative genes from the sweet orange genome and divided the citrus WOX genes into three clades (modern/WUSCHEL(WUS), intermediate, and ancient). Subsequently, we performed syntenic relationship, intron-exon organization, motif composition, and cis-element analysis. Co-expression analysis based on RNA-seq and tissue-specific expression patterns revealed that CsWOX gene expression has multiple intrinsic functions. CsWUS homolog of AtWUS functions as a transcriptional activator and binds to specific DNA. Overexpression of CsWUS in tobacco revealed dramatic phenotypic changes, including malformed leaves and reduced gynoecia with no seed development. Silencing of CsWUS in lemon using the virus-induced gene silencing (VIGS) system implied the involvement of CsWUS in cells of the plant stem. In addition, CsWUS was found to interact with CsCYCD3, an ortholog in Arabidopsis (AtCYCD3,1). Yeast one-hybrid screening and dual luciferase activity revealed that two TFs (CsRAP2.12 and CsHB22) bind to the promoter of CsWUS and regulate its expression. Altogether, these results extend our knowledge of the WOX gene family along with CsWUS function and provide valuable findings for future study on development regulation and comprehensive data of WOX members in citrus.

Impacts of allopolyploidization and structural variation on intraspecific diversification in Brassica rapa.

BACKGROUND: Despite the prevalence and recurrence of polyploidization in the speciation of flowering plants, its impacts on crop intraspecific genome diversification are largely unknown. Brassica rapa is a mesopolyploid species that is domesticated into many subspecies with distinctive morphotypes. RESULTS: Herein, we report the consequences of the whole-genome triplication (WGT) on intraspecific diversification using a pan-genome analysis of 16 de novo assembled and two reported genomes. Among the genes that derive from WGT, 13.42% of polyploidy-derived genes accumulate more transposable elements and non-synonymous mutations than other genes during individual genome evolution. We denote such genes as being "flexible." We construct the Brassica rapa ancestral genome and observe the continuing influence of the dominant subgenome on intraspecific diversification in B. rapa. The gene flexibility is biased to the more fractionated subgenomes (MFs), in contrast to the more intact gene content of the dominant LF (least fractionated) subgenome. Furthermore, polyploidy-derived flexible syntenic genes are implicated in the response to stimulus and the phytohormone auxin; this may reflect adaptation to the environment. Using an integrated graph-based genome, we investigate the structural variation (SV) landscapes in 524 B. rapa genomes. We observe that SVs track morphotype domestication. Four out of 266 candidate genes for Chinese cabbage domestication are speculated to be involved in the leafy head formation. CONCLUSIONS: This pan-genome uncovers the possible contributions of allopolyploidization on intraspecific diversification and the possible and underexplored role of SVs in favorable trait domestication. Collectively, our work serves as a rich resource for genome-based B. rapa improvement.

Subtelomeric assembly of a multi-gene pathway for antimicrobial defense compounds in cereals.

Non-random gene organization in eukaryotes plays a significant role in genome evolution. Here, we investigate the origin of a biosynthetic gene cluster for production of defence compounds in oat-the avenacin cluster. We elucidate the structure and organisation of this 12-gene cluster, characterise the last two missing pathway steps, and reconstitute the entire pathway in tobacco by transient expression. We show that the cluster has formed de novo since the divergence of oats in a subtelomeric region of the genome that lacks homology with other grasses, and that gene order is approximately colinear with the biosynthetic pathway. We speculate that the positioning of the late pathway genes furthest away from the telomere may mitigate against a 'self-poisoning' scenario in which toxic intermediates accumulate as a result of telomeric gene deletions. Our investigations reveal a striking example of adaptive evolution underpinned by remarkable genome plasticity.

Comparison of Brassica Genomes reveals asymmetrical gene retention between functional groups of genes in recurrent polyploidizations.

KEY MESSAGE: We provided a study on homeologous gene evolution of homeologous genes by comparing Brassica genomes. Polyploidy has played fundamental roles during the evolution of plants. Following polyploidization, many duplicated genes are diversified or lost in a process termed diploidization. Understanding the retention and diversification of homeologs after polyploidization will help elucidate the process of diploidization. Here, we investigated the evolution of homeologous genes in Brassica genomes and observed similarly asymmetrical gene retention among different functional groups and consistent retention after recurrent polyploidizations. In the comparative analysis of Brassica diploid genomes, we found that preferentially retained genes show different patterns on sequence and expression divergence: genes with the function of 'biosynthetic process' and 'transport' were under much stronger purifying selection, while transcriptional regulatory genes diverged much faster than other genes. Duplicate pairs of the former two functional groups show conserved high expression patterns, while most of transcriptional regulatory genes are simultaneously lowly expressed. Furthermore, homeologs in diploids and allotetraploids showed similar loss and retention patterns: duplicates in progenitor genomes were more likely to be retained and accumulated fewer substitutions. However, transcriptional regulation is also enriched in the genes that do not have any non-synonymous mutations in the Brassica allotetraploids, indicating that some of these genes were under strong purifying selection. Overall, our study provided insight into the evolution of homeologs genes during diploidization process.

Evolution of sex hormone binding globulins reveals early gene duplication at the root of vertebrates.

Sex hormone-binding globulin (Shbg) is an important vertebrate blood carrier protein synthetized in the liver and involved in the transport and local regulation of sex steroids in target tissues. A novel shbg gene (shbgb) with a predominant ovarian expression was recently characterized. Being initially found only in salmonids, this shbgb was originally thought to result from the Salmonid-specific whole genome duplication. Using updated transcriptomic and genomic resources we identified Shbgb orthologs in non-salmonid teleosts (European eel, arowana), holosteans (spotted gar, bowfin), polypteriformes (reedfish), agnatha (sea lamprey) and in amphibians, and found that the classical Shbg gene (Shbga) displays a predominant hepatic expression whereas Shbgb has a predominant gonadal expression. Together, these results indicate that these two Shgb genes most likely originate from a whole genome duplication event at the root of vertebrate evolution, followed by numerous and independent losses and by tissue expression specialization of Shbga and Shbgb paralogs.

Comprehensive analysis of polygalacturonase gene family highlights candidate genes related to pollen development and male fertility in wheat (Triticum aestivum L.).

MAIN CONCLUSION: Four polygalacturonase gene family members were highlighted that contribute to elucidate the roles of polygalacturonase during the fertility conversion process in male-sterile wheat. Polygalacturonase (PG) belongs to a large family of hydrolases with important functions in cell separation during plant growth and development via the degradation of pectin. Specific expressed PGs in anthers may be significant for male sterility research and hybrid wheat breeding, but they have not been characterized in wheat (Triticum aestivum L.). In this study, we systematically studied the PG gene family using the latest published wheat reference genomic information. In total, 113 wheat PG genes were identified, which could be classified into six categories A-F according to their structure characteristics and phylogenetic comparisons with Arabidopsis and rice. Polyploidy and segmental duplications in wheat were proved to be mainly responsible for the expansion of the wheat PG gene family. RNA-seq showed that TaPGs have specific temporal and spatial expression characteristics, in which 12 TaPGs with spike-specific expression patterns were detected by qRT-PCR in different fertility anthers of KTM3315A, a thermo-sensitive cytoplasmic male-sterile wheat. Four of them specific upregulated (TaPG09, TaPG95, and TaPG93) or downregulated (TaPG87) at trinucleate stage of fertile anthers, and further aligning with the homologous in Arabidopsis revealed that they may undertake functions such as anther dehiscence, separation of pollen, pollen development, and pollen tube elongation, thereby inducing male fertility conversion in KTM3315A. These findings facilitate function investigations of the wheat PG gene family and provide new insights into the fertility conversion mechanism in male-sterile wheat.

Intestinal microbiota domination under extreme selective pressures characterized by metagenomic read cloud sequencing and assembly.

BACKGROUND: Low diversity of the gut microbiome, often progressing to the point of intestinal domination by a single species, has been linked to poor outcomes in patients undergoing hematopoietic cell transplantation (HCT). Our ability to understand how certain organisms attain intestinal domination over others has been restricted in part by current metagenomic sequencing technologies that are typically unable to reconstruct complete genomes for individual organisms present within a sequenced microbial community. We recently developed a metagenomic read cloud sequencing and assembly approach that generates improved draft genomes for individual organisms compared to conventional short-read sequencing and assembly methods. Herein, we applied metagenomic read cloud sequencing to four stool samples collected longitudinally from an HCT patient preceding treatment and over the course of heavy antibiotic exposure. RESULTS: Characterization of microbiome composition by taxonomic classification of reads reveals that that upon antibiotic exposure, the subject's gut microbiome experienced a marked decrease in diversity and became dominated by Escherichia coli. While diversity is restored at the final time point, this occurs without recovery of the original species and strain-level composition. Draft genomes for individual organisms within each sample were generated using both read cloud and conventional assembly. Read clouds were found to improve the completeness and contiguity of genome assemblies compared to conventional assembly. Moreover, read clouds enabled the placement of antibiotic resistance genes present in multiple copies both within a single draft genome and across multiple organisms. The occurrence of resistance genes associates with the timing of antibiotics administered to the patient, and comparative genomic analysis of the various intestinal E. coli strains across time points as well as the bloodstream isolate showed that the subject's E. coli bloodstream infection likely originated from the intestine. The E. coli genome from the initial pre-transplant stool sample harbors 46 known antimicrobial resistance genes, while all other species from the pre-transplant sample each contain at most 5 genes, consistent with a model of heavy antibiotic exposure resulting in selective outgrowth of the highly antibiotic-resistant E. coli. CONCLUSION: This study demonstrates the application and utility of metagenomic read cloud sequencing and assembly to study the underlying strain-level genomic factors influencing gut microbiome dynamics under extreme selective pressures in the clinical context of HCT.

Distinguishing successive ancient polyploidy levels based on genome-internal syntenic alignment.

BACKGROUND: A basic tool for studying the polyploidization history of a genome, especially in plants, is the distribution of duplicate gene similarities in syntenically aligned regions of a genome. This distribution can usually be decomposed into two or more components identifiable by peaks, or local maxima, each representing a different polyploidization event. The distributions may be generated by means of a discrete time branching process, followed by a sequence divergence model. The branching process, as well as the inference of fractionation rates based on it, requires knowledge of the ploidy level of each event, which cannot be directly inferred from the pair similarity distribution. RESULTS: For a sequence of two events of unknown ploidy, either tetraploid, giving rise to whole genome doubling (WGD), or hexaploid, giving rise to whole genome tripling (WGT), we base our analysis on triples of similar genes. We calculate the probability of the four triplet types with origins in one or the other event, or both, and impose a mutational model so that the distribution resembles the original data. Using a ML transition point in the similarities between the two events as a discriminator for the hypothesized origin of each similarity, we calculate the predicted number of triplets of each type for each model combining WGT and/or WGD. This yields a predicted profile of triplet types for each model. We compare the observed and predicted triplet profiles for each model to confirm the polyploidization history of durian, poplar and cabbage. CONCLUSIONS: We have developed a way of inferring the ploidy of up to three successive WGD and/or WGT events by estimating the time of origin of each of the similarities in triples of genes. This may be generalized to a larger number of events and to higher ploidies.

Avian Expression Patterns and Genomic Mapping Implicate Leptin in Digestion and TNF in Immunity, Suggesting That Their Interacting Adipokine Role Has Been Acquired Only in Mammals

In mammals, leptin and tumor-necrosis factor (TNF) are prominent interacting adipokines mediating appetite control and insulin sensitivity. While TNF pleiotropically functions in immune defense and cell survival, leptin is largely confined to signaling energy stores in adipocytes. Knowledge about the function of avian leptin and TNF is limited and they are absent or lowly expressed in adipose, respectively. Employing radiation-hybrid mapping and FISH-TSA, we mapped TNF and its syntenic genes to chicken chromosome 16 within the major histocompatibility complex (MHC) region. This mapping position suggests that avian TNF has a role in regulating immune response. To test its possible interaction with leptin within the immune system and beyond, we compared the transcription patterns of TNF, leptin and their cognate receptors obtained by meta-analysis of GenBank RNA-seq data. While expression of leptin and its receptor (LEPR) were detected in the brain and digestive tract, TNF and its receptor mRNAs were primarily found in viral-infected and LPS-treated leukocytes. We confirmed leptin expression in the duodenum by immunohistochemistry staining. Altogether, we suggest that whereas leptin and TNF interact as adipokines in mammals, in birds, they have distinct roles. Thus, the interaction between leptin and TNF may be unique to mammals.

Syntenies Between Cohosted Mitochondrial, Chloroplast, and Phycodnavirus Genomes: Functional Mimicry and/or Common Ancestry?

Recent analyses suggest bacterial and/or mitochondrion-like ancestry for giant viruses (Megavirales sensu latu): amoeban mitochondrial gene arrangements resemble those of their candidate homologs in megaviral genomes. This presumed ancestral synteny decreases with genome size across megaviral families at large and within Poxviridae. In this study, analyses focus on Phycodnaviridae, a polyphyletic group of giant viruses infecting Haplophyta, Stramenopiles, and other algae, using syntenies between algal mitogene arrangements and chloroplast genomes and Rickettsia prowazekii as positive controls. Mitogene alignment qualities with Rickettsia are much higher than with viral genomes. Mitogenome synteny with some viruses is higher, for others lower than with Rickettsia, despite lower alignments qualities. In some algae, syntenies among cohosted chloroplast, virus, and mitochondrion are higher, in others lower than expected. This suggests gene order coevolution in cohosted genomes, different coregulations of organelle metabolisms for different algae, and viral mitogenome mimicry, to hijack organelle-committed cellular resources and/or escape cellular defenses/genetic immunity systems. This principle might explain high synteny between human mitochondria and the pathogenic endocellular alphaproteobacterium R. prowazekii beyond common ancestry. Results indicate that putative bacteria/mitochondrion-like genomic ancestors of Phycodnaviridae originated before or at the mitochondrion-bacteria split, and ulterior functional constraints on gene arrangements of cohosted genomes.

Integrative Meta-Assembly Pipeline (IMAP): Chromosome-level genome assembler combining multiple de novo assemblies.

BACKGROUND: Genomic data have become major resources to understand complex mechanisms at fine-scale temporal and spatial resolution in functional and evolutionary genetic studies, including human diseases, such as cancers. Recently, a large number of whole genomes of evolving populations of yeast (Saccharomyces cerevisiae W303 strain) were sequenced in a time-dependent manner to identify temporal evolutionary patterns. For this type of study, a chromosome-level sequence assembly of the strain or population at time zero is required to compare with the genomes derived later. However, there is no fully automated computational approach in experimental evolution studies to establish the chromosome-level genome assembly using unique features of sequencing data. METHODS AND RESULTS: In this study, we developed a new software pipeline, the integrative meta-assembly pipeline (IMAP), to build chromosome-level genome sequence assemblies by generating and combining multiple initial assemblies using three de novo assemblers from short-read sequencing data. We significantly improved the continuity and accuracy of the genome assembly using a large collection of sequencing data and hybrid assembly approaches. We validated our pipeline by generating chromosome-level assemblies of yeast strains W303 and SK1, and compared our results with assemblies built using long-read sequencing and various assembly evaluation metrics. We also constructed chromosome-level sequence assemblies of S. cerevisiae strain Sigma1278b, and three commonly used fungal strains: Aspergillus nidulans A713, Neurospora crassa 73, and Thielavia terrestris CBS 492.74, for which long-read sequencing data are not yet available. Finally, we examined the effect of IMAP parameters, such as reference and resolution, on the quality of the final assembly of the yeast strains W303 and SK1. CONCLUSIONS: We developed a cost-effective pipeline to generate chromosome-level sequence assemblies using only short-read sequencing data. Our pipeline combines the strengths of reference-guided and meta-assembly approaches. Our pipeline is available online at http://github.com/jkimlab/IMAP including a Docker image, as well as a Perl script, to help users install the IMAP package, including several prerequisite programs. Users can use IMAP to easily build the chromosome-level assembly for the genome of their interest.

A comprehensive genomic scan reveals gene dosage balance impacts on quantitative traits in Populus trees.

Gene dosage variation and the associated changes in gene expression influence a wide variety of traits, ranging from cancer in humans to yield in plants. It is also expected to affect important traits of ecological and agronomic importance in forest trees, but this variation has not been systematically characterized or exploited. Here we performed a comprehensive scan of the Populus genome for dosage-sensitive loci affecting quantitative trait variation for spring and fall phenology and biomass production. The study population was a large collection of clonally propagated F1 hybrid lines of Populus that saturate the genome 10-fold with deletions and insertions (indels) of known sizes and positions. As a group, the phenotypic means of the indel lines consistently differed from control nonindel lines, with an overall negative effect of both insertions and deletions on all biomass-related traits but more diverse effects and an overall wider phenotypic distribution of the indel lines for the phenology-related traits. We also investigated the correlation between gene dosage at specific chromosomal locations and phenotype, to identify dosage quantitative trait loci (dQTL). Such dQTL were detected for most phenotypes examined, but stronger effect dQTL were identified for the phenology-related traits than for the biomass traits. Our genome-wide screen for dosage sensitivity in a higher eukaryote demonstrates the importance of global genomic balance and the impact of dosage on life history traits.

Impact of Chromosomal Rearrangements on the Interpretation of Lupin Karyotype Evolution.

Plant genome evolution can be very complex and challenging to describe, even within a genus. Mechanisms that underlie genome variation are complex and can include whole-genome duplications, gene duplication and/or loss, and, importantly, multiple chromosomal rearrangements. Lupins (Lupinus) diverged from other legumes approximately 60 mya. In contrast to New World lupins, Old World lupins show high variability not only for chromosome numbers (2n = 32⁻52), but also for the basic chromosome number (x = 5⁻9, 13) and genome size. The evolutionary basis that underlies the karyotype evolution in lupins remains unknown, as it has so far been impossible to identify individual chromosomes. To shed light on chromosome changes and evolution, we used comparative chromosome mapping among 11 Old World lupins, with Lupinusangustifolius as the reference species. We applied set of L.angustifolius-derived bacterial artificial chromosome clones for fluorescence in situ hybridization. We demonstrate that chromosome variations in the species analyzed might have arisen from multiple changes in chromosome structure and number. We hypothesize about lupin karyotype evolution through polyploidy and subsequent aneuploidy. Additionally, we have established a cytogenomic map of L.angustifolius along with chromosome markers that can be used for related species to further improve comparative studies of crops and wild lupins.

Studies Into ß-Glucan Recognition in Fish Suggests a Key Role for the C-Type Lectin Pathway.

Immune-modulatory effects of ß-glucans are generally considered beneficial to fish health. Despite the frequent application of ß-glucans in aquaculture practice, the exact receptors and downstream signalling remains to be described for fish. In mammals, Dectin-1 is a member of the C-type lectin receptor (CLR) family and the best-described receptor for ß-glucans. In fish genomes, no clear homologue of Dectin-1 could be identified so far. Yet, in previous studies we could activate carp macrophages with curdlan, considered a Dectin-1-specific ß-(1,3)-glucan ligand in mammals. It was therefore proposed that immune-modulatory effects of ß-glucan in carp macrophages could be triggered by a member of the CLR family activating the classical CLR signalling pathway, different from Dectin-1. In the current study, we used primary macrophages of common carp to examine immune modulation by ß-glucans using transcriptome analysis of RNA isolated 6 h after stimulation with two different ß-glucan preparations. Pathway analysis of differentially expressed genes (DEGs) showed that both ß-glucans regulate a comparable signalling pathway typical of CLR activation. Carp genome analysis identified 239 genes encoding for proteins with at least one C-type Lectin Domains (CTLD). Narrowing the search for candidate ß-glucan receptors, based on the presence of a conserved glucan-binding motif, identified 13 genes encoding a WxH sugar-binding motif in their CTLD. These genes, however, were not expressed in macrophages. Instead, among the ß-glucan-stimulated DEGs, a total of six CTLD-encoding genes were significantly regulated, all of which were down-regulated in carp macrophages. Several candidates had a protein architecture similar to Dectin-1, therefore potential conservation of synteny of the mammalian Dectin-1 region was investigated by mining the zebrafish genome. Partial conservation of synteny with a region on the zebrafish chromosome 16 highlighted two genes as candidate ß-glucan receptor. Altogether, the regulation of a gene expression profile typical of a signalling pathway associated with CLR activation and, the identification of several candidate ß-glucan receptors, suggest that immune-modulatory effects of ß-glucan in carp macrophages could be a result of signalling mediated by a member of the CLR family.

Characterization of channel catfish (Ictalurus punctatus) melanocortin-3 receptor reveals a potential network in regulation of energy homeostasis.

The melanocortin-3 receptor (MC3R) is known to be involved in regulation of energy homeostasis, regulating feed efficiency and nutrient partitioning in mammals. Its physiological roles in non-mammalian vertebrates, especially economically important aquaculture species, are not well understood. Channel catfish (Ictalurus punctatus) is the main freshwater aquaculture species in North America. In this study, we characterized the channel catfish MC3R. The mc3r of channel catfish encoded a putative protein (ipMC3R) of 367 amino acids. We transfected HEK293T cells with ipMC3R plasmid for functional studies. Five agonists, including adrenocorticotropin, α-melanocyte stimulating hormone (α-MSH), ß-MSH, [Nle4, D-Phe7]-α-MSH, and D-Trp8-γ-MSH, were used in the pharmacological studies. Our results showed that ipMC3R bound ß-MSH with higher affinity and D-Trp8-γ-MSH with lower affinity compared with human MC3R. All agonists could stimulate ipMC3R and increase intracellular cAMP production with sub-nanomolar potencies. The extracellular signal-regulated kinases 1 and 2 (ERK1/2) activation could also be triggered by ipMC3R. The ipMC3R exhibited constitutive activities in both cAMP and ERK1/2 pathways, and Agouti-related protein served as an inverse agonist at ipMC3R, potently inhibiting the high basal cAMP level. Moreover, we showed that melanocortin receptor accessory protein 2 (MRAP2) preferentially modulated ipMC3R in cAMP production rather than ERK1/2 activation. Our study will assist further investigation of the physiological roles of the ipMC3R, especially in energy homeostasis, in channel catfish.