Fishing Innate Immune System Properties through the Transcriptomic Single-Cell Data of Teleostei

The innate immune system is the first line of defense in multicellular organisms. Danio rerio is widely considered a promising model for IIS-related research, with the most amount of scRNAseq data available among Teleostei. We summarized the scRNAseq and spatial transcriptomics experiments related to the IIS for zebrafish and other Teleostei from the GEO NCBI and the Single-Cell Expression Atlas. We found a considerable number of scRNAseq experiments at different stages of zebrafish development in organs such as the kidney, liver, stomach, heart, and brain. These datasets could be further used to conduct large-scale meta-analyses and to compare the IIS of zebrafish with the mammalian one. However, only a small number of scRNAseq datasets are available for other fish (turbot, salmon, cavefish, and dark sleeper). Since fish biology is very diverse, it would be a major mistake to use zebrafish alone in fish immunology studies. In particular, there is a special need for new scRNAseq experiments involving nonmodel Teleostei, e.g., long-lived species, cancer-resistant fish, and various fish ecotypes.

Isolation and transcriptome analysis of a biotechnologically promising Black Sea protist, Thraustochytrium aureum ssp. strugatskii

Background

Marine protists are an important part of the ocean ecosystem. They may possess unique sets of biosynthetic pathways and, thus, be promising model organisms for metabolic engineering for producing substances for the pharmaceutical, cosmetic, and perfume industries. Currently, full-genome data are available just for a limited number of protists hampering their use in biotechnology.

Methods

We characterized the morphology of a new cultured strain of Thraustochytriaceae isolated from the Black Sea ctenophore Beroe ovata using phase-contrast microscopy. Cell culture was performed in the FAND culture medium based on fetal bovine serum and DMEM. Phylogenetic analysis was performed using the 18S rRNA sequence. We also conducted a transcriptome assembly and compared the data with the closest species.

Results

The protist belongs to the genus Thraustochytrium based on the 18S rRNA sequence analysis. We designated the isolated protist as T. aureum ssp. strugatskii. The closest species with the genome assembly is Schizochytrium aggregatum. Transcriptome analysis revealed the majority of the fatty acid synthesis enzymes.

Conclusion

Our findings suggest that the T. aureum ssp. strugatskii is a promising candidate for biotechnological use. Together with the previously available, our data would allow the establishment of an accurate phylogeny of the family Thraustochytriaceae. Also, it could be a reference point for studying the evolution of the enzyme families.

Comparative transcriptome profiling of a resistant vs susceptible bread wheat (Triticum aestivum L.) cultivar in response to water deficit and cold stress

Bread wheat (Triticum aestivum L.) is one of the most important agricultural plants wearing abiotic stresses, such as water deficit and cold, that cause its productivity reduction. Since resistance to abiotic factors is a multigenic trait, therefore modern genome-wide approaches can help to involve various genetic material in breeding. One technique is full transcriptome analysis that reveals groups of stress response genes serving marker-assisted selection markers. Comparing transcriptome profiles of the same genetic material under several stresses is essential and makes the whole picture. Here, we addressed this by studying the transcriptomic response to water deficit and cold stress for two evolutionarily distant bread wheat varieties: stress-resistant cv. Saratovskaya 29 (S29) and stress-sensitive cv. Yanetzkis Probat (YP). For the first time, transcriptomes for these cultivars grown under abiotic stress conditions were obtained using Illumina based MACE technology. We identified groups of genes involved in response to cold and water deficiency stresses, including responses to each stress factor and both factors simultaneously that may be candidates for resistance genes. We discovered a core group of genes that have a similar pattern of stress-induced expression changes. The particular expression pattern was revealed not only for the studied varieties but also for the published transcriptomic data on cv. Jing 411 and cv. Fielder. Comparative transcriptome profiling of cv. S29 and cv. YP in response to water deficit and cold stress confirmed the hypothesis that stress-induced expression change is unequal within a homeologous gene group. As a rule, at least one changed significantly while the others had a relatively lower expression. Also, we found several SNPs distributed throughout the genomes of cv. S29 and cv. YP and distinguished the studied varieties from each other and the reference cv. Chinese Spring. Our results provide new data for genomics-assisted breeding of stress-tolerant wheat cultivars.

Comparative Expression Analysis of Stress-Inducible Candidate Genes in Response to Cold and Drought in Tea Plant [Camellia sinensis (L.) Kuntze]

Cold and drought are two of the most severe threats affecting the growth and productivity of the tea plant, limiting its global spread. Both stresses cause osmotic changes in the cells of the tea plant by decreasing their water potential. To develop cultivars that are tolerant to both stresses, it is essential to understand the genetic responses of tea plant to these two stresses, particularly in terms of the genes involved. In this study, we combined literature data with interspecific transcriptomic analyses (using Arabidopsis thaliana and Solanum lycopersicum) to choose genes related to cold tolerance. We identified 45 stress-inducible candidate genes associated with cold and drought responses in tea plants based on a comprehensive homologous detection method. Of these, nine were newly characterized by us, and 36 had previously been reported. The gene network analysis revealed upregulated expression in ICE1-related cluster of bHLH factors, HSP70/BAM5 connected genes (hexokinases, galactinol synthases, SnRK complex, etc.) indicating their possible co-expression. Using qRT-PCR we revealed that 10 genes were significantly upregulated in response to both cold and drought in tea plant: HSP70, GST, SUS1, DHN1, BMY5, bHLH102, GR-RBP3, ICE1, GOLS1, and GOLS3. SnRK1.2, HXK1/2, bHLH7/43/79/93 were specifically upregulated in cold, while RHL41, CAU1, Hydrolase22 were specifically upregulated in drought. Interestingly, the expression of CIP was higher in the recovery stage of both stresses, indicating its potentially important role in plant recovery after stress. In addition, some genes, such as DHN3, bHLH79, PEI54, SnRK1.2, SnRK1.3, and Hydrolase22, were significantly positively correlated between the cold and drought responses. CBF1, GOLS1, HXK2, and HXK3, by contrast, showed significantly negative correlations between the cold and drought responses. Our results provide valuable information and robust candidate genes for future functional analyses intended to improve the stress tolerance of the tea plant and other species.

LSM-W2: laser scanning microscopy worker for wheat leaf surface morphology

Background

Microscopic images are widely used in plant biology as an essential source of information on morphometric characteristics of the cells and the topological characteristics of cellular tissue pattern due to modern computer vision algorithms. High-resolution 3D confocal images allow extracting quantitative characteristics describing the cell structure of leaf epidermis. For some issues in the study of cereal leaves development, it is required to apply the staining techniques with fluorescent dyes and to scan rather large fragments consisting of several frames. We aimed to develop a tool for processing multi-frame multi-channel 3D images obtained from confocal laser scanning microscopy and taking into account the peculiarities of the cereal leaves staining.

Results

We elaborated an ImageJ-plugin LSM-W² that allows extracting data on Leaf Surface Morphology from Laser Scanning Microscopy images. The plugin is a crucial link in a workflow for obtaining data on structural properties of leaf epidermis and morphological properties of epidermal cells. It allows converting large lsm-files (laser scanning microscopy) into segmented 2D/3D images or tables with data on cells and/or nuclei sizes. In the article, we also represent some case studies showing the plugin application for solving biological tasks. Namely the plugin is applied in the following cases: defining parameters of jigsaw-puzzle pattern for maize leaf epidermal cells, analysis of the pavement cells morphological parameters for the mature wheat leaf grown under control and water deficit conditions, initiation of cell longitudinal rows, and detection of guard mother cells emergence at the initial stages of the stomatal morphogenesis in the growth zone of a wheat leaf.

Conclusion

The proposed plugin is efficient for high-throughput analysis of cellular architecture for cereal leaf epidermis. The workflow implies using inexpensive and rapid sample preparation and does not require the applying of transgenesis and reporter genetic structures expanding the range of species and varieties to study. Obtained characteristics of the cell structure and patterns further could act as a basis for the development and verification for spatial models of plant tissues formation mechanisms accounting for structural features of cereal leaves.

Availability

The implementation of this workflow is available as an ImageJ plugin distributed as a part of the Fiji project (FijiisjustImageJ: https://fiji.sc/). The plugin is freely available at https://imagej.net/LSM_Worker, https://github.com/JmanJ/LSM_Worker and http://pixie.bionet.nsc.ru/LSM_WORKER/.

Electronic supplementary material

The online version of this article (10.1186/s12918-019-0689-8) contains supplementary material, which is available to authorized users.

Quantitative characteristics of pubescence in wheat (Triticum aestivum L.) are associated with photosynthetic parameters under conditions of normal and limited water supply

Density and length of leaf pubescence are important factors of diversity in the response to water deficiency among wheat genotypes. Many studies evidence an important protective value of leaf hairiness in plants, especially under the conditions of drought, thermal loads and increased solar radiation. However, the physiological and adaptive roles of such traits in cereals, including cultivated plants, have not been sufficiently studied to date. The aim of this work was to study the association of morphological characteristics of leaves with parameters of gas exchange and chlorophyll fluorescence in wheat lines carrying a genetically different leaf hairiness. Isogenic and inter-varietal substitution wheat lines were used, carrying various combinations of dominant and recessive alleles of the known genes. A quantitative assessment of the pubescence was carried out in contrasting watering conditions to establish the physiological role of this trait in adaptation to drought. With the help of a portable system for studying the gas exchange and chlorophyll fluorescence, ten parameters of photosynthesis were studied, as well as morphological features of leaves and shoot biomass. It was found that gas exchange parameters are inversely proportional to the density and length of trichomes. In drought conditions, the trichome density increased and the length of trichomes decreased under the observed decrease in the level of gas exchange. A similar dependence was observed for the level of non-photochemical quenching of chlorophyll fluorescence. Under optimal conditions, the poorly haired cultivars exhibited a higher biomass than the densely haired. However, under water deficiency they significantly reduced the biomass and showed a low value of the tolerance index.

The evolution of gene regulatory networks controlling Arabidopsis thaliana L. trichome development

BACKGROUND

The variation in structure and function of gene regulatory networks (GRNs) participating in organisms development is a key for understanding species-specific evolutionary strategies. Even the tiniest modification of developmental GRN might result in a substantial change of a complex morphogenetic pattern. Great variety of trichomes and their accessibility makes them a useful model for studying the molecular processes of cell fate determination, cell cycle control and cellular morphogenesis. Nowadays, a large number of genes regulating the morphogenesis of A. thaliana trichomes are described. Here we aimed at a study the evolution of the GRN defining the trichome formation, and evaluation its importance in other developmental processes.

RESULTS

In study of the evolution of trichomes formation GRN we combined classical phylogenetic analysis with information on the GRN topology and composition in major plants taxa. This approach allowed us to estimate both times of evolutionary emergence of the GRN components which are mainly proteins, and the relative rate of their molecular evolution. Various simplifications of protein structure (based on the position of amino acid residues in protein globula, secondary structure type, and structural disorder) allowed us to demonstrate the evolutionary associations between changes in protein globules and speciations/duplications events. We discussed their potential involvement in protein-protein interactions and GRN function.

CONCLUSIONS

We hypothesize that the divergence and/or the specialization of the trichome-forming GRN is linked to the emergence of plant taxa. Information about the structural targets of the protein evolution in the GRN may predict switching points in gene networks functioning in course of evolution. We also propose a list of candidate genes responsible for the development of trichomes in a wide range of plant species.

Identification of nuclear genes controlling chlorophyll synthesis in barley by RNA-seq

BACKGROUND

Albinism in plants is characterized by lack of chlorophyll and results in photosynthesis impairment, abnormal plant development and premature death. These abnormalities are frequently encountered in interspecific crosses and tissue culture experiments. Analysis of albino mutant phenotypes with full or partial chlorophyll deficiency can shed light on genetic determinants and molecular mechanisms of albinism. Here we report analysis of RNA-seq transcription profiling of barley (Hordeum vulgare L.) near-isogenic lines, one of which is a carrier of mutant allele of the Alm gene for albino lemma and pericarp phenotype (line i:BwAlm).

RESULTS

1221 genome fragments have statistically significant changes in expression levels between lines i:BwAlm and Bowman, with 148 fragments having increased expression levels in line i:BwAlm, and 1073 genome fragments, including 42 plastid operons, having decreased levels of expression in line i:BwAlm. We detected functional dissimilarity between genes with higher and lower levels of expression in i:BwAlm line. Genes with lower level of expression in the i:BwAlm line are mostly associated with photosynthesis and chlorophyll synthesis, while genes with higher expression level are functionally associated with vesicle transport. Differentially expressed genes are shown to be involved in several metabolic pathways; the largest fraction of such genes was observed for the Calvin-Benson-Bassham cycle. Finally, de novo assembly of transcriptome contains several transcripts, not annotated in current H. vulgare genome version.

CONCLUSIONS

Our results provide the new information about genes which could be involved in formation of albino lemma and pericarp phenotype. They demonstrate the interplay between nuclear and chloroplast genomes in this physiological process.

[Methods of high-throughput plant phenotyping for large-scale breeding and genetic experiments]

Phenomics is a field of science at the junction of biology and informatics which solves the problems of rapid, accurate estimation of the plant phenotype; it was rapidly developed because of the need to analyze phenotypic characteristics in large scale genetic and breeding experiments in plants. It is based on using the methods of computer image analysis and integration of biological data. Owing to automation, new approaches make it possible to considerably accelerate the process of estimating the characteristics of a phenotype, to increase its accuracy, and to remove a subjectivism (inherent to humans). The main technologies of high-throughput plant phenotyping in both controlled and field conditions, their advantages and disadvantages, and also the prospects of their use for the efficient solution of problems of plant genetics and breeding are presented in the review.

[Genetic analysis of leaf pubescence in isogenic lines of conmmon wheat Novosibirsk 67]

In this study, genetic and monosomic analyses of the leaf pubescence of ANK 7A, ANK 7B, and ANK 7C wheat isogenic lines were carried out based on the Novosibirsk 67 wheat variety. According to visual analysis, the variety-recipient has a soft, uniform pubescence, and lines have trichomes on the surfaces of their leaves inherited from the two. Chinese varieties and one Soviet variety. Using the high throughput phenotyping method LHDetect2, which allows one to allocate the phenotypic classes of offspring in crosses based on the quantitative characteristics of leaf pubescence, it was found that chromosome 7B of the isogenic lines has a gene that determines the presence of long trichomes, and chromosome 7D of the Novosibirsk 67 variety has a gene that increases the density of pubescence. The obtained data allowed for the formulation of a hypothesis for the existence of a homoallelic series of genes that control leaf pubescence in the chromosomes of the seventh homeologous group of common wheat.

Extraction of quantitative characteristics describing wheat leaf pubescence with a novel image-processing technique

Leaf pubescence (hairiness) in wheat plays an important biological role in adaptation to the environment. However, this trait has always been methodologically difficult to phenotype. An important step forward has been taken with the use of computer technologies. Computer analysis of a photomicrograph of a transverse fold line of a leaf is proposed for quantitative evaluation of wheat leaf pubescence. The image-processing algorithm is implemented in the LHDetect2 software program accessible as a Web service at http://wheatdb.org/lhdetect2 . The results demonstrate that the proposed method is rapid, adequately assesses leaf pubescence density and the length distribution of trichomes and the data obtained using this method are significantly correlated with the density of trichomes on the leaf surface. Thus, the proposed method is efficient for high-throughput analysis of leaf pubescence morphology in cereal genetic collections and mapping populations.

Combined in silico/in vivo analysis of mechanisms providing for root apical meristem self-organization and maintenance

BACKGROUND AND AIMS

The root apical meristem (RAM) is the plant stem cell niche which provides for the formation and continuous development of the root. Auxin is the main regulator of RAM functioning, and auxin maxima coincide with the sites of RAM initiation and maintenance. Auxin gradients are formed due to local auxin biosynthesis and polar auxin transport. The PIN family of auxin transporters plays a critical role in polar auxin transport, and two mechanisms of auxin maximum formation in the RAM based on PIN-mediated auxin transport have been proposed to date: the reverse fountain and the reflected flow mechanisms.

METHODS

The two mechanisms are combined here in in silico studies of auxin distribution in intact roots and roots cut into two pieces in the proximal meristem region. In parallel, corresponding experiments were performed in vivo using DR5::GFP Arabidopsis plants.

KEY RESULTS

The reverse fountain and the reflected flow mechanism naturally cooperate for RAM patterning and maintenance in intact root. Regeneration of the RAM in decapitated roots is provided by the reflected flow mechanism. In the excised root tips local auxin biosynthesis either alone or in cooperation with the reverse fountain enables RAM maintenance.

CONCLUSIONS

The efficiency of a dual-mechanism model in guiding biological experiments on RAM regeneration and maintenance is demonstrated. The model also allows estimation of the concentrations of auxin and PINs in root cells during development and under various treatments. The dual-mechanism model proposed here can be a powerful tool for the study of several different aspects of auxin function in root.

[Morphological and genetic characteristics of leaf hairiness in wheat (Triticum aestivum L.) as analyzed by computer-aided phenotyping]

Computer-aided image processing was used to study the morphology of leaf hairiness in the wheat cultivars Saratovskaya 29 and Golubka, as well as the introgressed strain 102/00i of the cultivar Rodina carrying the hairiness control gene introgressed from Aegilops speltoides. Morphological differences in leaf hairiness were detected and described in detail. The genetic control ofhairiness was studied in two cultivars (Golubka and Saratovskaya 29) with similar hairiness patterns. Crossing these cultivars with the cultivar Rodina showed a monogenic inheritance in the cultivar Golubka and a digenic inheritance in the cultivar Saratovskaya 29, which has a denser hairiness. In the strain 102/00i and the cultivar Golubka, the number of trichomes was positively correlated with their mean length. The cultivar Golubka was used as an example to study the effect of environmental conditions on the formation of hairiness. Plants of these cultivars were found to form more but shorter trichomes.