Development of an automated 3D high content cell screening platform for organoid phenotyping.

The use of organoid models in biomedical research has grown substantially since their inception. As they gain popularity among scientists seeking more complex and biologically relevant systems, there is a direct need to expand and clarify potential uses of such systems in diverse experimental contexts. Herein we outline a high-content screening (HCS) platform that allows researchers to screen drugs or other compounds against three-dimensional (3D) cell culture systems in a multi-well format (384-well). Furthermore, we compare the quality of robotic liquid handling with manual pipetting and characterize and contrast the phenotypic effects detected by confocal imaging and biochemical assays in response to drug treatment. We show that robotic liquid handling is more consistent and amendable to high throughput experimental designs when compared to manual pipetting due to improved precision and automated randomization capabilities. We also show that image-based techniques are more sensitive to detecting phenotypic changes within organoid cultures than traditional biochemical assays that evaluate cell viability, supporting their integration into organoid screening workflows. Finally, we highlight the enhanced capabilities of confocal imaging in this organoid screening platform as they relate to discerning organoid drug responses in single-well co-cultures of organoids derived from primary human biopsies and patient-derived xenograft (PDX) models. Altogether, this platform enables automated, imaging-based HCS of 3D cellular models in a non-destructive manner, opening the path to complementary analysis through integrated downstream methods.

Stratification of Lung Cancer Risk with Thoracic Imaging Phenotypes.

In lung cancer screening, estimation of future lung cancer risk is usually guided by demographics and smoking status. The role of constitutional profiles of human body, a.k.a. body habitus, is increasingly understood to be important, but has not been integrated into risk models. Chest low dose computed tomography (LDCT) is the standard imaging study in lung cancer screening, with the capability to discriminate differences in body composition and organ arrangement in the thorax. We hypothesize that the primary phenotypes identified using lung screening chest LDCT can form a representation of body habitus and add predictive power for lung cancer risk stratification. In this pilot study, we evaluated the feasibility of body habitus image-based phenotyping on a large lung screening LDCT dataset. A thoracic imaging manifold was estimated based on an intensity-based pairwise (dis)similarity metric for pairs of spatial normalized chest LDCT images. We applied the hierarchical clustering method on this manifold to identify the primary phenotypes. Body habitus features of each identified phenotype were evaluated and associated with future lung cancer risk using time-to-event analysis. We evaluated the method on the baseline LDCT scans of 1,200 male subjects sampled from National Lung Screening Trial. Five primary phenotypes were identified, which were associated with highly distinguishable clinical and body habitus features. Time-to-event analysis against future lung cancer incidences showed two of the five identified phenotypes were associated with elevated future lung cancer risks (HR=1.61, 95% CI = [1.08, 2.38], p=0.019; HR=1.67, 95% CI = [0.98, 2.86], p=0.057). These results indicated that it is feasible to capture the body habitus by image-base phenotyping using lung screening LDCT and the learned body habitus representation can potentially add value for future lung cancer risk stratification.

A Lycopene ε-Cyclase TILLING Allele Enhances Lycopene and Carotenoid Content in Fruit and Improves Drought Stress Tolerance in Tomato Plants.

In the scenario of climate change, the availability of genetic resources for tomato cultivation that combine improved nutritional properties and more tolerance to water deficiency is highly desirable. Within this context, the molecular screenings of the Red Setter cultivar-based TILLING platform led to the isolation of a novel lycopene ε-cyclase gene (SlLCY-E) variant (G/3378/T) that produces modifications in the carotenoid content of tomato leaves and fruits. In leaf tissue, the novel G/3378/T SlLCY-E allele enhances ß,ß-xanthophyll content at the expense of lutein, which decreases, while in ripe tomato fruit the TILLING mutation induces a significant increase in lycopene and total carotenoid content. Under drought stress conditions, the G/3378/T SlLCY-E plants produce more abscisic acid (ABA) and still conserve their leaf carotenoid profile (reduction of lutein and increase in ß,ß-xanthophyll content). Furthermore, under said conditions, the mutant plants grow much better and are more tolerant to drought stress, as revealed by digital-based image analysis and in vivo monitoring of the OECT (Organic Electrochemical Transistor) sensor. Altogether, our data indicate that the novel TILLING SlLCY-E allelic variant is a valuable genetic resource that can be used for developing new tomato varieties, improved in drought stress tolerance and enriched in fruit lycopene and carotenoid content.

Image-based time series analysis to establish differential disease progression for two Fusarium head blight pathogens in oat spikelets with variable resistance.

Oat-based value-added products have increased their value as healthy foodstuff. Fusarium head blight (FHB) infections and the mycotoxins accumulated to the oat seeds, however, pose a challenge to oat production. The FHB infections are predicted to become more prevalent in the future changing climates and under more limited use of fungicides. Both these factors increase the pressure for breeding new resistant cultivars. Until now, however, genetic links in oats against FHB infection have been difficult to identify. Therefore, there is a great need for more effective breeding efforts, including improved phenotyping methods allowing time series analysis and the identification of molecular markers during disease progression. To these ends, dissected spikelets of several oat genotypes with different resistance profiles were studied by image-based methods during disease progression by Fusarium culmorum or F. langsethiae species. The chlorophyll fluorescence of each pixel in the spikelets was recorded after inoculation by the two Fusarium spp., and the progression of the infections was analyzed by calculating the mean maximum quantum yield of PSII (Fv/Fm) values for each spikelet. The recorded values were (i) the change in the photosynthetically active area of the spikelet as percentage of its initial size, and (ii) the mean of Fv/Fm values of all fluorescent pixels per spikelet post inoculation, both indicative of the progression of the FHB disease. The disease progression was successfully monitored, and different stages of the infection could be defined along the time series. The data also confirmed the differential rate of disease progression by the two FHB causal agents. In addition, oat varieties with variable responses to the infections were indicated.

Deep learning-based phenotyping for genome wide association studies of sudden death syndrome in soybean.

Using a reliable and accurate method to phenotype disease incidence and severity is essential to unravel the complex genetic architecture of disease resistance in plants, and to develop disease resistant cultivars. Genome-wide association studies (GWAS) involve phenotyping large numbers of accessions, and have been used for a myriad of traits. In field studies, genetic accessions are phenotyped across multiple environments and replications, which takes a significant amount of labor and resources. Deep Learning (DL) techniques can be effective for analyzing image-based tasks; thus DL methods are becoming more routine for phenotyping traits to save time and effort. This research aims to conduct GWAS on sudden death syndrome (SDS) of soybean [Glycine max L. (Merr.)] using disease severity from both visual field ratings and DL-based (using images) severity ratings collected from 473 accessions. Images were processed through a DL framework that identified soybean leaflets with SDS symptoms, and then quantified the disease severity on those leaflets into a few classes with mean Average Precision of 0.34 on unseen test data. Both visual field ratings and image-based ratings identified significant single nucleotide polymorphism (SNP) markers associated with disease resistance. These significant SNP markers are either in the proximity of previously reported candidate genes for SDS or near potentially novel candidate genes. Four previously reported SDS QTL were identified that contained a significant SNPs, from this study, from both a visual field rating and an image-based rating. The results of this study provide an exciting avenue of using DL to capture complex phenotypic traits from images to get comparable or more insightful results compared to subjective visual field phenotyping of traits for disease symptoms.

Genetic analysis reveals three novel QTLs underpinning a butterfly egg-induced hypersensitive response-like cell death in Brassica rapa.

BACKGROUND: Cabbage white butterflies (Pieris spp.) can be severe pests of Brassica crops such as Chinese cabbage, Pak choi (Brassica rapa) or cabbages (B. oleracea). Eggs of Pieris spp. can induce a hypersensitive response-like (HR-like) cell death which reduces egg survival in the wild black mustard (B. nigra). Unravelling the genetic basis of this egg-killing trait in Brassica crops could improve crop resistance to herbivory, reducing major crop losses and pesticides use. Here we investigated the genetic architecture of a HR-like cell death induced by P. brassicae eggs in B. rapa. RESULTS: A germplasm screening of 56 B. rapa accessions, representing the genetic and geographical diversity of a B. rapa core collection, showed phenotypic variation for cell death. An image-based phenotyping protocol was developed to accurately measure size of HR-like cell death and was then used to identify two accessions that consistently showed weak (R-o-18) or strong cell death response (L58). Screening of 160 RILs derived from these two accessions resulted in three novel QTLs for Pieris brassicae-induced cell death on chromosomes A02 (Pbc1), A03 (Pbc2), and A06 (Pbc3). The three QTLs Pbc1-3 contain cell surface receptors, intracellular receptors and other genes involved in plant immunity processes, such as ROS accumulation and cell death formation. Synteny analysis with A. thaliana suggested that Pbc1 and Pbc2 are novel QTLs associated with this trait, while Pbc3 also contains an ortholog of LecRK-I.1, a gene of A. thaliana previously associated with cell death induced by a P. brassicae egg extract. CONCLUSIONS: This study provides the first genomic regions associated with the Pieris egg-induced HR-like cell death in a Brassica crop species. It is a step closer towards unravelling the genetic basis of an egg-killing crop resistance trait, paving the way for breeders to further fine-map and validate candidate genes.

Climate Change Impact on Wheat Performance-Effects on Vigour, Plant Traits and Yield from Early and Late Drought Stress in Diverse Lines.

Global climate change is threatening wheat productivity; improved yield under drought conditions is urgent. Here, diverse spring-wheat lines (modern, old and wheat-rye introgressions) were examined in an image-based early-vigour assay and a controlled-conditions (Biotron) trial that evaluated 13 traits until maturity. Early root vigour was significantly higher in the old Swedish lines (root length 8.50 cm) and introgressed lines with 1R (11.78 cm) and 1RS (9.91 cm) than in the modern (4.20 cm) and 2R (4.67 cm) lines. No significant correlation was noted between early root and shoot vigour. A higher yield was obtained under early drought stress in the 3R genotypes than in the other genotype groups, while no clear patterns were noted under late drought. Evaluating the top 10% of genotypes in terms of the stress-tolerance index for yield showed that root biomass, grains and spikes per plant were accountable for tolerance to early drought, while 1000-grain weight and flag-leaf area were accountable for tolerance to late drought. Early root vigour was determined as an important focus trait of wheat breeding for tolerance to climate-change-induced drought. The responsible genes for the trait should be searched for in these diverse lines. Additional drought-tolerance traits determined here need further elaboration to identify the responsible genes.

Genome-wide association studies dissect the genetic architecture of seed shape and size in common bean.

Seed weight and size are important yield components. Thus, selecting for large seeds has been a key objective in crop domestication and breeding. In common bean, seed shape is also important since it influences industrial processing and plays a vital role in determining the choices of consumers and farmers. In this study, we performed genome-wide association studies on a core collection of common bean accessions to dissect the genetic architecture and identify genomic regions associated with seed morphological traits related to weight, size, and shape. Phenotypic data were collected by high-throughput image-based approaches, and utilized to test associations with 10,362 single-nucleotide polymorphism markers using multilocus mixed models. We searched within genome-associated regions for candidate genes putatively involved in seed phenotypic variation. The collection exhibited high variability for the entire set of seed traits, and the Andean gene pool was found to produce larger, heavier seeds than the Mesoamerican gene pool. Strong pairwise correlations were verified for most seed traits. Genome-wide association studies identified marker-trait associations accounting for a considerable amount of phenotypic variation in length, width, projected area, perimeter, and circularity in 4 distinct genomic regions. Promising candidate genes were identified, e.g. those encoding an AT-hook motif nuclear-localized protein 8, type 2C protein phosphatases, and a protein Mei2-like 4 isoform, known to be associated with seed size and weight regulation. Moreover, the genes that were pinpointed are also good candidates for functional analysis to validate their influence on seed shape and size in common bean and other related crops.

Radiomics: a primer on high-throughput image phenotyping.

Radiomics is a high-throughput approach to image phenotyping. It uses computer algorithms to extract and analyze a large number of quantitative features from radiological images. These radiomic features collectively describe unique patterns that can serve as digital fingerprints of disease. They may also capture imaging characteristics that are difficult or impossible to characterize by the human eye. The rapid development of this field is motivated by systems biology, facilitated by data analytics, and powered by artificial intelligence. Here, as part of Abdominal Radiology's special issue on Quantitative Imaging, we provide an introduction to the field of radiomics. The technique is formally introduced as an advanced application of data analytics, with illustrating examples in abdominal radiology. Artificial intelligence is then presented as the main driving force of radiomics, and common techniques are defined and briefly compared. The complete step-by-step process of radiomic phenotyping is then broken down into five key phases. Potential pitfalls of each phase are highlighted, and recommendations are provided to reduce sources of variation, non-reproducibility, and error associated with radiomics.

Review: Application of Artificial Intelligence in Phenomics.

Plant phenomics has been rapidly advancing over the past few years. This advancement is attributed to the increased innovation and availability of new technologies which can enable the high-throughput phenotyping of complex plant traits. The application of artificial intelligence in various domains of science has also grown exponentially in recent years. Notably, the computer vision, machine learning, and deep learning aspects of artificial intelligence have been successfully integrated into non-invasive imaging techniques. This integration is gradually improving the efficiency of data collection and analysis through the application of machine and deep learning for robust image analysis. In addition, artificial intelligence has fostered the development of software and tools applied in field phenotyping for data collection and management. These include open-source devices and tools which are enabling community driven research and data-sharing, thereby availing the large amounts of data required for the accurate study of phenotypes. This paper reviews more than one hundred current state-of-the-art papers concerning AI-applied plant phenotyping published between 2010 and 2020. It provides an overview of current phenotyping technologies and the ongoing integration of artificial intelligence into plant phenotyping. Lastly, the limitations of the current approaches/methods and future directions are discussed.

Genetic improvement of wheat early vigor promote weed-competitiveness under Mediterranean climate.

Chemical weed-control is the most effective practice for wheat, however, rapid evolution of herbicide-resistant weeds threat food-security and calls for integration of non-chemical practices. We hypothesis that integration of alternative GA-responsive dwarfing genes into elite wheat cultivars can promote early vigor and weed-competitiveness under Mediterranean climate. We develop near-isogenic lines of bread wheat cultivars with GAR dwarfing genes and evaluate them for early vigor and weed-competitiveness under various environmental and management conditions to identify promising NIL for weed-competitiveness and grain yield. While all seven NILs responded to external gibberellic acid application, they exhibited differences in early vigor. Greenhouse and field evaluations highlighted NIL OC1 (Rht8andRht12) as a promising line, with significant advantage in canopy early vigor over its parental. To facilitate accurate and continuous early vigor data collection, we applied non-destructive image-based phenotyping approaches which offers non-expensive and end-user friendly solution for selection. NIL OC1 was tested under different weed density level, infestation waves, and temperatures and highlight the complex genotypic × environmental × management interactions. Our findings demonstrate the potential of genetic modification of dwarfing genes as promising approach to improve weed-competitiveness, and serve as basis for future breeding efforts to support sustainable wheat production under semi-arid Mediterranean climate.

Drought Tolerant Near Isogenic Lines of Pusa 44 Pyramided With qDTY2.1 and qDTY3.1, Show Accelerated Recovery Response in a High Throughput Phenomics Based Phenotyping.

Reproductive stage drought stress (RSDS) is a major challenge in rice production worldwide. Cultivar development with drought tolerance has been slow due to the lack of precise high throughput phenotyping tools to quantify drought stress-induced effects. Most of the available techniques are based on destructive sampling and do not assess the progress of the plant's response to drought. In this study, we have used state-of-the-art image-based phenotyping in a phenomics platform that offers a controlled environment, non-invasive phenotyping, high accuracy, speed, and continuity. In rice, several quantitative trait loci (QTLs) which govern grain yield under drought determine RSDS tolerance. Among these, qDTY2.1 and qDTY3.1 were used for marker-assisted breeding. A set of 35 near-isogenic lines (NILs), introgressed with these QTLs in the popular variety, Pusa 44 were used to assess the efficiency of image-based phenotyping for RSDS tolerance. NILs offered the most reliable contrast since they differed from Pusa 44 only for the QTLs. Four traits, namely, the projected shoot area (PSA), water use (WU), transpiration rate (TR), and red-green-blue (RGB) and near-infrared (NIR) values were used. Differential temporal responses could be seen under drought, but not under unstressed conditions. NILs showed significant level of RSDS tolerance as compared to Pusa 44. Among the traits, PSA showed strong association with yield (80%) as well as with two drought tolerances indices, stress susceptibility index (SSI) and tolerance index (TOL), establishing its ability in identifying the best drought tolerant NILs. The results revealed that the introgression of QTLs helped minimize the mean WU per unit of biomass per day, suggesting the potential role of these QTLs in improving WU-efficiency (WUE). We identified 11 NILs based on phenomics traits as well as performance under imposed drought in the field. The study emphasizes the use of phenomics traits as selection criteria for RSDS tolerance at an early stage, and is the first report of using phenomics parameters in RSDS selection in rice.

Challenges and Opportunities in Machine-Augmented Plant Stress Phenotyping.

Plant stress phenotyping is essential to select stress-resistant varieties and develop better stress-management strategies. Standardization of visual assessments and deployment of imaging techniques have improved the accuracy and reliability of stress assessment in comparison with unaided visual measurement. The growing capabilities of machine learning (ML) methods in conjunction with image-based phenotyping can extract new insights from curated, annotated, and high-dimensional datasets across varied crops and stresses. We propose an overarching strategy for utilizing ML techniques that methodically enables the application of plant stress phenotyping at multiple scales across different types of stresses, program goals, and environments.

Time-resolved mapping of genetic interactions to model rewiring of signaling pathways.

Context-dependent changes in genetic interactions are an important feature of cellular pathways and their varying responses under different environmental conditions. However, methodological frameworks to investigate the plasticity of genetic interaction networks over time or in response to external stresses are largely lacking. To analyze the plasticity of genetic interactions, we performed a combinatorial RNAi screen in Drosophila cells at multiple time points and after pharmacological inhibition of Ras signaling activity. Using an image-based morphology assay to capture a broad range of phenotypes, we assessed the effect of 12768 pairwise RNAi perturbations in six different conditions. We found that genetic interactions form in different trajectories and developed an algorithm, termed MODIFI, to analyze how genetic interactions rewire over time. Using this framework, we identified more statistically significant interactions compared to end-point assays and further observed several examples of context-dependent crosstalk between signaling pathways such as an interaction between Ras and Rel which is dependent on MEK activity. Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

An Automated Image Analysis Pipeline Enables Genetic Studies of Shoot and Root Morphology in Carrot (Daucus carota L.).

Carrot is a globally important crop, yet efficient and accurate methods for quantifying its most important agronomic traits are lacking. To address this problem, we developed an automated image analysis platform that extracts components of size and shape for carrot shoots and roots, which are necessary to advance carrot breeding and genetics. This method reliably measured variation in shoot size and shape, petiole number, petiole length, and petiole width as evidenced by high correlations with hundreds of manual measurements. Similarly, root length and biomass were accurately measured from the images. This platform also quantified shoot and root shapes in terms of principal components, which do not have traditional, manually measurable equivalents. We applied the pipeline in a study of a six-parent diallel population and an F2 mapping population consisting of 316 individuals. We found high levels of repeatability within a growing environment, with low to moderate repeatability across environments. We also observed co-localization of quantitative trait loci for shoot and root characteristics on chromosomes 1, 2, and 7, suggesting these traits are controlled by genetic linkage and/or pleiotropy. By increasing the number of individuals and phenotypes that can be reliably quantified, the development of a rapid, automated image analysis pipeline to measure carrot shoot and root morphology will expand the scope and scale of breeding and genetic studies.

Genomic and Phenomic Screens for Flower Related RING Type Ubiquitin E3 Ligases in Arabidopsis.

Flowering time control integrates endogenous as well as environmental signals to promote flower development. The pathways and molecular networks involved are complex and integrate many modes of signal transduction. In plants ubiquitin mediated protein degradation pathway has been proposed to be as important mode of signaling as phosphorylation and transcription. To systematically study the role of ubiquitin signaling in the molecular regulation of flowering we have taken a genomic approach to identify flower related Ubiquitin Proteasome System components. As a large and versatile gene family the RING type ubiquitin E3 ligases were chosen as targets of the genomic screen. The complete list of Arabidopsis RING E3 ligases were retrieved and verified in the Arabidopsis genome v11 and their differential expression was used for their categorization into flower organs or developmental stages. Known regulators of flowering time or floral organ development were identified in these categories through literature search and representative mutants for each category were purchased for functional characterization by growth and morphological phenotyping. To this end, a workflow was developed for high throughput phenotypic screening of growth, morphology and flowering of nearly a thousand Arabidopsis plants in one experimental round.

A map of directional genetic interactions in a metazoan cell.

Gene-gene interactions shape complex phenotypes and modify the effects of mutations during development and disease. The effects of statistical gene-gene interactions on phenotypes have been used to assign genes to functional modules. However, directional, epistatic interactions, which reflect regulatory relationships between genes, have been challenging to map at large-scale. Here, we used combinatorial RNA interference and automated single-cell phenotyping to generate a large genetic interaction map for 21 phenotypic features of Drosophila cells. We devised a method that combines genetic interactions on multiple phenotypes to reveal directional relationships. This network reconstructed the sequence of protein activities in mitosis. Moreover, it revealed that the Ras pathway interacts with the SWI/SNF chromatin-remodelling complex, an interaction that we show is conserved in human cancer cells. Our study presents a powerful approach for reconstructing directional regulatory networks and provides a resource for the interpretation of functional consequences of genetic alterations.