Bulk segregant RNA-seq reveals expression and positional candidate genes and allele-specific expression for disease resistance against enteric septicemia of catfish

Archetypal Analysis of Kidney Allograft Biopsies Using Next-generation Sequencing Technology

BACKGROUND

In kidney transplantation, molecular diagnostics may be a valuable approach to improve the precision of the diagnosis. Using next-generation sequencing (NGS), we aimed to identify clinically relevant archetypes.

METHODS

We conducted an Illumina bulk RNA sequencing on 770 kidney biopsies (540 kidney recipients) collected between 2006 and 2021 from 11 European centers. Differentially expressed genes were determined for 11 Banff lesions. An ElasticNet model was used for feature selection, and 4 machine learning classifiers were trained to predict the probability of presence of the lesions. NGS-based classifiers were used in an unsupervised archetypal analysis to different archetypes. The association of the archetypes with allograft survival was assessed using the iBox risk prediction score.

RESULTS

The ElasticNet feature selection reduced the number of the genes from a range of 859-10 830 to a range of 52-867 genes. NGS-based classifiers demonstrated robust performances (precision-recall area under the curves 0.708-0.980) in predicting the Banff lesions. Archetypal analysis revealed 8 distinct phenotypes, each characterized by distinct clinical, immunological, and histological features. Although the archetypes confirmed the well-defined Banff rejection phenotypes for T cell-mediated rejection and antibody-mediated rejection, equivocal histologic antibody-mediated rejection, and borderline diagnoses were reclassified into different archetypes based on their molecular signatures. The 8 NGS-based archetypes displayed distinct allograft survival profiles with incremental graft loss rates between archetypes, ranging from 90% to 56% rates 7 y after evaluation ( P < 0.0001).

CONCLUSIONS

Using molecular phenotyping, 8 archetypes were identified. These NGS-based archetypes might improve disease characterization, reclassify ambiguous Banff diagnoses, and enable patient-specific risk stratification.

Advancing genetic improvement in the omics era: status and priorities for United States aquaculture

BACKGROUND

The innovations of the "Omics Era" have ushered in significant advancements in genetic improvement of agriculturally important animal species through transforming genetics, genomics and breeding strategies. These advancements were often coordinated, in part, by support provided over 30 years through the 1993-2023 National Research Support Project 8 (NRSP8, National Animal Genome Research Program, NAGRP) and affiliate projects focused on enabling genomic discoveries in livestock, poultry, and aquaculture species. These significant and parallel advances demand strategic planning of future research priorities. This paper, as an output from the May 2023 Aquaculture Genomics, Genetics, and Breeding Workshop, provides an updated status of genomic resources for United States aquaculture species, highlighting major achievements and emerging priorities.

MAIN TEXT

Finfish and shellfish genome and omics resources enhance our understanding of genetic architecture and heritability of performance and production traits. The 2023 Workshop identified present aims for aquaculture genomics/omics research to build on this progress: (1) advancing reference genome assembly quality; (2) integrating multi-omics data to enhance analysis of production and performance traits; (3) developing resources for the collection and integration of phenomics data; (4) creating pathways for applying and integrating genomics information across animal industries; and (5) providing training, extension, and outreach to support the application of genome to phenome. Research focuses should emphasize phenomics data collection, artificial intelligence, identifying causative relationships between genotypes and phenotypes, establishing pathways to apply genomic information and tools across aquaculture industries, and an expansion of training programs for the next-generation workforce to facilitate integration of genomic sciences into aquaculture operations to enhance productivity, competitiveness, and sustainability.

CONCLUSION

This collective vision of applying genomics to aquaculture breeding with focus on the highlighted priorities is intended to facilitate the continued advancement of the United States aquaculture genomics, genetics and breeding research community and industries. Critical challenges ahead include the practical application of genomic tools and analytical frameworks beyond academic and research communities that require collaborative partnerships between academia, government, and industry. The scope of this review encompasses the use of omics tools and applications in the study of aquatic animals cultivated for human consumption in aquaculture settings throughout their life-cycle.

Impacts of Working Hours, Wages, and Regular Employment Opportunity on Suicide Mortalities of Employed and Unemployed Individuals before and during the COVID-19 Pandemic in Japan

Standardized suicide mortality rates per 100,000 population (SMRs) in Japan consistently decreased from 2009-2019, but these decreasing trends were reversed to increase in 2020. To clarify the mechanisms of recent increasing suicide in Japan, temporal fluctuations of SMRs disaggregated by sex and employment status (employed and unemployed individuals) and labor indices such as working hours, wages, and regular employment opportunity index (REO) from January 2012 to June 2023 were analyzed using interrupted time-series analysis. Additionally, temporal causalities from labor indices to SMRs were analyzed using vector autoregressive and non-linear auto-regressive distributed lag analyses. Decreasing trends among employed SMRs of both sexes were attenuated after the enactment of the "Work Style Reform Program" in 2018, but male SMRs were unaffected by the COVID-19 pandemic. However, female employed SMRs sharply increased, synchronized with the "Work Style Reform Act" and the COVID-19 pandemic outbreak (the COVID-19 impact was greater than the "Work Style Reform Act"). Additionally, unemployed SMRs of both sexes sharply increased with the revision and scale-down of countermeasures against economic deterioration caused by COVID-19 ("revision of economic supportive countermeasures against economic deterioration caused by COVID-19"). Unexpectedly, after enacting the "Work Style Reform Act", wages decreased due to possibly decreasing working hours. Increasing REO, which consistently increased, was a protective factor for male suicides, but unemployed SMRs were not affected by any labor indices. It has been established that controlling a heavy workload plays an important role in suppressing the deterioration of physical and mental conditions, including suicide; however, this study suggested that, at least within appropriate ranges of working hours, decreasing working hours due to excessive management probably contributes to increasing suicides of some vulnerable individuals via de-creasing their wages. Although governmental welfare and economic support measures had to be revised according to rapidly changing situations during the COVID-19 pandemic, this study also suggested that temporal gaps among a part of revisions of several welfare and economic support measures were unexpectedly involved in drastically/sharply increasing suicides of unemployed individuals in 2022.

Transcriptome analysis of Macrobrachium rosenbergii hemocytes in response to Staphylococcus aureus infection

The aquaculture industry has suffered significant financial losses as a result of disease outbreaks. In particular, disease outbreaks have become a major problem that can seriously affect the sustainable development of the Macrobrachium rosenbergii aquaculture industry. It is crucial to determine the defense mechanism of the host after pathogenic invasion in order to provide effective defense measures after disease outbreaks. Shrimp, like other invertebrates, primarily depend on their innate immune systems to defend against pathogens, and recognize and resist pathogens through humoral and cellular immune responses. In this investigation, we used RNA-seq technology to investigate the transcriptome of hemocytes from M. rosenbergii induced by Staphylococcus aureus. Our main targets were immune pathways and genes related to innate immunity. RNA-seq identified 209,069 and 204,775 unigenes in the control and experimental groups, respectively. In addition, we identified 547 and 1734 differentially expressed genes (DEGs) following S. aureus challenge after 6 and 12 h (h), respectively. GO and KEGG enrichment analysis revealed that the DEGs were significantly enriched in several biological signalling pathways, including NOD-like receptor, PI3K-Akt, Toll and Imd, IL-17, TGF-beta, RIG-I-like receptor, cAMP, apoptosis, and C-type lectin receptor. Sixteen DEGs were chosen at random for qPCR verification; these results concurred with those from sequencing. Our findings revealed that immune-related genes play an important role in antibacterial activities and have specific functions for gram-positive bacteria. These results provide more data for the prevention of M. rosenbergii diseases and offer a basis for the better prevention of diseases.

Applying genetic technologies to combat infectious diseases in aquaculture

Disease and parasitism cause major welfare, environmental and economic concerns for global aquaculture. In this review, we examine the status and potential of technologies that exploit genetic variation in host resistance to tackle this problem. We argue that there is an urgent need to improve understanding of the genetic mechanisms involved, leading to the development of tools that can be applied to boost host resistance and reduce the disease burden. We draw on two pressing global disease problems as case studies-sea lice infestations in salmonids and white spot syndrome in shrimp. We review how the latest genetic technologies can be capitalised upon to determine the mechanisms underlying inter- and intra-species variation in pathogen/parasite resistance, and how the derived knowledge could be applied to boost disease resistance using selective breeding, gene editing and/or with targeted feed treatments and vaccines. Gene editing brings novel opportunities, but also implementation and dissemination challenges, and necessitates new protocols to integrate the technology into aquaculture breeding programmes. There is also an ongoing need to minimise risks of disease agents evolving to overcome genetic improvements to host resistance, and insights from epidemiological and evolutionary models of pathogen infestation in wild and cultured host populations are explored. Ethical issues around the different approaches for achieving genetic resistance are discussed. Application of genetic technologies and approaches has potential to improve fundamental knowledge of mechanisms affecting genetic resistance and provide effective pathways for implementation that could lead to more resistant aquaculture stocks, transforming global aquaculture.

Identification of genes for seed isoflavones based on bulk segregant analysis sequencing in soybean natural population

We identified four hub genes for isoflavone biosynthesis based on BSA-seq and WGCNA methods and validated that GmIE3-1 positively contribute to isoflavone accumulation in soybean. Soybean isoflavones are secondary metabolites of great interest owing to their beneficial impact on human health. Herein, we profiled the seed isoflavone content by HPLC in 1551 soybean accessions grown in two locations for two years and constructed two extreme pools with high (4065.1 µg g^-1) and low (1427.23 µg g^-1) isoflavone contents to identify candidate genes involved in isoflavone biosynthesis pathways using bulk segregant analysis sequencing (BSA-seq) approach. The results showed that the average sequencing depths were 50.3× and 65.7× in high and low pools, respectively. A total of 23,626 polymorphic SNPs and 5299 InDels were detected between both pools and 1492 genes with different variations were identified. Based on differential genes in BSA-seq and weighted gene co-expression network analysis (WGCNA), four hub genes, Glyma.06G290400 (designated as GmIE3-1), Glyma.01G239200, Glyma.01G241500, Glyma.13G256100 were identified, encoding E3 ubiquitin-protein ligase, arm repeat protein interacting with ABF2, zinc metallopeptidase EGY3, and dynamin-related protein 3A, respectively. The allelic variation in GmIE3-1 showed a significant influence on isoflavone accumulation. The virus-induced gene silencing (VIGS) and RNAi hairy root transformation of GmIE3-1 revealed partial suppression of this gene could cause a significant decrease (P < 0.0001) of total isoflavone content, suggesting GmIE3-1 is a positive regulator for isoflavones. The present study demonstrated that the BSA-seq approach combined with WGCNA, VIGS and hairy root transformation can efficiently identify isoflavone candidate genes in soybean natural population.

Transcript Characteristics on the Susceptibility Difference of Bovine Respiratory Disease

Bovine respiratory disease (BRD) is one of the major health issues in the cattle industry, resulting in significant financial crises globally. There is currently no good treatment, and cattle are made resistant to pneumonia through disease-resistant breeding. The serial blood samples from six Xinjiang brown (XJB) calves were collected for the RNA sequencing (RNA-seq). The obtained six samples were grouped into two groups, in each group as infected with BRD and healthy calves, respectively. In our study, the differential expression mRNAs were detected by using RNA-seq and constructed a protein-protein interaction (PPI) network related to the immunity in cattle. The key genes were identified by protein interaction network analysis, and the results from RNA-seq were verified using reverse transcription-quantitative polymerase chain reaction (RT-qPCR). A total of 488 differentially expressed (DE) mRNAs were identified. Importantly, the enrichment analysis of these identified DEGs classified them as mainly enriched in the regulation and immune response processes. The 16 hub genes were found to be related to immune pathways categorized by PPIs analysis. Results revealed that many hub genes were related to the immune response to respiratory disease. These results will provide the basis for a better understanding of the molecular mechanism of bovine resistance to BRD.

Harnessing the potential of bulk segregant analysis sequencing and its related approaches in crop breeding

Most plant traits are governed by polygenes including both major and minor genes. Linkage mapping and positional cloning have contributed greatly to mapping genomic loci controlling important traits in crop species. However, they are low-throughput, time-consuming, and have low resolution due to which their efficiency in crop breeding is reduced. In this regard, the bulk segregant analysis sequencing (BSA-seq) and its related approaches, viz., quantitative trait locus (QTL)-seq, bulk segregant RNA-Seq (BSR)-seq, and MutMap, have emerged as efficient methods to identify the genomic loci/QTLs controlling specific traits at high resolution, accuracy, reduced time span, and in a high-throughput manner. These approaches combine BSA with next-generation sequencing (NGS) and enable the rapid identification of genetic loci for qualitative and quantitative assessments. Many previous studies have shown the successful identification of the genetic loci for different plant traits using BSA-seq and its related approaches, as discussed in the text with details. However, the efficiency and accuracy of the BSA-seq depend upon factors like sequencing depth and coverage, which enhance the sequencing cost. Recently, the rapid reduction in the cost of NGS together with the expected cost reduction of third-generation sequencing in the future has further increased the accuracy and commercial applicability of these approaches in crop improvement programs. This review article provides an overview of BSA-seq and its related approaches in crop breeding together with their merits and challenges in trait mapping.

Bulked Segregant RNA Sequencing Revealed Difference Between Virulent and Avirulent Brown Planthoppers

The brown planthopper (Nilaparvata lugens Stål, BPH) is one of the most devastating insect pests of rice (Oryza sativa L.), but BPH populations have varying degrees of virulence to rice varieties carrying different resistance genes. To help efforts to characterize these variations we applied bulked segregant RNA sequencing (BSR-seq) to identify differentially expressed genes (DEGs) and genetic loci associated with BPH virulence to YHY15 rice plants carrying the resistance gene Bph15. BPHs that are highly virulent or avirulent to these plants were selected from an F2 population to form two contrasting bulks, and BSR-seq identified 751 DEGs between the bulks. Genes associated with carbohydrate, amino acid and nucleotide metabolism, the endocrine system, and signal transduction were upregulated in the avirulent insects when they fed on these plants. The results also indicated that shifts in lipid metabolism and digestive system pathways were crucial for the virulent BPHs' adaptation to the resistant rice. We identified 24 single-nucleotide polymorphisms (SNPs) in 21 genes linked with BPH virulence. Possible roles of genes apparently linked to BPH virulence are discussed. Our results provide potentially valuable information for further studies of BPH virulence mechanisms and development of robust control strategies.

Identification of Candidate Genes Associated With Hypoxia Tolerance in Trachinotus blochii Using Bulked Segregant Analysis and RNA-Seq

Golden Pompano (Trachinotus blochii) has rapidly developed into the one of the main valuable fish species in Chinese marine aquaculture. Due to its rapid growth, active metabolism, and high oxygen consumption, hypoxia will increase its mortality and cause serious economic losses. We constructed two experimental groups of fish with different degrees of tolerance to hypoxia, used BSR-Seq analysis based on genome and genetic linkage groups to locate SNPs and genes that were related to the differences in hypoxia tolerance. The results showed that hypoxia tolerance SNPs of golden pompano may be jointly determined by multiple linkage groups, especially linkage groups 18 and 22. There were 768 and 348 candidate genes located in the candidate regions of the brain and liver, respectively. These genes were mainly involved in anaerobic energy metabolism, stress response, immune response, waste discharge, and cell death. The prostaglandin-endoperoxide synthase 2 (PTGS2) on LG8, which is involved in the metabolism of arachidonic acid, has a G/A nonsynonymous mutation at position 20641628, and the encoded amino acid was changed from hydrophobic aspartic acid to asparaginate. The specific pathway of the RIG-I-like receptor signaling pathway in the liver may mediate the metabolic system and the immune system, linking glucose metabolism with immune regulation. The death of the hypoxia-intolerant group may be due to the accumulation of lactic acid caused by the activation of anaerobic glycolysis during the early stage of hypoxia stress, and the activation of type I interferon was inhibited, which resulted in decreased immunity. Among the genes involved in the RIG-I-like receptor signaling pathway, the CYLD Lysine 63 Deubiquitinase (CYLD) located on LG16 had a G/T nonsynonymous mutation at position 13629651, and the encoded amino acid was changed from alanine acid to valine. The interferon induced with helicase C domain 1 (Ifih1) located on LG18 has a G/C nonsynonymous mutation at position 16153700, and the encoded hydrophilic glycine was changed to hydrophobic alanine. Our findings suggest these SNPs may assist in the molecular breeding of hypoxia-tolerant golden pompano, and speculate that the balance of glucose and lipid metabolism plays a key role in Trachinotus blochii under acute hypoxia.

Construction of a chromosome-level genome assembly for genome-wide identification of growth-related quantitative trait loci in Sinocyclocheilus grahami (Cypriniformes, Cyprinidae)

The Dianchi golden-line barbel, Sinocyclocheilus grahami (Regan, 1904), is one of the “Four Famous Fishes” of Yunnan Province, China. Given its economic value, this species has been artificially bred successfully since 2007, with a nationally selected breed (“S. grahami, Bayou No. 1”) certified in 2018. For the future utilization of this species, its growth rate, disease resistance, and wild adaptability need to be improved, which could be achieved with the help of molecular marker-assisted selection (MAS). In the current study, we constructed the first chromosome-level genome of S. grahami, assembled 48 pseudo-chromosomes, and obtained a genome assembly of 1.49 Gb. We also performed QTL-seq analysis of S. grahami using the highest and lowest bulks (i.e., largest and smallest size) in both a sibling and random population. We screened two quantitative trait loci (QTLs) (Chr3, 14.9–39.1 Mb and Chr17, 4.1–27.4 Mb) as the major growth-related locations. Several candidate genes (e.g., map2k5, stat1, phf21a, sox6, and smad6) were also identified, with functions related to growth, such as cell differentiation, neuronal development, skeletal muscle development, chondrogenesis, and immunity. These results built a solid foundation for in-depth MAS studies on the growth traits of S. grahami.

The loss of function of HEL, which encodes a cellulose synthase interactive protein, causes helical and vine-like growth of tomato

Helical growth is an economical way for plant to obtain resources. The classic microtubule-microfibril alignment model of Arabidopsis helical growth involves restriction of the appropriate orientation of cellulose microfibrils appropriately in the cell walls. However, the molecular mechanism underlying tomato helical growth remains unknown. Here, we identified a spontaneous tomato helical (hel) mutant with right-handed helical cotyledons and petals but left-handed helical stems and true leaves. Genetic analysis revealed that the hel phenotype was controlled by a single recessive gene. Using map-based cloning, we cloned the HEL gene, which encodes a cellulose interacting protein homologous to CSI1 of Arabidopsis. We identified a 27 bp fragment replacement that generated a premature stop codon. Transgenic experiments showed that the helical growth phenotype could be restored by the allele of this gene from wild-type Pyriforme. In contrast, the knockout mutation of HEL in Pyriforme via CRISPR/Cas9 resulted in helical growth. These findings shed light on the molecular control of the helical growth of tomato.

Weeping candidate genes screened using comparative transcriptomic analysis of weeping and upright progeny in an F1 population of Prunus mume

Weeping is a specific plant architecture with high ornamental value. Despite the considerable importance of the weeping habit to landscaping applications and knowledge of plant architecture biology, little is known regarding the underlying molecular mechanisms. In this study, growth and phytohormone content were analyzed among the progeny of different branch types in an F1 mapping population of Prunus mume with varying plant architecture. Bulked segregant RNA sequencing was conducted to compare differences among progeny at a transcriptional level. The weeping habit appears to be a complex process regulated by a series of metabolic pathways, with photosynthesis and flavonoid biosynthesis highly enriched in differentially expressed genes between weeping and upright progeny. Based on functional annotation and homologous analyses, we identified 30 candidate genes related to weeping that merit further analysis, including 10 genes related to IAA and GA₃ biosynthesis, together with 6 genes related to secondary branch growth. The results of this study will facilitate further studies of the molecular mechanisms underlying the weeping habit in P. mume.

Bulked Segregant RNA-Seq Reveals Distinct Expression Profiling in Chinese Wheat Cultivar Jimai 23 Responding to Powdery Mildew

Wheat powdery mildew, caused by Blumeria graminis f. sp. tritici (Bgt), is one of the most destructive fungal diseases threatening global wheat production. Host resistance is well known to be the most efficient method to control this disease. However, the molecular mechanism of wheat powdery mildew resistance (Pm) is still unclear. To analyze the molecular mechanism of Pm, we used the resistant wheat cultivar Jimai 23 to investigate its potential resistance components and profiled its expression in response to powdery mildew infection using bulked segregant RNA-Seq (BSR-Seq). We showed that the Pm of Jimai 23 was provided by a single dominant gene, tentatively designated PmJM23, and assigned it to the documented Pm2 region of chromosome 5DS. 3,816 consistently different SNPs were called between resistant and susceptible parents and the bulked pools derived from the combinations between the resistant parent Jimai23 and the susceptible parent Tainong18. 58 of the SNPs were assigned to the candidate region of PmJM23. Subsequently, 3,803 differentially expressed genes (DEGs) between parents and bulks were analyzed by GO, COG and KEGG pathway enrichment. The temporal expression patterns of associated genes following Bgt inoculation were further determined by RT-qPCR. Expression of six disease-related genes was investigated during Bgt infection and might serve as valuable genetic resources for the improvement of durable resistance to Bgt.

Transcriptome Profiling of Pacu (Piaractus mesopotamicus) Challenged With Pathogenic Aeromonas hydrophila: Inference on Immune Gene Response

Pacu (Piaractus mesopotamicus) is a Neotropical fish of major importance for South American aquaculture. Septicemia caused by Aeromonas hydrophila bacteria is currently considered a substantial threat for pacu aquaculture that have provoked infectious disease outbreaks with high economic losses. The understanding of molecular aspects on progress of A. hydrophila infection and pacu immune response is scarce, which have limited the development of genomic selection for resistance to this infection. The present study aimed to generate information on transcriptome of pacu in face of A. hydrophila infection, and compare the transcriptomic responses between two groups of time-series belonging to a disease resistance challenge, peak mortality (HM) and mortality plateau (PM) groups of individuals. Nine RNA sequencing (RNA-Seq) libraries were prepared from liver tissue of challenged individuals, generating ∼160 million 150 bp pair-end reads. After quality trimming/cleanup, these reads were assembled de novo generating 211,259 contigs. When the expression of genes from individuals of HM group were compared to individuals from control group, a total of 4,413 differentially expressed transcripts were found (2,000 upregulated and 2,413 downregulated candidate genes). Additionally, 433 transcripts were differentially expressed when individuals from MP group were compared with those in the control group (155 upregulated and 278 downregulated candidate genes). The resulting differentially expressed transcripts were clustered into the following functional categories: cytokines and signaling, epithelial protection, antigen processing and presentation, apoptosis, phagocytosis, complement system cascades and pattern recognition receptors. The proposed results revealing relevant differential gene expression on HM and PM groups which will contribute to a better understanding of the molecular defense mechanisms during A. hydrophila infection.

Transcriptome analysis of grass carp (Ctenopharyngodon idella) between fast- and slow-growing fish

Grass carp is one of the most important freshwater aquaculture species in China. However, the mechanisms underlying the growth of muscle tissue in the fish are unclear. High-throughput RNA-Seq was used to analyze the transcriptome of grass carp muscle tissue between fast- and slow-growing fish family groups. Twenty-four individuals each from 4 fast-growing families and 4 slow-growing families were used to reduce background noise. 71 up-regulated and 35 down-regulated genes were identified in the differentially expressed genes (DEGs). GO and KEGG enrichment analyses revealed the DEGs were involved in the GH/IGF axis, calcium metabolism, protein and glycogen synthesis, oxygen transport, cytoskeletal and myofibrillar components. IGFBP1 was up-regulated in big fish while GHR2 was down-regulated. Glutamic pyruvate transaminase 2, an indicator of liver tissue damage, was down-regulated in big grass carp, which indicates that the fish was better adapted to an artificially formulated diet. GAPDH, the rate-limiting enzyme in glycolytic flux was highly expressed in fast-growing grass carp, reflecting enhanced carbohydrate metabolism. Higher expression of ALAS2 and myoglobin 1 in big grass carp, related to oxygen transport might promote aerobic exercise along with food intake and muscle growth. Genes for cytoskeletal and myofibrillar components such as tropomyosin, meromyosin, and troponin I were also up-regulated in big grass carp. These results provide valuable information about the key genes for use as biomarkers of growth in selective breeding programs for grass carp and contribute to our understanding of the molecular mechanisms and regulative pathways regulating growth in fish.

Sex determination in the GIFT strain of tilapia is controlled by a locus in linkage group 23

Background

Tilapias (Family Cichlidae) are the second most important group of aquaculture species in the world. They have been the subject of much research on sex determination due to problems caused by early maturation in culture and their complex sex-determining systems. Different sex-determining loci (linkage group 1, 20 and 23) have been detected in various tilapia stocks. The ‘genetically improved farmed tilapia’ (GIFT) stock, founded from multiple Nile tilapia (Oreochromis niloticus) populations, with some likely to have been introgressed with O. mossambicus, is a key resource for tilapia aquaculture. The sex-determining mechanism in the GIFT stock was unknown, but potentially complicated due to its multiple origins.

Results

A bulk segregant analysis (BSA) version of double-digest restriction-site associated DNA sequencing (BSA-ddRADseq) was developed and used to detect and position sex-linked single nucleotide polymorphism (SNP) markers in 19 families from the GIFT strain breeding nucleus and two Stirling families as controls (a single XY locus had been previously mapped to LG1 in the latter). About 1500 SNPs per family were detected across the genome. Phenotypic sex in Stirling families showed strong association with LG1, whereas only SNPs located in LG23 showed clear association with sex in the majority of the GIFT families. No other genomic regions linked to sex determination were apparent. This region was validated using a series of LG23-specific DNA markers (five SNPs with highest association to sex from this study, the LG23 sex-associated microsatellite UNH898 and ARO172, and the recently isolated amhy marker for individual fish (n = 284).

Conclusions

Perhaps surprisingly given its multiple origins, sex determination in the GIFT strain breeding nucleus was associated only with a locus in LG23. BSA-ddRADseq allowed cost-effective analysis of multiple families, strengthening this conclusion. This technique has potential to be applied to other complex traits. The sex-linked SNP markers identified will be useful for potential marker-assisted selection (MAS) to control sex-ratio in GIFT tilapia to suppress unwanted reproduction during growout.

A predicted NEDD8 conjugating enzyme gene identified as a Capsicum candidate Rf gene using bulk segregant RNA sequencing

Cytoplasmic male sterility (CMS) is an important tool for producing F₁ hybrids, which can exhibit heterosis. The companion system, restorer-of-fertility (Rf), is poorly understood at the molecular level and would be valuable in producing restorer lines for hybrid seed production. The identity of the Rf gene in Capsicum (pepper) is currently unclear. In this study, using bulked segregant RNA sequencing (BSR-seq), a strong candidate Rf gene, Capana06g002866, which is annotated as a NEDD8 conjugating enzyme E2, was identified. Capana06g002866 has an ORF of 555 bp in length encoding 184 amino acids; it can be cloned from F₁ plants from the hybridization of the CMS line 8A and restorer line R1 but is not found in CMS line 8A. With qRT-PCR validation, Capana06g002866 was found to be upregulated in restorer accessions compared to sterile accessions. The relative expression in flower buds increased with the developmental stage in F₁ plants, while the expression was very low in all flower bud stages of the CMS lines. These results provide new insights into the Rf gene in pepper and will be useful for other crops utilizing the CMS system.

Transcriptional Regulation Analysis of Alzheimer's Disease Based on FastNCA Algorithm

Background:

Understanding the relationship between genetic variation and gene expression is a central issue in genetics. Although many studies have identified genetic variations associated with gene expression, it is unclear how they perturb the underlying regulatory network of gene expression.

Objective:

To explore how genetic variations perturb potential transcriptional regulation networks of Alzheimer’s disease (AD) to paint a more complete picture of the complex landscape of transcription regulation.

Methods:

Fast network component analysis (FastNCA), which can capture the genetic variations in the form of single nucleotide polymorphisms (SNPs), is applied to analyse the expression activities of TFs and their regulatory strengths on TGs using microarray and RNA-seq data of AD. Then, multi-data fusion analysis was used to analyze the different TGs regulated by the same TFs in the different data by constructing the transcriptional regulatory networks of differentially expressed genes.

Results:

the common TF regulating TGs are not necessarily identical in different data, they may be involved in the same pathways that are closely related to the pathogenesis of AD, such as immune response, signal transduction and cytokine-cytokine receptor interaction pathways. Even if they are involved in different pathways, these pathways are also confirmed to have a potential link with AD.

Conclusion:

The study shows that the pathways of different TGs regulated by the same TFs in different data are all closely related to AD. Multi-data fusion analysis can form a certain complement to some extent and get more comprehensive results in the process of exploring the pathogenesis of AD.

Comparative transcriptional analysis of Capsicum flower buds between a sterile flower pool and a restorer flower pool provides insight into the regulation of fertility restoration

BACKGROUND

Cytoplasmic male sterility (CMS) and its restoration of fertility (Rf) system is an important mechanism to produce F₁ hybrid seeds. Understanding the interaction that controls restoration at a molecular level will benefit plant breeders. The CMS is caused by the interaction between mitochondrial and nuclear genes, with the CMS phenotype failing to produce functional anthers, pollen, or male gametes. Thus, understanding the complex processes of anther and pollen development is a prerequisite for understanding the CMS system. Currently it is accepted that the Rf gene in the nucleus restores the fertility of CMS, however the Rf gene has not been cloned. In this study, CMS line 8A and the Rf line R1, as well as a sterile pool (SP) of accessions and a restorer pool (RP) of accessions analyzed the differentially expressed genes (DEGs) between CMS and its fertility restorer using the conjunction of RNA sequencing and bulk segregation analysis.

RESULTS

A total of 2274 genes were up-regulated in R1 as compared to 8A, and 1490 genes were up-regulated in RP as compared to SP. There were 891 genes up-regulated in both restorer accessions, R1 and RP, as compared to both sterile accessions, 8A and SP. Through annotation and expression analysis of co-up-regulated expressed genes, eight genes related to fertility restoration were selected. These genes encode putative fructokinase, phosphatidylinositol 4-phosphate 5-kinase, pectate lyase, exopolygalacturonase, pectinesterase, cellulose synthase, fasciclin-like arabinogalactan protein and phosphoinositide phospholipase C. In addition, a phosphatidylinositol signaling system and an inositol phosphate metabolism related to the fertility restorer of CMS were ranked as the most likely pathway for affecting the restoration of fertility in pepper.

CONCLUSIONS

Our study revealed that eight genes were related to the restoration of fertility, which provides new insight into understanding the molecular mechanism of fertility restoration of CMS in Capsicum.

Screening for differential expression of genes for resistance to Sitodiplosis mosellana in bread wheat via BSR-seq analysis

Five putative candidate genes for OWBM resistance in Chinese winter wheat 'Jimai 24' were identified via BSR-seq and differential expression analyses. Orange wheat blossom midge (OWBM), Sitodiplosis mosellana, is one of the most serious threats to wheat production worldwide. Conventional gene mapping methods to identify genes require significant amounts of financial support and time. Here, bulked segregant RNA-seq (BSR-seq) was applied to profile candidate genes and develop associated markers for OWBM resistance. Previously, we identified a major QTL (QSm.hebau-4A) for OWBM resistance on the long arm of chromosome 4A. In this study, we aimed at screening differentially expressed resistance genes associated with this QTL. Twelve differentially expressed genes (DEGs) were obtained based on BSR-seq and differential expression analyses. Among them, four were confirmed to be associated with OWBM resistance via quantitative reverse transcription PCR, using an additional set of wheat samples subjected to OWBM invasion. One SPI-like gene and one Malectin-like gene were revealed by gene annotation, respectively. Sequencing results confirmed that the four DEGs and the SPI gene had SNP polymorphisms between wheat parents. All these five resistance-related genes for OWBM were located in the same genomic region with QSm.hebau-4A. Furthermore, six new markers developed based on sequences of the five genes were also mapped in the same genomic region using genetic population. These five genes may be the candidate genes for OWBM resistance in Chinese wheat 'Jimai 24' and should be the targets for further positional isolation.

Transcriptome sequencing reveals phagocytosis as the main immune response in the pathogen-challenged sea urchin Strongylocentrotus intermedius

The clarification of host immune responses to causative bacteria of spotting disease in the sea urchin Strongylocentrotus intermedius is vital to preventing and controlling this disease, especially to selective breeding for disease resistance. For this purpose, sea urchins were challenged with the causative bacterium Vibrio sp. to obtain spotting diseased and undiseased samples. We conducted next-generation sequencing to assess the key genes/pathways in control (CG), diseased (DG), and undiseased (UG) groups. A total of 454.1 million clean reads were obtained and assembled into 23,899 UniGenes with an N50 of 1359 bp, with 86.11% of them matching the genome sequence of the sea urchin S. purpuratus. A total of 8415 UniGenes were mapped to the non-redundant database. Salmon expression analysis revealed 725 significantly differentially expressed genes (DEGs) among CG, DG, and UG. These DEGs were enriched into 72 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, including a core set of immune correlated pathways notably in the phagosome, vitamin digestion and absorption, Wnt signaling, and Notch signaling pathways. DG was evidently upregulated in these immune pathways and could enhance phagocytosis directly or indirectly. Thus, phagocytosis was the main coelomic cellular immune response in S. intermedius challenged by spotting disease causative bacterium. The expression patterns of 10 DEGs were confirmed via RT-qPCR, and the expression levels were consistent with the results of RNA-seq. Furthermore, 9899 SSRs were identified, and 123,692, 151,827, and 114,368 candidate SNPs were identified from CG, DG, and UG, respectively. These results provide basic information for our understanding of sea urchin antibacterial immunity.

Whole-genome resequencing from bulked-segregant analysis reveals gene set based association analyses for the Vibrio anguillarum resistance of turbot (Scophthalmus maximus)

Many achievements have been made to develop quantitative trait loci (QTLs) and gene-associated single nucleotide polymorphisms (SNPs) to facilitate practical marker-assisted selection (MAS) in aquatic animals. However, the systematic studies of SNPs associated with extreme threshold traits were poor in populations lacking of parental genomic information. Coupling next generation sequencing with bulked segregant analysis (BSA) should allow identification of numerous associated SNPs with extreme phenotypes. In the present study, using combination of SNP frequency difference and Euclidean distance, we conducted linkage analysis of SNPs located in genes involved in immune responses, and identified markers associated with Vibrio anguillarum resistance in turbot (Scophthalmus maximus). A total of 221 SNPs was found as candidate SNPs between resistant and susceptible individuals. Among these SNPs, 35 loci located in immune related genes were genotyped in verification population and 7 of them showed significant association with V. anguillarum resistance in both alleles and genotypes (P < 0.05). Among these 7 genes, PIK3CA-like, CYLD, VCAM1, RhoB and RhoGEF are involved in PI3K/Akt/mTOR pathway and NF-κB pathway, which influence the efficiency of bacteria entering the host and inflammation. SNP-SNP interaction analysis was performed by generalized multifactor dimensionality reduction (GMDR). The combination of SNP loci in RhoB, PIK3CA-like and ADCY3 showed a significant effect on V. anguillarum resistance with the verification rate in the sequencing population up to 70.8%. Taken all, our findings demonstrated the feasibility of BSA-seq approach in identifying genes responsible for the extreme phenotypes and will aid in performing MAS in turbot.

Transcriptomic profiles of striped snakehead cells (SSN-1) infected with snakehead vesiculovirus (SHVV) identifying IFI35 as a positive factor for SHVV replication

Snakehead vesiculovirus (SHVV) has caused great economic loss in snakehead fish culture in China. However, there is no effective strategy to prevent the epidemic of the virus. Understanding the host factors in response to virus infection is the basis for the prevention of viral disease. In this study, the transcriptomic profiles of SHVV-infected and mock-infected SSN-1 cells (derived from striped snakehead, Channa striatus) at 3 and 24 h (h) post of infection (poi) were obtained using high-throughput sequencing technique. A total of 93,372 unigenes were obtained. The differently expressed genes (DEGs) of SSN-1 cells upon SHVV infection were thereby identified, including 3668 and 3536 DEGs at 3 and 24 h poi, respectively. These DEGs were involved in many pathways of viral pathogenesis, including retinoic acid-inducible gene I (RIG-I) like receptors pathway, Toll-like receptor signaling pathway, NF-kappa B signaling pathway, PI3K-Akt signaling pathway and MAPK signaling pathway. Therefore, several immune-related DEGs were randomly selected and confirmed by quantitative real-time PCR (qRT-PCR). In addition, the effects of the interferon inducible protein 35 (IFI35) on SHVV replication were further investigated. Over-expression or inhibition of IFI35 significantly promoted or reduced SHVV replication at the level of viral gene expression, which indicated that IFI35 might be a positive factor for SHVV replication in SSN-1 cells. Our findings presented some valuable information, which will benefit for future study on SHVV-host interactions.

Differentially Expressed Genes Associated with the Cabbage Yellow-Green-Leaf Mutant in the ygl-1 Mapping Interval with Recombination Suppression

Although the genetics and preliminary mapping of the cabbage yellow-green-leaf mutant YL-1 has been extensively studied, transcriptome profiling associated with the yellow-green-leaf mutant of YL-1 has not been discovered. Positional mapping with two populations showed that the yellow-green-leaf gene ygl-1 is located in a recombination-suppressed genomic region. Then, a bulk segregant RNA-seq (BSR) was applied to identify differentially expressed genes (DEGs) using an F₃ population (YL-1 × 11-192) and a BC₂ population (YL-1 × 01-20). Among the 37,286 unique genes, 5730 and 4118 DEGs were detected between the yellow-leaf and normal-leaf pools from the F₃ and BC₂ populations. BSR analysis with four pools greatly reduced the number of common DEGs from 4924 to 1112. In the ygl-1 gene mapping region with suppressed recombination, 43 common DEGs were identified. Five of the DEGs were related to chloroplasts, including the down-regulated Bo1g087310, Bo1g094360, and Bo1g098630 and the up-regulated Bo1g059170 and Bo1g098440. The Bo1g098440 and Bo1g098630 genes were excluded by qRT-PCR. Hence, we inferred that these three DEGs (Bo1g094360, Bo1g087310, and Bo1g059170) in the mapping interval may be tightly associated with the development of the yellow-green-leaf mutant phenotype.

A Review of Molecular Responses of Catfish to Bacterial Diseases and Abiotic Stresses

Catfish is one of the major aquaculture species in the United States. However, the catfish industry is threatened by several bacterial diseases such as enteric septicemia of catfish (ESC), columnaris disease and Aeromonas disease, as well as by abiotic stresses such as high temperature and low oxygen. Research has been conducted for several decades to understand the host responses to these diseases and abiotic stresses. With the development of sequencing technologies, and the application of genome-wide association studies in aquaculture species, significant progress has been made. This review article summarizes recent progress in understanding the molecular responses of catfish after bacterial infection and stress challenges, and in understanding of genomic and genetic basis for disease resistance and stress tolerance.

Genome-wide association analysis of intra-specific QTL associated with the resistance for enteric septicemia of catfish

Disease resistance is one of the most important traits for aquaculture industry. For catfish industry, enteric septicemia of catfish (ESC), caused by the bacterial pathogen Edwardsiella ictaluri, is the most severe disease, causing enormous economic losses every year. In this study, we used three channel catfish families with 900 individuals (300 fish per family) and the 690K catfish SNP array, and conducted a genome-wide association study to detect the quantitative trait loci (QTL) associated with ESC resistance. Three significant QTL, with two of located on LG1 and one on LG26, and three suggestive QTL located on LG1, LG3, and LG21, respectively, were identified to be associated with ESC resistance. With a well-assembled- and -annotated reference genome sequence, genes around the involved QTL regions were identified. Among these genes, 37 genes had known functions in immunity, which may be involved in ESC resistance. Notably, nlrc3 and nlrp12 identified here were also found in QTL regions of ESC resistance in the channel catfish × blue catfish interspecific hybrid system, suggesting this QTL was operating within both intra-specific channel catfish populations and interspecific hybrid backcross populations. Many of the genes of the Class I MHC pathway, for mediated antigen processing and presentation, were found in the QTL regions. The positional correlation found in this study and the expressional correlation found in previous studies indicated that Class I MHC pathway was significantly associated with ESC resistance. This study validated one QTL previously identified using the second and fourth generation of the interspecific hybrid backcross progenies, and identified five additional QTL among channel catfish families. Taken together, it appears that there are only a few major QTL for ESC disease resistance, making marker-assisted selection an effective approach for genetic improvements of ESC resistance.

Applications of transcriptomics and proteomics in understanding fish immunity

With the development of intensive aquaculture, economic losses increasingly result from fish mortality due to pathogen infection. In recent years, a growing number of researchers have used transcriptomic and proteomic analyses to study fish immune responses to exogenous pathogen infection. Integrating transcriptomic and proteomic analyses provides a better understanding of the fish immune system including gene expression, regulation, and the intricate biological processes underlying immune responses against infection. This review focuses on the recent advances in the fields of transcriptomics and proteomics, which have contributed to our understanding of fish immunity to exogenous pathogens.

GWAS analysis using interspecific backcross progenies reveals superior blue catfish alleles responsible for strong resistance against enteric septicemia of catfish

Infectious diseases pose significant threats to the catfish industry. Enteric septicemia of catfish (ESC) caused by Edwardsiella ictaluri is the most devastating disease for catfish aquaculture, causing huge economic losses annually. Channel catfish and blue catfish exhibit great contrast in resistance against ESC, with channel catfish being highly susceptible and blue catfish being highly resistant. As such, the interspecific backcross progenies provide an ideal system for the identification of quantitative trait locus (QTL). We previously reported one significant QTL on linkage group (LG) 1 using the third-generation backcrosses, but the number of founders used to make the second- and third-generation backcross progenies was very small. Although the third-generation backcross progenies provided a greater power for fine mapping than the first-generation backcrosses, some major QTL for disease resistance may have been missing due to the small numbers of founders used to produce the higher generation backcrosses. In this study, we performed a genome-wide association study using first-generation backcrosses with the catfish 690 K SNP arrays to identify additional ESC disease resistance QTL, especially those at the species level. Two genomic regions on LG1 and LG23 were determined to be significantly associated with ESC resistance as revealed by a mixed linear model and family-based association test. Examination of the resistance alleles indicated their origin from blue catfish, indicating that at least two major disease resistance loci exist among blue catfish populations. Upon further validation, markers linked with major ESC disease resistance QTL should be useful for marker-assisted introgression, allowing development of highly ESC resistant breeds of catfish.

JAK and STAT members in channel catfish: Identification, phylogenetic analysis and expression profiling after Edwardsiella ictaluri infection

The Janus kinase/signal transducers and activators of transcription (JAK/STAT) signaling pathway is one of the main pleiotropic cascades used to transmit information from extracellular receptors to the nucleus, which results in DNA transcription and expression of genes involved in immunity, proliferation, differentiation, migration, apoptosis, and cell survival. Members of JAK family and STAT family have been extensively studied in different mammalian species because of their important roles in innate and adaptive immune responses. However, they have not been systematically studied among teleost fish species. In this study, five JAK family members and eight STAT family members were identified and characterized from channel catfish. Phylogenetic analysis was conducted to properly annotate these genes. Syntenic analysis was also conducted to establish orthology, and confirm the results from phylogenetic analysis. Compared to mammals, more members of the JAK and STAT family were identified in channel catfish genome. Expression of JAK and STAT family members was detected in healthy catfish tissues, but was induced in gill, liver, and intestine after bacterial challenge. Notably, the significant upregulation of STAT1b gene in catfish liver, gill and intestine after Edwardsiella ictaluri infection supported the notion that high STAT1 expression are involved in defense against pathogens. Collectively, the increased expression of JAK and STAT members in tested tissues suggested their crucial function in defending the host against pathogen invasion.

Transcriptome Analysis Based on RNA-Seq in Understanding Pathogenic Mechanisms of Diseases and the Immune System of Fish: A Comprehensive Review

In recent years, with the advent of next-generation sequencing along with the development of various bioinformatics tools, RNA sequencing (RNA-Seq)-based transcriptome analysis has become much more affordable in the field of biological research. This technique has even opened up avenues to explore the transcriptome of non-model organisms for which a reference genome is not available. This has made fish health researchers march towards this technology to understand pathogenic processes and immune reactions in fish during the event of infection. Recent studies using this technology have altered and updated the previous understanding of many diseases in fish. RNA-Seq has been employed in the understanding of fish pathogens like bacteria, virus, parasites, and oomycetes. Also, it has been helpful in unraveling the immune mechanisms in fish. Additionally, RNA-Seq technology has made its way for future works, such as genetic linkage mapping, quantitative trait analysis, disease-resistant strain or broodstock selection, and the development of effective vaccines and therapies. Until now, there are no reviews that comprehensively summarize the studies which made use of RNA-Seq to explore the mechanisms of infection of pathogens and the defense strategies of fish hosts. This review aims to summarize the contemporary understanding and findings with regard to infectious pathogens and the immune system of fish that have been achieved through RNA-Seq technology.

The CC and CXC chemokine receptors in channel catfish (Ictalurus punctatus) and their involvement in disease and hypoxia responses

Chemokines are vital regulators of cell mobilization for immune surveillance, inflammation, and development. Chemokines signal through binding to their receptors that are a superfamily of seven-transmembrane domain G-coupled receptors. Recently, a complete repertoire of both CC and CXC chemokines have been identified in channel catfish, but nothing is known about their receptors. In this study, a set of 29 CC chemokine receptor (CCR) genes and 8 CXC chemokine receptor (CXCR) genes were identified and annotated from the channel catfish genome. Extensive phylogenetic and comparative genomic analyses were conducted to annotate these genes, revealing fish-specific CC chemokine receptors, and lineage-specific tandem duplications of chemokine receptors in the teleost genomes. With 29 genes, the channel catfish genome harbors the largest numbers of CC chemokine receptors among all the genomes characterized. Analysis of gene expression after bacterial infections indicated that the chemokine receptors were regulated in a gene-specific manner. Most differentially expressed chemokine receptors were up-regulated after Edwardsiella ictaluri and Flavobacterium columnare infection. Among which, CXCR3 and CXCR4 were observed to participate in immune responses to both bacterial infections, indicating their potential roles in catfish immune activities. In addition, CXCR3.2 was significantly up-regulated in ESC-susceptible fish, and CXCR4b was mildly induced in ESC-resistant fish, further supporting the significant roles of CXCR3 and CXCR4 in catfish immune responses. CXCR4b and CCR9a were both up-regulated not only after bacterial infection, but also after hypoxia stress, providing the linkage between bacterial infection and low oxygen stresses. These results should be valuable for comparative immunological studies and provide insights into their roles in disease and stress responses.

Identification of SNPs associated with muscle yield and quality traits using allelic-imbalance analyses of pooled RNA-Seq samples in rainbow trout

Background

Coding/functional SNPs change the biological function of a gene and, therefore, could serve as “large-effect” genetic markers. In this study, we used two bioinformatics pipelines, GATK and SAMtools, for discovering coding/functional SNPs with allelic-imbalances associated with total body weight, muscle yield, muscle fat content, shear force, and whiteness. Phenotypic data were collected for approximately 500 fish, representing 98 families (5 fish/family), from a growth-selected line, and the muscle transcriptome was sequenced from 22 families with divergent phenotypes (4 low- versus 4 high-ranked families per trait).

Results

GATK detected 59,112 putative SNPs; of these SNPs, 4798 showed allelic imbalances (>2.0 as an amplification and <0.5 as loss of heterozygosity). SAMtools detected 87,066 putative SNPs; and of them, 4962 had allelic imbalances between the low- and high-ranked families. Only 1829 SNPs with allelic imbalances were common between the two datasets, indicating significant differences in algorithms. The two datasets contained 7930 non-redundant SNPs of which 4439 mapped to 1498 protein-coding genes (with 6.4% non-synonymous SNPs) and 684 mapped to 295 lncRNAs. Validation of a subset of 92 SNPs revealed 1) 86.7–93.8% success rate in calling polymorphic SNPs and 2) 95.4% consistent matching between DNA and cDNA genotypes indicating a high rate of identifying SNPs with allelic imbalances. In addition, 4.64% SNPs revealed random monoallelic expression. Genome distribution of the SNPs with allelic imbalances exhibited high density for all five traits in several chromosomes, especially chromosome 9, 20 and 28. Most of the SNP-harboring genes were assigned to important growth-related metabolic pathways.

Conclusion

These results demonstrate utility of RNA-Seq in assessing phenotype-associated allelic imbalances in pooled RNA-Seq samples. The SNPs identified in this study were included in a new SNP-Chip design (available from Affymetrix) for genomic and genetic analyses in rainbow trout.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3992-z) contains supplementary material, which is available to authorized users.

The NCK and ABI adaptor genes in catfish and their involvement in ESC disease response

Adaptor proteins non-catalytic region of tyrosine kinase (NCK) and Abelson interactor (ABI) are crucial for disease response. NCK1 was identified to be a candidate gene for enteric septicemia of catfish (ESC) disease resistance, and was speculated to play similar roles during ESC and enteropathogenic Escherichia coli (EPEC) pathogenicity. ABI1 was reported as a positional candidate gene for bacterial cold water disease (BCWD) resistance in rainbow trout. In this study, three NCK genes and six ABI genes were identified in the channel catfish (Ictalurus punctatus) genome and blue catfish (I. furcatus) transcriptome, and annotated by domain structures, phylogenetic and syntenic analyses. Their expression patterns were examined in the intestine and liver of catfish after challenge with Edwardsiella ictaluri. In the intestine, NCK1, ABI2a, ABI2b, ABI3a were differentially expressed after E. ictaluri infection. In the liver, NCK2a, NCK2b, ABI1b, ABI2a, ABI2b were significantly upregulated in ESC susceptible fish. In general, the NCK and ABI genes, with exception of ABI3a gene and NCK1 gene, were expressed at higher levels in susceptible fish after infection than in control fish, but were expressed at lower levels in resistant fish than in the control fish. Taken together, these results support the notion that NCK and ABI genes are involved in disease processes facilitating pathogenesis of the E. ictaluri bacteria.

The chemokinome superfamily: II. The 64 CC chemokines in channel catfish and their involvement in disease and hypoxia responses

Chemokines are a superfamily of structurally related chemotactic cytokines exerting significant roles in regulating cell migration and activation. Based on the arrangement of the first four cysteine residues, they are classified into CC, CXC, C and CX3C subfamilies. In this study, a complete set of 64 CC chemokine ligand (CCL) genes was systematically identified, annotated, and characterized from the channel catfish genome. Extensive phylogenetic and comparative genomic analyses supported their annotations, allowing establishment of their orthologies, revealing fish-specific CC chemokines and the expansion of CC chemokines in the teleost genomes through lineage-specific tandem duplications. With 64 genes, the channel catfish genome harbors the largest numbers of CC chemokines among all the genomes characterized to date, however, they fall into 11 distinct CC chemokine groups. Analysis of gene expression after bacterial infections indicated that the CC chemokines were regulated in a gene-specific and time-dependent manner. While only one member of CCL19 (CCL19a.1) was significantly up-regulated after Edwardsiella ictaluri infection, all CCL19 members (CCL19a.1, CCL19a.2 and CCL19b) were significantly induced after Flavobacterium columnare infection. In addition, CCL19a.1, CCL19a.2 and CCL19b were also drastically up-regulated in ESC-susceptible fish, but not in resistant fish, suggesting potential significant roles of CCL19 in catfish immune responses. High expression levels of certain CC appeared to be correlated with susceptibility to diseases and intolerance to hypoxia.

Bulked Segregant RNA-seq Reveals Differential Expression and SNPs of Candidate Genes Associated with Waterlogging Tolerance in Maize

Waterlogging has increasingly become one of the major constraints to maize productivity in some maize production zones because it causes serious yield loss. Bulked segregant RNA-seq (BSR-seq) has been widely applied to profile candidate genes and map associated Single Nucleotide Polymorphism (SNP) markers in many species. In this study, 10 waterlogging sensitive and eight tolerant inbred lines were selected from 60 maize inbred lines with waterlogging response determined and preselected by the International Maize and Wheat Improvement Center (CIMMYT) from over 400 tropical maize inbred lines. BSR-seq was performed to identify differentially expressed genes and SNPs associated with waterlogging tolerance. Upon waterlogging stress, 354 and 1094 genes were differentially expressed in the tolerant and sensitive pools, respectively, compared to untreated controls. When tolerant and sensitive pools were compared, 593 genes were differentially expressed under untreated and 431 genes under waterlogged conditions, of which 122 genes overlapped. To validate the BSR-seq results, the expression levels of six genes were determined by qRT-PCR. The qRT-PCR results were consistent with BSR-seq results. Comparison of allelic polymorphism in mRNA sequences between tolerant and sensitive pools revealed 165 (normal condition) and 128 (waterlogged condition) high-probability SNPs. We found 18 overlapping SNPs with genomic positions mapped. Eighteen SNPs were contained in 18 genes, and eight and nine of 18 genes were responsive to waterlogging stress in tolerant and sensitive lines, respectively. Six alleles of the 18 originated from tolerant pool were significantly up-regulated under waterlogging, but not those from sensitive pool. Importantly, one allele (GRMZM2G055704) of the six genes was mapped between umc1619 and umc1948 on chromosome 1 where a QTL associated with waterlogging tolerance was identified in a previous research, strongly indicating that GRMZM2G055704 is a candidate gene responsive to waterlogging. Our research contributes to the knowledge of the molecular mechanism for waterlogging tolerance in maize.

Functional natural allelic variants of flavonoid 3',5'-hydroxylase gene governing catechin traits in tea plant and its relatives

Functional allelic variants of the flavonoid 3',5'-hydroxylase (F3'5'H) gene provides new information of F3'5'H function of tea plant and its relatives. This insight may serve as the foundation upon which to advance molecular breeding in the tea plant. Catechins are the active components of tea that determine its quality and health attributes. This study established the first integrated genomic strategy for deciphering the genetic basis of catechin traits of tea plant. With the RNA-sequencing analysis of bulked segregants representing the tails of a F1 population segregated for total catechin content, we identified a flavonoid 3',5'-hydroxylase (F3'5'H) gene. F3'5'H had one copy in the genomic DNA of tea plant. Among 202 tea accessions, we identified 120 single nucleotide polymorphisms (SNPs) at F3'5'H locus. Seventeen significant marker-trait associations were identified by association mapping in multiple environments, which were involved in 10 SNP markers, and the traits including the ratio of di/tri-hydroxylated catechins and catechin contents. The associated individual and combination of SNPs explained 4.5-25.2 and 53.0-63.0% phenotypic variations, respectively. In the F1 population (validation population), the catechin trait variation percentages explained by F3'5'H diplotype were 6.9-74.3%. The genotype effects of ten functional SNPs in the F1 population were all consistent with the association population. Furthermore, the function of SNP-711/-655 within F3'5'H was validated by gene expression analysis. Altogether, our work indicated functional SNP allelic variants within F3'5'H governing the ratio of di/tri-hydroxylated catechins and catechin contents. The strong catechin-associated SNPs identified in this study can be used for future marker-assisted selection to improve tea quality.

Aquaculture genomics, genetics and breeding in the United States: current status, challenges, and priorities for future research

Advancing the production efficiency and profitability of aquaculture is dependent upon the ability to utilize a diverse array of genetic resources. The ultimate goals of aquaculture genomics, genetics and breeding research are to enhance aquaculture production efficiency, sustainability, product quality, and profitability in support of the commercial sector and for the benefit of consumers. In order to achieve these goals, it is important to understand the genomic structure and organization of aquaculture species, and their genomic and phenomic variations, as well as the genetic basis of traits and their interrelationships. In addition, it is also important to understand the mechanisms of regulation and evolutionary conservation at the levels of genome, transcriptome, proteome, epigenome, and systems biology. With genomic information and information between the genomes and phenomes, technologies for marker/causal mutation-assisted selection, genome selection, and genome editing can be developed for applications in aquaculture. A set of genomic tools and resources must be made available including reference genome sequences and their annotations (including coding and non-coding regulatory elements), genome-wide polymorphic markers, efficient genotyping platforms, high-density and high-resolution linkage maps, and transcriptome resources including non-coding transcripts. Genomic and genetic control of important performance and production traits, such as disease resistance, feed conversion efficiency, growth rate, processing yield, behaviour, reproductive characteristics, and tolerance to environmental stressors like low dissolved oxygen, high or low water temperature and salinity, must be understood. QTL need to be identified, validated across strains, lines and populations, and their mechanisms of control understood. Causal gene(s) need to be identified. Genetic and epigenetic regulation of important aquaculture traits need to be determined, and technologies for marker-assisted selection, causal gene/mutation-assisted selection, genome selection, and genome editing using CRISPR and other technologies must be developed, demonstrated with applicability, and application to aquaculture industries.Major progress has been made in aquaculture genomics for dozens of fish and shellfish species including the development of genetic linkage maps, physical maps, microarrays, single nucleotide polymorphism (SNP) arrays, transcriptome databases and various stages of genome reference sequences. This paper provides a general review of the current status, challenges and future research needs of aquaculture genomics, genetics, and breeding, with a focus on major aquaculture species in the United States: catfish, rainbow trout, Atlantic salmon, tilapia, striped bass, oysters, and shrimp. While the overall research priorities and the practical goals are similar across various aquaculture species, the current status in each species should dictate the next priority areas within the species. This paper is an output of the USDA Workshop for Aquaculture Genomics, Genetics, and Breeding held in late March 2016 in Auburn, Alabama, with participants from all parts of the United States.

Identification, annotation and expression analysis of 29 Rho GTPase genes from channel catfish (Ictalurus punctatus) after bacterial infections

The Rho family GTPases are a group of small monomeric G proteins, which are molecular switches in signaling pathways. They have been known to regulate a diverse range of cellular processes including actin cytoskeleton rearrangement and microtubule dynamics. In particular, their participations in immune responses are also significant. However, little information of the Rho GTPases is available in teleost including channel catfish, an economically important species and one of the best teleost models forimmunological research. In this study, Rho GTPase genes were identified from channel catfish and well annotated by phylogenetic and syntenic analyses. Their expression profiles were determined in channel catfish healthy tissues and infected tissues. Altogether seven Rho GTPase genes were significantly regulated after bacterial infection, with six genes in the gill after Flavobacterium columnare challenge and two genes in the intestine in response to Edwardsiella ictaluri. All the differentially expressed genes were up-regulated soon after bacterial infection. Different expression patterns between the two experiments were observed, which may be attributed to tissue-specific regulation or pathogen-specific regulation. These results suggested that Rho GTPases play important roles in immune responses to bacterial pathogens, setting a foundation for future investigation on Rho GTPases.

Development of a 690 K SNP array in catfish and its application for genetic mapping and validation of the reference genome sequence

Single nucleotide polymorphisms (SNPs) are capable of providing the highest level of genome coverage for genomic and genetic analysis because of their abundance and relatively even distribution in the genome. Such a capacity, however, cannot be achieved without an efficient genotyping platform such as SNP arrays. In this work, we developed a high-density SNP array with 690,662 unique SNPs (herein 690 K array) that were relatively evenly distributed across the entire genome, and covered 98.6% of the reference genome sequence. Here we also report linkage mapping using the 690 K array, which allowed mapping of over 250,000 SNPs on the linkage map, the highest marker density among all the constructed linkage maps. These markers were mapped to 29 linkage groups (LGs) with 30,591 unique marker positions. This linkage map anchored 1,602 scaffolds of the reference genome sequence to LGs, accounting for over 97% of the total genome assembly. A total of 1,007 previously unmapped scaffolds were placed to LGs, allowing validation and in few instances correction of the reference genome sequence assembly. This linkage map should serve as a valuable resource for various genetic and genomic analyses, especially for GWAS and QTL mapping for genes associated with economically important traits.

The chemokinome superfamily in channel catfish: I. CXC subfamily and their involvement in disease defense and hypoxia responses

Chemokines are a superfamily of structurally related chemotactic cytokines exerting significant roles in regulating cell migration and activation. They are defined by the presence of four conserved cysteine residues and are divided into four subfamilies depending on the arrangement of the first two conserved cysteines residues: CXC, CC, C and CX3C. In this study, a complete set of 17 CXC chemokine ligand (CXCL) genes was systematically identified and characterized from channel catfish genome through data mining of existing genomic resources. Phylogenetic analysis allowed annotation of the 17 CXC chemokines. Extensive comparative genomic analyses supported their annotations and orthologies, revealing the existence of fish-specific CXC chemokines and the expansion of CXC chemokines in the teleost genomes. The analysis of gene expression after bacterial infection indicated the CXC chemokines were expressed in a gene-specific manner. CXCL11.3 and CXCL20.3 were expressed significantly higher in resistant fish than in susceptible fish after ESC infection, while CXCL20.2 were expressed significantly higher in resistant fish than in susceptible fish after columnaris infection. The expression of those CXC chemokines, therefore can be a useful indicator of disease resistance. A similar pattern of expression was observed between resistant and susceptible fish with biotic and abiotic stresses, ESC, columnaris and hypoxia, suggesting that high levels of expression of the majority of CXC chemokines, with exception of CXC11 and CXC20, are detrimental to the host.

Genetic architecture of the maize kernel row number revealed by combining QTL mapping using a high-density genetic map and bulked segregant RNA sequencing

BACKGROUND

The maize kernel row number (KRN) is a key component that contributes to grain yield and has high broad-sense heritability (H ² ). Quantitative trait locus/loci (QTL) mapping using a high-density genetic map is a powerful approach to detecting loci that are responsible for traits of interest. Bulked segregant ribonucleic acid (RNA) sequencing (BSR-seq) is another rapid and cost-effective strategy to identify QTL. Combining QTL mapping using a high-density genetic map and BSR-seq may dissect comprehensively the genetic architecture underlying the maize KRN.

RESULTS

A panel of 300 F₂ individuals derived from inbred lines abe2 and B73 were genotyped using the specific-locus amplified fragment sequencing (SLAF-seq) method. A total of 4,579 high-quality polymorphic SLAF markers were obtained and used to construct a high-density genetic map with a total length of 2,123 centimorgan (cM) and an average distance between adjacent markers of 0.46 cM. Combining the genetic map and KRN of F₂ individuals, four QTL (qKRN1, qKRN2, qKRN5, and qKRN8-1) were identified on chromosomes 1, 2, 5, and 8, respectively. The physical intervals of these four QTL ranged from 4.36 Mb for qKRN8-1 to 7.11 Mb for qKRN1 with an average value of 6.08 Mb. Based on high-throughput sequencing of two RNA pools bulked from leaves of plants with extremely high and low KRNs, two QTL were detected on chromosome 8 in the 10-25 Mb (BSR_QTL1) and 60-150 Mb (BSR_QTL2) intervals. According to the physical positions of these QTL, qKRN8-1 was included by BSR_QTL2. In addition, qKRN8-1 was validated using QTL mapping with a recombinant inbred lines population that was derived from inbred lines abe2 and B73.

CONCLUSIONS

In this study, we proved that combining QTL mapping using a high-density genetic map and BSR-seq is a powerful and cost-effective approach to comprehensively revealing genetic architecture underlying traits of interest. The QTL for the KRN detected in this study, especially qKRN8-1, can be used for performing fine mapping experiments and marker-assisted selection in maize breeding.

Gene expression comparison of resistant and susceptible Atlantic salmon fry challenged with Infectious Pancreatic Necrosis virus reveals a marked contrast in immune response

BACKGROUND

Infectious Pancreatic Necrosis (IPN) is a highly contagious birnavirus disease of farmed salmonid fish, which often causes high levels of morbidity and mortality. A large host genetic component to resistance has been previously described for Atlantic salmon (Salmo salar L.), which mediates high mortality rates in some families and zero mortality in others. However, the molecular and immunological basis for this resistance is not yet fully known. This manuscript describes a global comparison of the gene expression profiles of resistant and susceptible Atlantic salmon fry following challenge with the IPN virus.

RESULTS

Salmon fry from two IPNV-resistant and two IPNV-susceptible full sibling families were challenged with the virus and sampled at 1 day, 7 days and 20 days post-challenge. Significant viral titre was observed in both resistant and susceptible fish at all timepoints, although generally at higher levels in susceptible fish. Gene expression profiles combined with gene ontology and pathway analyses demonstrated that while a clear immune response was observed in both resistant and susceptible fish, there were striking differences between the two phenotypes. The susceptible fish showed marked up-regulation of genes related to cytokine activity and inflammatory response that evidently failed to protect against the virus. In contrast, the resistant fish demonstrated a less pronounced immune response including up-regulation of genes relating to the M2 macrophage system.

CONCLUSIONS

While only the susceptible phenotype shows appreciable mortality levels, both resistant and susceptible fish can become infected with IPNV. Susceptible fish are characterized by a much larger, yet ineffective, immune response, largely related to cytokine and inflammatory systems. Resistant fish demonstrate a more moderate, putative macrophage-mediated inflammatory response, which may contribute to their survival.

Identification of Immune-Related Genes and Development of SSR/SNP Markers from the Spleen Transcriptome of Schizothorax prenanti

Schizothorax prenanti (S. prenanti) is mainly distributed in the upstream regions of the Yangtze River and its tributaries in China. This species is indigenous and commercially important. However, in recent years, wild populations and aquacultures have faced the serious challenges of germplasm variation loss and an increased susceptibility to a range of pathogens. Currently, the genetics and immune mechanisms of S. prenanti are unknown, partly due to a lack of genome and transcriptome information. Here, we sought to identify genes related to immune functions and to identify molecular markers to study the function of these genes and for trait mapping. To this end, the transcriptome from spleen tissues of S. prenanti was analyzed and sequenced. Using paired-end reads from the Illumina Hiseq2500 platform, 48,517 transcripts were isolated from the spleen transcriptome. These transcripts could be clustered into 37,785 unigenes with an N50 length of 2,539 bp. The majority of the unigenes (35,653, 94.4%) were successfully annotated using non-redundant nucleotide sequence analysis (nt), and the non-redundant protein (nr), Swiss-Prot, Gene Ontology (GO), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. KEGG pathway assignment identified more than 500 immune-related genes. Furthermore, 7,545 putative simple sequence repeats (SSRs), 857,535 single nucleotide polymorphisms (SNPs), and 53,481 insertion/deletion (InDels) were detected from the transcriptome. This is the first reported high-throughput transcriptome analysis of S. prenanti, and it provides valuable genetic resources for the investigation of immune mechanisms, conservation of germplasm, and molecular marker-assisted breeding of S. prenanti.

Identification of genes regulated by histone acetylation during root development in Populus trichocarpa

Background

Histone deacetylases (HDACs) are key enzymes catalyzing the removal of acetyl groups from histones. HDACs act in concert with histone acetyltransferases (HATs) to regulate histone acetylation status, which modifies chromatin structure, affecting gene transcription and thus regulating multiple biological processes such as plant growth and development. Over a decade, certain HDACs in herbaceous plants have been deeply studied. However, functions of HDACs in woody plants are not well understood.

Results

Histone deacetylase specific inhibitor trichostatin A (TSA) was used to investigate the role of HDACs in organogenesis of roots and root development in Populus trochocarpa. The adventitious roots were regenerated and grown on medium supplemented with 0, 1, and 2.5 μM TSA. TSA treatment delayed root regeneration and inhibited primary root growth. To examine the genes modified by TSA in the regenerated roots, tag-based digital gene expression (DGE) analysis was performed using Illumina HiSeqTM 2000. Approximately 4.5 million total clean tags were mapped per library. The distinct clean tags for the three libraries corresponding to 0, 1 and 2.5 μM TSA treatment were 166167, 143103 and 153507, from which 38.45 %, 31.84 % and 38.88 % were mapped unambiguously to the unigene database, respectively. Most of the tags were expressed at similar levels, showing a < 5-fold difference after 1 μM and 2.5 μM TSA treatments and the maximum fold-change of the tag copy number was around 20. The expression levels of many genes in roots were significantly altered by TSA. A total of 36 genes were up-regulated and 1368 genes were down-regulated after 1 μM TSA treatment, while 166 genes were up-regulated and 397 genes were down-regulated after 2.5 μM TSA treatment. Gene ontology (GO) and pathway analyses indicated that the differentially expressed genes were related to many kinds of molecular functions and biological processes. The genes encoding key enzymes catalyzing gibberellin biosynthesis were significantly down-regulated in the roots exposed to 2.5 μM TSA and their expression changes were validated by using real-time PCR.

Conclusions

HDACs were required for de novo organogenesis and normal growth of populus roots. DGE data provides the gene profiles in roots probably regulated by histone acetylation during root growth and development, which will lead to a better understanding of the mechanism controlling root development.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-2407-x) contains supplementary material, which is available to authorized users.

Evaluation of potential candidate genes involved in salinity tolerance in striped catfish (Pangasianodon hypophthalmus) using an RNA-Seq approach

Increasing salinity levels in freshwater and coastal environments caused by sea level rise linked to climate change is now recognized to be a major factor that can impact fish growth negatively, especially for freshwater teleost species. Striped catfish (Pangasianodon hypophthalmus) is an important freshwater teleost that is now widely farmed across the Mekong River Delta in Vietnam. Understanding the basis for tolerance and adaptation to raised environmental salinity conditions can assist the regional culture industry to mitigate predicted impacts of climate change across this region. Attempt of next generation sequencing using the ion proton platform results in more than 174 million raw reads from three tissue libraries (gill, kidney and intestine). Reads were filtered and de novo assembled using a variety of assemblers and then clustered together to generate a combined reference transcriptome. Downstream analysis resulted in a final reference transcriptome that contained 60,585 transcripts with an N50 of 683 bp. This resource was further annotated using a variety of bioinformatics databases, followed by differential gene expression analysis that resulted in 3062 transcripts that were differentially expressed in catfish samples raised under two experimental conditions (0 and 15 ppt). A number of transcripts with a potential role in salinity tolerance were then classified into six different functional gene categories based on their gene ontology assignments. These included; energy metabolism, ion transportation, detoxification, signal transduction, structural organization and detoxification. Finally, we combined the data on functional salinity tolerance genes into a hypothetical schematic model that attempted to describe potential relationships and interactions among target genes to explain the molecular pathways that control adaptive salinity responses in P. hypophthalmus. Our results indicate that P. hypophthalmus exhibit predictable plastic regulatory responses to elevated salinity by means of characteristic gene expression patterns, providing numerous candidate genes for future investigations.

Patterns of Piscirickettsia salmonis load in susceptible and resistant families of Salmo salar

The pathogen Piscirickettsia salmonis produces a systemic aggressive infection that involves several organs and tissues in salmonids. In spite of the great economic losses caused by this pathogen in the Atlantic salmon (Salmo salar) industry, very little is known about the resistance mechanisms of the host to this pathogen. In this paper, for the first time, we aimed to identify the bacterial load in head kidney and muscle of Atlantic salmon exhibiting differential familiar mortality. Furthermore, in order to assess the patterns of gene expression of immune related genes in susceptible and resistant families, a set of candidate genes was evaluated using deep sequencing of the transcriptome. The results showed that the bacterial load was significantly lower in resistant fish, when compared with the susceptible individuals. Based on the candidate genes analysis, we infer that the resistant hosts triggered up-regulation of specific genes (such as for example the LysC), which may explain a decrease in the bacterial load in head kidney, while the susceptible fish presented an exacerbated innate response, which is unable to exert an effective response against the bacteria. Interestingly, we found a higher bacterial load in muscle when compared with head kidney. We argue that this is possible due to the availability of an additional source of iron in muscle. Besides, the results show that the resistant fish could not be a likely reservoir of the bacteria.