A Workflow for the Performance of the Differential Ovarian Follicle Count Using Deep Neuronal Networks.

In order to automate the counting of ovarian follicles required in multigeneration reproductive studies performed in the rat according to Organization for Economic Co-operation and Development guidelines 443 and 416, the application of deep neural networks was tested. The manual evaluation of the differential ovarian follicle count is a tedious and time-consuming task that requires highly trained personnel. In this regard, deep learning outputs provide overlay pictures for a more detailed documentation, together with an increased reproducibility of the counts. To facilitate the planned good laboratory practice (GLP) validation a workflow was set up using MLFlow to make all steps from generating of scans, training of the neural network, uploading of study images to the neural network, generation and storage of the results in a compliant manner controllable and reproducible. PyTorch was used as main framework to build the Faster region-based convolutional neural network for the training. We compared the performances of different depths of ResNet models with specific regard to the sensitivity, specificity, accuracy of the models. In this paper, we describe all steps from data labeling, training of networks, and the performance metrics chosen to evaluate different network architectures. We also make recommendation on steps, which should be taken into consideration when GLP validation is aimed for.

Genome-wide association mapping in a diverse spring barley collection reveals the presence of QTL hotspots and candidate genes for root and shoot architecture traits at seedling stage.

BACKGROUND: Adaptation to drought-prone environments requires robust root architecture. Genotypes with a more vigorous root system have the potential to better adapt to soils with limited moisture content. However, root architecture is complex at both, phenotypic and genetic level. Customized mapping panels in combination with efficient screenings methods can resolve the underlying genetic factors of root traits. RESULTS: A mapping panel of 233 spring barley genotypes was evaluated for root and shoot architecture traits under non-stress and osmotic stress. A genome-wide association study elucidated 65 involved genomic regions. Among them were 34 root-specific loci, eleven hotspots with associations to up to eight traits and twelve stress-specific loci. A list of candidate genes was established based on educated guess. Selected genes were tested for associated polymorphisms. By this, 14 genes were identified as promising candidates, ten remained suggestive and 15 were rejected. The data support the important role of flowering time genes, including HvPpd-H1, HvCry2, HvCO4 and HvPRR73. Moreover, seven root-related genes, HERK2, HvARF04, HvEXPB1, PIN5, PIN7, PME5 and WOX5 are confirmed as promising candidates. For the QTL with the highest allelic effect for root thickness and plant biomass a homologue of the Arabidopsis Trx-m3 was revealed as the most promising candidate. CONCLUSIONS: This study provides a catalogue of hotspots for seedling growth, root and stress-specific genomic regions along with candidate genes for future potential incorporation in breeding attempts for enhanced yield potential, particularly in drought-prone environments. Root architecture is under polygenic control. The co-localization of well-known major genes for barley development and flowering time with QTL hotspots highlights their importance for seedling growth. Association analysis revealed the involvement of HvPpd-H1 in the development of the root system. The co-localization of root QTL with HERK2, HvARF04, HvEXPB1, PIN5, PIN7, PME5 and WOX5 represents a starting point to explore the roles of these genes in barley. Accordingly, the genes HvHOX2, HsfA2b, HvHAK2, and Dhn9, known to be involved in abiotic stress response, were located within stress-specific QTL regions and await future validation.

Genetic variation of growth dynamics in maize (Zea mays L.) revealed through automated non-invasive phenotyping.

Hitherto, most quantitative trait loci of maize growth and biomass yield have been identified for a single time point, usually the final harvest stage. Through this approach cumulative effects are detected, without considering genetic factors causing phase-specific differences in growth rates. To assess the genetics of growth dynamics, we employed automated non-invasive phenotyping to monitor the plant sizes of 252 diverse maize inbred lines at 11 different developmental time points; 50 k SNP array genotype data were used for genome-wide association mapping and genomic selection. The heritability of biomass was estimated to be over 71%, and the average prediction accuracy amounted to 0.39. Using the individual time point data, 12 main effect marker-trait associations (MTAs) and six pairs of epistatic interactions were detected that displayed different patterns of expression at various developmental time points. A subset of them also showed significant effects on relative growth rates in different intervals. The detected MTAs jointly explained up to 12% of the total phenotypic variation, decreasing with developmental progression. Using non-parametric functional mapping and multivariate mapping approaches, four additional marker loci affecting growth dynamics were detected. Our results demonstrate that plant biomass accumulation is a complex trait governed by many small effect loci, most of which act at certain restricted developmental phases. This highlights the need for investigation of stage-specific growth affecting genes to elucidate important processes operating at different developmental phases.

Dissecting the phenotypic components of crop plant growth and drought responses based on high-throughput image analysis.

Significantly improved crop varieties are urgently needed to feed the rapidly growing human population under changing climates. While genome sequence information and excellent genomic tools are in place for major crop species, the systematic quantification of phenotypic traits or components thereof in a high-throughput fashion remains an enormous challenge. In order to help bridge the genotype to phenotype gap, we developed a comprehensive framework for high-throughput phenotype data analysis in plants, which enables the extraction of an extensive list of phenotypic traits from nondestructive plant imaging over time. As a proof of concept, we investigated the phenotypic components of the drought responses of 18 different barley (Hordeum vulgare) cultivars during vegetative growth. We analyzed dynamic properties of trait expression over growth time based on 54 representative phenotypic features. The data are highly valuable to understand plant development and to further quantify growth and crop performance features. We tested various growth models to predict plant biomass accumulation and identified several relevant parameters that support biological interpretation of plant growth and stress tolerance. These image-based traits and model-derived parameters are promising for subsequent genetic mapping to uncover the genetic basis of complex agronomic traits. Taken together, we anticipate that the analytical framework and analysis results presented here will be useful to advance our views of phenotypic trait components underlying plant development and their responses to environmental cues.

Non-canonical structure, function and phylogeny of the Bsister MADS-box gene OsMADS30 of rice (Oryza sativa).

Bsister MADS-box genes play key roles in female reproductive organ and seed development throughout seed plants. This view is supported by their high conservation in terms of sequence, expression and function. In grasses, there are three subclades of Bsister genes: the OsMADS29-, the OsMADS30- and the OsMADS31-like genes. Here, we report on the evolution of the OsMADS30-like genes. Our analyses indicate that these genes evolved under relaxed purifying selection and are rather weakly expressed. OsMADS30, the representative of the OsMADS30-like genes from rice (Oryza sativa), shows strong sequence deviations in its 3' region when compared to orthologues from other grass species. We show that this is due to a 2.4-kbp insertion, possibly of a hitherto unknown helitron, which confers a heterologous C-terminal domain to OsMADS30. This putative helitron is not present in the OsMADS30 orthologues from closely related wild rice species, pointing to a relatively recent insertion event. Unlike other Bsister mutants O. sativa plants carrying a T-DNA insertion in the OsMADS30 gene do not show aberrant seed phenotypes, indicating that OsMADS30 likely does not have a canonical 'Bsister function'. However, imaging-based phenotyping of the T-DNA carrying plants revealed alterations in shoot size and architecture. We hypothesize that sequence deviations that accumulated during a period of relaxed selection in the gene lineage that led to OsMADS30 and the alteration of the C-terminal domain might have been a precondition for a potential neo-functionalization of OsMADS30 in O. sativa.

Dissecting the chromatin interactome of microRNA genes.

Our knowledge of the role of higher-order chromatin structures in transcription of microRNA genes (MIRs) is evolving rapidly. Here we investigate the effect of 3D architecture of chromatin on the transcriptional regulation of MIRs. We demonstrate that MIRs have transcriptional features that are similar to protein-coding genes. RNA polymerase II-associated ChIA-PET data reveal that many groups of MIRs and protein-coding genes are organized into functionally compartmentalized chromatin communities and undergo coordinated expression when their genomic loci are spatially colocated. We observe that MIRs display widespread communication in those transcriptionally active communities. Moreover, miRNA-target interactions are significantly enriched among communities with functional homogeneity while depleted from the same community from which they originated, suggesting MIRs coordinating function-related pathways at posttranscriptional level. Further investigation demonstrates the existence of spatial MIR-MIR chromatin interacting networks. We show that groups of spatially coordinated MIRs are frequently from the same family and involved in the same disease category. The spatial interaction network possesses both common and cell-specific subnetwork modules that result from the spatial organization of chromatin within different cell types. Together, our study unveils an entirely unexplored layer of MIR regulation throughout the human genome that links the spatial coordination of MIRs to their co-expression and function.

Integrated Analysis Platform: An Open-Source Information System for High-Throughput Plant Phenotyping.

High-throughput phenotyping is emerging as an important technology to dissect phenotypic components in plants. Efficient image processing and feature extraction are prerequisites to quantify plant growth and performance based on phenotypic traits. Issues include data management, image analysis, and result visualization of large-scale phenotypic data sets. Here, we present Integrated Analysis Platform (IAP), an open-source framework for high-throughput plant phenotyping. IAP provides user-friendly interfaces, and its core functions are highly adaptable. Our system supports image data transfer from different acquisition environments and large-scale image analysis for different plant species based on real-time imaging data obtained from different spectra. Due to the huge amount of data to manage, we utilized a common data structure for efficient storage and organization of data for both input data and result data. We implemented a block-based method for automated image processing to extract a representative list of plant phenotypic traits. We also provide tools for build-in data plotting and result export. For validation of IAP, we performed an example experiment that contains 33 maize (Zea mays 'Fernandez') plants, which were grown for 9 weeks in an automated greenhouse with nondestructive imaging. Subsequently, the image data were subjected to automated analysis with the maize pipeline implemented in our system. We found that the computed digital volume and number of leaves correlate with our manually measured data in high accuracy up to 0.98 and 0.95, respectively. In summary, IAP provides a multiple set of functionalities for import/export, management, and automated analysis of high-throughput plant phenotyping data, and its analysis results are highly reliable.

Dissecting spatiotemporal biomass accumulation in barley under different water regimes using high-throughput image analysis.

Phenotyping large numbers of genotypes still represents the rate-limiting step in many plant genetic experiments and in breeding. To address this issue, novel automated phenotyping technologies have been developed. We investigated for a core set of barley cultivars if high-throughput image analysis can help to dissect vegetative biomass accumulation in response to two different watering regimes under semi-controlled greenhouse conditions. We found that experiments, treatments, genotypes and genotype by environment interaction (G × E) can be characterized at any time point by certain digital traits. Biomass accumulation under control and stress conditions was highly heritable. Growth model-derived maximum vegetative biomass (K max), inflection point (I) and regrowth rate (k) were identified as promising candidate traits for genome-wide association studies. Drought stress symptoms can be visualized, dissected and modelled. Especially the highly heritable regrowth rate, which had the biggest influence on biomass accumulation in stress treatment, seems promising for future studies to improve drought tolerance in different crop species. A proof of concept study revealed potential correlations between digital traits obtained from pot experiments under greenhouse conditions and agronomic traits from field experiments. Overall, non-invasive, imaging-based phenotyping platforms under greenhouse conditions offer excellent possibilities for trait discovery, trait development and industrial applications.

Towards recommendations for metadata and data handling in plant phenotyping.

Recent methodological developments in plant phenotyping, as well as the growing importance of its applications in plant science and breeding, are resulting in a fast accumulation of multidimensional data. There is great potential for expediting both discovery and application if these data are made publicly available for analysis. However, collection and storage of phenotypic observations is not yet sufficiently governed by standards that would ensure interoperability among data providers and precisely link specific phenotypes and associated genomic sequence information. This lack of standards is mainly a result of a large variability of phenotyping protocols, the multitude of phenotypic traits that are measured, and the dependence of these traits on the environment. This paper discusses the current situation of standardization in the area of phenomics, points out the problems and shortages, and presents the areas that would benefit from improvement in this field. In addition, the foundations of the work that could revise the situation are proposed, and practical solutions developed by the authors are introduced.

The differential transcription network between embryo and endosperm in the early developing maize seed.

Transcriptome analysis of early-developing maize (Zea mays) seed was conducted using Illumina sequencing. We mapped 11,074,508 and 11,495,788 paired-end reads from endosperm and embryo, respectively, at 9 d after pollination to define gene structure and alternative splicing events as well as transcriptional regulators of gene expression to quantify transcript abundance in both embryo and endosperm. We identified a large number of novel transcribed regions that did not fall within maize annotated regions, and many of the novel transcribed regions were tissue-specifically expressed. We found that 50.7% (8,556 of 16,878) of multiexonic genes were alternatively spliced, and some transcript isoforms were specifically expressed either in endosperm or in embryo. In addition, a total of 46 trans-splicing events, with nine intrachromosomal events and 37 interchromosomal events, were found in our data set. Many metabolic activities were specifically assigned to endosperm and embryo, such as starch biosynthesis in endosperm and lipid biosynthesis in embryo. Finally, a number of transcription factors and imprinting genes were found to be specifically expressed in embryo or endosperm. This data set will aid in understanding how embryo/endosperm development in maize is differentially regulated.

Apomictic and sexual ovules of Boechera display heterochronic global gene expression patterns.

We have compared the transcriptomic profiles of microdissected live ovules at four developmental stages between a diploid sexual and diploid apomictic Boechera. We sequenced >2 million SuperSAGE tags and identified (1) heterochronic tags (n = 595) that demonstrated significantly different patterns of expression between sexual and apomictic ovules across all developmental stages, (2) stage-specific tags (n = 577) that were found in a single developmental stage and differentially expressed between the sexual and apomictic ovules, and (3) sex-specific (n = 237) and apomixis-specific (n = 1106) tags that were found in all four developmental stages but in only one reproductive mode. Most heterochronic and stage-specific tags were significantly downregulated during early apomictic ovule development, and 110 were associated with reproduction. By contrast, most late stage-specific tags were upregulated in the apomictic ovules, likely the result of increased gene copy number in apomictic (hexaploid) versus sexual (triploid) endosperm or of parthenogenesis. Finally, we show that apomixis-specific gene expression is characterized by a significant overrepresentation of transcription factor activity. We hypothesize that apomeiosis is associated with global downregulation at the megaspore mother cell stage. As the diploid apomict analyzed here is an ancient hybrid, these data are consistent with the postulated link between hybridization and asexuality and provide a hypothesis for multiple evolutionary origins of apomixis in the genus Boechera.

Creating views on integrated multidomain data.

MOTIVATION: Modern data acquisition methods in biology allow the procurement of different types of data in increasing quantity, facilitating a comprehensive view of biological systems. As data are usually gathered and interpreted by separate domain scientists, it is hard to grasp multidomain properties and structures. Consequently, there is a need for the integration of biological data from different sources and of different types in one application, providing various visualization approaches. RESULTS: In this article, methods for the integration and visualization of multimodal biological data are presented. This is achieved based on two graphs representing the meta-relations between biological data and the measurement combinations, respectively. Both graphs are linked and serve as different views of the integrated data with navigation and exploration possibilities. Data can be combined and visualized multifariously, resulting in views of the integrated biological data. AVAILABILITY: http://vanted.ipk-gatersleben.de/hive/. CONTACT: rohn@ipk-gatersleben.de.

MyBioNet: interactively visualize, edit and merge biological networks on the web.

SUMMARY: MyBioNet is a web-based application for biological network analysis, which provides user-friendly web interfaces to visualize, edit and merge biological networks. In addition, MyBioNet integrated KEGG metabolic network data from 1366 organisms and allows users to search and navigate interesting networks. AVAILABILITY AND IMPLEMENTATION: All KEGG metabolic network data are organized and stored in the MySQL database. MyBioNet is implemented in Flex/Actionscript and PHP languages and deployed on an Apache web server. MyBioNet is accessible through all the Flash-embedded browsers at http://bis.zju.edu.cn/mybionet/. CONTACT: mchen@zju.edu.cn.

Precision phenotyping across the life cycle to validate and decipher drought-adaptive QTLs of wild emmer wheat (Triticum turgidum ssp. dicoccoides) introduced into elite wheat varieties.

Drought events or the combination of drought and heat conditions are expected to become more frequent due to global warming, and wheat yields may fall below their long-term average. One way to increase climate-resilience of modern high-yielding varieties is by their genetic improvement with beneficial alleles from crop wild relatives. In the present study, the effect of two beneficial QTLs introgressed from wild emmer wheat and incorporated in the three wheat varieties BarNir, Zahir and Uzan was studied under well-watered conditions and under drought stress using non-destructive High-throughput Phenotyping (HTP) throughout the life cycle in a single pot-experiment. Plants were daily imaged with RGB top and side view cameras and watered automatically. Further, at two time points, the quantum yield of photosystem II was measured with a top view FluorCam. The QTL carrying near isogenic lines (NILs) were compared with their corresponding parents by t-test for all non-invasively obtained traits and for the manually determined agronomic and yield parameters. Data quality of phenotypic traits (repeatability) in the controlled HTP experiment was above 85% throughout the life cycle and at maturity. Drought stress had a strong effect on growth in all wheat genotypes causing biomass reduction from 2% up to 70% at early and late points in the drought period, respectively. At maturity, the drought caused 47-55% decreases in yield-related traits grain weight, straw weight and total biomass and reduced TKW by 10%, while water use efficiency (WUE) increased under drought by 29%. The yield-enhancing effect of the introgressed QTLs under drought conditions that were previously demonstrated under field/screenhouse conditions in Israel, could be mostly confirmed in a greenhouse pot experiment using HTP. Daily precision phenotyping enabled to decipher the mode of action of the QTLs in the different genetic backgrounds throughout the entire wheat life cycle. Daily phenotyping allowed a precise determination of the timing and size of the QTLs effect (s) and further yielded information about which image-derived traits are informative at which developmental stage of wheat during the entire life cycle. Maximum height and estimated biovolume were reached about a week after heading, so experiments that only aim at exploring these traits would not need a longer observation period. To obtain information on different onset and progress of senescence, the CVa curves represented best the ongoing senescence of plants. The QTL on 7A in the BarNir background was found to improve yield under drought by increased biomass growth, a higher photosynthetic performance, a higher WUE and a "stay green effect."

Analysis of Few-Shot Techniques for Fungal Plant Disease Classification and Evaluation of Clustering Capabilities Over Real Datasets.

Plant fungal diseases are one of the most important causes of crop yield losses. Therefore, plant disease identification algorithms have been seen as a useful tool to detect them at early stages to mitigate their effects. Although deep-learning based algorithms can achieve high detection accuracies, they require large and manually annotated image datasets that is not always accessible, specially for rare and new diseases. This study focuses on the development of a plant disease detection algorithm and strategy requiring few plant images (Few-shot learning algorithm). We extend previous work by using a novel challenging dataset containing more than 100,000 images. This dataset includes images of leaves, panicles and stems of five different crops (barley, corn, rape seed, rice, and wheat) for a total of 17 different diseases, where each disease is shown at different disease stages. In this study, we propose a deep metric learning based method to extract latent space representations from plant diseases with just few images by means of a Siamese network and triplet loss function. This enhances previous methods that require a support dataset containing a high number of annotated images to perform metric learning and few-shot classification. The proposed method was compared over a traditional network that was trained with the cross-entropy loss function. Exhaustive experiments have been performed for validating and measuring the benefits of metric learning techniques over classical methods. Results show that the features extracted by the metric learning based approach present better discriminative and clustering properties. Davis-Bouldin index and Silhouette score values have shown that triplet loss network improves the clustering properties with respect to the categorical-cross entropy loss. Overall, triplet loss approach improves the DB index value by 22.7% and Silhouette score value by 166.7% compared to the categorical cross-entropy loss model. Moreover, the F-score parameter obtained from the Siamese network with the triplet loss performs better than classical approaches when there are few images for training, obtaining a 6% improvement in the F-score mean value. Siamese networks with triplet loss have improved the ability to learn different plant diseases using few images of each class. These networks based on metric learning techniques improve clustering and classification results over traditional categorical cross-entropy loss networks for plant disease identification.

Editing, validating and translating of SBGN maps.

MOTIVATION: The recently proposed Systems Biology Graphical Notation (SBGN) provides a standard for the visual representation of biochemical and cellular processes. It aims to support more efficient and accurate communication of biological knowledge between different research communities in the life sciences. However, to increase the use of SBGN, tools for editing, validating and translating SBGN maps are desirable. RESULTS: We present SBGN-ED, a tool which allows the creation of all three types of SBGN maps from scratch or the editing of existing maps, the validation of these maps for syntactical and semantical correctness, the translation of networks from the KEGG and MetaCrop databases into SBGN and the export of SBGN maps into several file and image formats. AVAILABILITY: SBGN-ED is freely available from http://vanted.ipk-gatersleben.de/addons/sbgn-ed. The web site contains also tutorials and example files.

Lessons from GNE-deficient embryonic stem cells: sialic acid biosynthesis is involved in proliferation and gene expression.

Sialic acids are widely expressed as terminal carbohydrates on glycoconjugates of eukaryotic cells. They are involved in a variety of cellular functions, such as cell adhesion or signal recognition. The key enzyme of sialic acid biosynthesis is the bifunctional UDP-N-acetylglucosamine-2-epimerase/N-acetylmannosamine kinase (GNE), which catalyzes the first two steps of sialic acid biosynthesis in the cytosol. Inactivation of GNE causes early embryonic lethality. In this study, we analyzed wild-type and GNE-deficient embryonic stem cells from mice. We found for the first time that proliferation is directly correlated with GNE-expression and the cellular sialic acid concentration. Furthermore, we identified growth-related genes that are differentially expressed in GNE-deficient embryonic stem cells compared to wild-type embryonic stem cells.

FBA-SimVis: interactive visualization of constraint-based metabolic models.

UNLABELLED: FBA-SimVis is a VANTED plug-in for the constraint-based analysis of metabolic models with special focus on the visual exploration of metabolic flux data resulting from model analysis. The program provides a user-friendly environment for model reconstruction, constraint-based model analysis, and interactive visualization of the simulation results. With the ability to quantitatively analyse metabolic fluxes in an interactive and visual manner, FBA-SimVis supports a comprehensive understanding of constraint-based metabolic flux models in both overview and detail. AVAILABILITY: Software with manual and tutorials are freely available at http://fbasimvis.ipk-gatersleben.de/

Mental-Map Preserving Visualisation of Partitioned Networks in Vanted.

Biological networks can be large and complex, often consisting of different sub-networks or parts. Separation of networks into parts, network partitioning and layouts of overview and sub-graphs are of importance for understandable visualisations of those networks. This article presents NetPartVis to visualise non-overlapping clusters or partitions of graphs in the Vanted framework based on a method for laying out overview graph and several sub-graphs (partitions) in a coordinated, mental-map preserving way.