Nothing Special   »   [go: up one dir, main page]

You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Sign in to use this feature.

Years

Between: -

Subjects

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Journals

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Article Types

Countries / Regions

remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline
remove_circle_outline

Search Results (880)

Search Parameters:
Keywords = high-throughput phenotyping

Order results
Result details
Results per page
Select all
Export citation of selected articles as:
19 pages, 4707 KiB  
Article
Chlorophyll Content Estimation of Ginkgo Seedlings Based on Deep Learning and Hyperspectral Imagery
by Zilong Yue, Qilin Zhang, Xingzhou Zhu and Kai Zhou
Forests 2024, 15(11), 2010; https://doi.org/10.3390/f15112010 - 14 Nov 2024
Viewed by 442
Abstract
Accurate estimation of chlorophyll content is essential for understanding the growth status and optimizing the cultivation practices of Ginkgo, a dominant multi-functional tree species in China. Traditional methods based on chemical analysis for determining chlorophyll content are labor-intensive and time-consuming, making them [...] Read more.
Accurate estimation of chlorophyll content is essential for understanding the growth status and optimizing the cultivation practices of Ginkgo, a dominant multi-functional tree species in China. Traditional methods based on chemical analysis for determining chlorophyll content are labor-intensive and time-consuming, making them unsuitable for large-scale dynamic monitoring and high-throughput phenotyping. To accurately quantify chlorophyll content in Ginkgo seedlings under different nitrogen levels, this study employed a hyperspectral imaging camera to capture canopy hyperspectral images of seedlings throughout their annual growth periods. Reflectance derived from pure leaf pixels of Ginkgo seedlings was extracted to construct a set of spectral parameters, including original reflectance, logarithmic reflectance, and first derivative reflectance, along with spectral index combinations. A one-dimensional convolutional neural network (1D-CNN) model was then developed to estimate chlorophyll content, and its performance was compared with four common machine learning methods, including Gaussian Process Regression (GPR), Partial Least Squares Regression (PLSR), Support Vector Regression (SVR), and Random Forest (RF). The results demonstrated that the 1D-CNN model outperformed others with the first derivative spectra, achieving higher CV-R2 and lower RMSE values (CV-R2 = 0.80, RMSE = 3.4). Furthermore, incorporating spectral index combinations enhanced the model’s performance, with the 1D-CNN model achieving the best performance (CV-R2 = 0.82, RMSE = 3.3). These findings highlight the potential of the 1D-CNN model in strengthening the chlorophyll estimations, providing strong technical support for the precise cultivation and the fertilization management of Ginkgo seedlings. Full article
(This article belongs to the Section Forest Inventory, Modeling and Remote Sensing)
Show Figures

Figure 1

Figure 1
<p>The workflow for estimating chlorophyll content in <span class="html-italic">Ginkgo</span> canopies based on hyperspectral imaging and 1D-CNN.</p>
Full article ">Figure 2
<p>Schematic representation of the experimental layout for <span class="html-italic">Ginkgo biloba</span> seedlings under five nitrogen treatments (N0–N4). Each treatment was replicated three times (R1–R3), resulting in 15 experimental units in total. Nitrogen was applied as a topdressing in three equal doses during the growing season.</p>
Full article ">Figure 3
<p>Hyperspectral images of <span class="html-italic">Ginkgo biloba</span> seedlings under different nitrogen levels (N0–N4) across growth stages (T1–T5). T1 corresponds to April (early bud development stage), T2 corresponds to May (early rapid growth stage), T3 corresponds to June (middle rapid growth stage), T4 corresponds to July (late rapid growth stage), and T5 corresponds to August (plant maturity stage).</p>
Full article ">Figure 4
<p>A suitable 1D-CNN model for spectral reflectance is proposed.</p>
Full article ">Figure 5
<p>Changes in canopy reflectance and SPAD chlorophyll content of <span class="html-italic">Ginkgo</span> seedlings across different growth stages and nitrogen fertilization levels (<b>A</b>) Spectral reflectance curves of the <span class="html-italic">Ginkgo</span> canopy at different SPAD levels. The figure includes three forms of reflectance spectra: (i) original reflectance, (ii) logarithmic reflectance, and (iii) first derivative reflectance. SPAD chlorophyll content is divided into low (SPAD_low), medium (SPAD_medium), and high (SPAD_high) levels. Low SPAD corresponds to values from 27 to 45, medium SPAD ranges between 45 and 55, and high SPAD corresponds to values from 55 to 65. Reflectance across the 400 to 900 nm range varies with SPAD levels, reflecting the sensitivity of different spectral regions to chlorophyll absorption and canopy structure. (<b>B</b>) Changes in SPAD chlorophyll content of <span class="html-italic">Ginkgo</span> seedlings at different growth stages (T1–T5). T1 represents April (early bud development stage), T2 represents May (early rapid growth stage), T3 represents June (middle rapid growth stage), T4 represents July (late rapid growth stage), and T5 represents August (plant maturity stage). SPAD content fluctuates across the different growth stages. (<b>C</b>) Changes in SPAD chlorophyll content of <span class="html-italic">Ginkgo</span> seedlings under different nitrogen fertilization treatments (N0–N4). SPAD content shows significant variation across the different nitrogen levels.</p>
Full article ">Figure 6
<p>Correlation coefficient curve between leaf SPAD-chlorophyll content in <span class="html-italic">Ginkgo</span> seedlings and original or transformed reflectance spectra.</p>
Full article ">Figure 7
<p>The correlation between leaf chlorophyll content and the three best-performing indices: SR<sub>708,775</sub>, GNDVI, and mCI<sub>Green</sub>.</p>
Full article ">Figure 8
<p>Optimal reflectance datasets for the correlation of DVI<sub>log</sub> ((<b>A</b>), logarithmic reflectance), RVI<sub>FD</sub> ((<b>B</b>), first derivative of reflectance), NDVI<sub>FD</sub> ((<b>C</b>), first derivative of reflectance), and mRVI<sub>log</sub> ((<b>D</b>), logarithmic reflectance) with chlorophyll content. The white arrow indicates the optimal band combination.</p>
Full article ">Figure 9
<p>Comparison of predicted versus measured SPAD values for <span class="html-italic">Ginkgo</span> seedling canopies using various modeling approaches. Best Spectrum-Orgi, Best Spectrum-log and Best Spectrum-FD represent the best-performing spectral data (Orgi: original spectra; log: logarithmic spectra; FD: first-derivative spectra) for each regression method. VI represents vegetation indices.</p>
Full article ">
22 pages, 5474 KiB  
Article
Comparative Transcriptome Analysis of Sexual Differentiation in Male and Female Gonads of Nao-Zhou Stock Large Yellow Croaker (Larimichthys crocea)
by Haojie Wang, Zirui Wen, Eric Amenyogbe, Jinghui Jin, Yi Lu, Zhongliang Wang and Jiansheng Huang
Animals 2024, 14(22), 3261; https://doi.org/10.3390/ani14223261 - 13 Nov 2024
Viewed by 427
Abstract
The Nao-zhou stock large yellow croaker (Larimichthys crocea) is a unique economic seawater fish species in China and exhibits significant dimorphism in both male and female phenotypes. Cultivating all-female seedlings can significantly improve breeding efficiency. To accelerate the cultivation process of [...] Read more.
The Nao-zhou stock large yellow croaker (Larimichthys crocea) is a unique economic seawater fish species in China and exhibits significant dimorphism in both male and female phenotypes. Cultivating all-female seedlings can significantly improve breeding efficiency. To accelerate the cultivation process of all female seedlings of this species, it is necessary to deeply understand the regulatory mechanisms of sexual differentiation and gonadal development. This study used Illumina high-throughput sequencing to sequence the transcriptome of the testes and ovaries of Nao-zhou stock large yellow croaker to identify genes and molecular functions related to sex determination. A total of 10,536 differentially expressed genes were identified between males and females, including 5682 upregulated and 4854 downregulated genes. Functional annotation screened out 70 important candidate genes related to sex, including 34 genes highly expressed in the testis (including dmrt1, foxm1, and amh) and 36 genes highly expressed in the ovary (including gdf9, hsd3b1, and sox19b). Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis found that differentially expressed genes were significantly enriched in nine signaling pathways related to sex determination and gonadal development, including steroid hormone biosynthesis, MAPK signaling pathway, and the TGF-beta signaling pathway. By screening sex-related differentially expressed genes and mapping protein–protein interaction networks, hub genes such as dmrt1, amh, and cyp19a1a were found to be highly connected. The expression levels of 15 sex-related genes, including amh, dmrt1, dmrt2a, foxl1, and zp3b, were determined by qRT–PCR and RNA sequencing. This study screened for differentially expressed genes related to sex determination and differentiation of Nao-zhou stock large yellow croaker and revealed the signaling pathways involved in gonad development of male and female individuals. The results provide important data for future research on sex determination and differentiation mechanisms, thereby providing a scientific basis for the cultivation of all-female seedlings. Full article
(This article belongs to the Section Animal Physiology)
Show Figures

Figure 1

Figure 1
<p>Histological characteristics of testes (<b>a</b>) and ovaries (<b>b</b>) of Nao-zhou stock large yellow croaker. Notes: Sp: sperm; Spe: sperm cell; Sl: sperm lobule; Yg: yolk granule; N: nucleus; Yv: yolk vesicle; Nu: Nucleolus.</p>
Full article ">Figure 2
<p>Volcano map of differentially expressed genes in Nao-zhou stock large yellow croaker. Note: The horizontal axis shows the log<sub>2</sub> value (fold change), the vertical axis is the −log<sub>10</sub> value (<span class="html-italic">p</span> value), green dots represent upregulated genes, red dots represent downregulated genes, and blue dots represent genes with no significance. The dotted lines represent the threshold of log<sub>2</sub>(FC) values.</p>
Full article ">Figure 3
<p>Violin plot and cluster heat map of 6 samples. Note: (<b>a</b>) represents the correlation of samples between and within groups. (<b>b</b>) shows cluster results of DEGs. The color indicates the expression amount (logarithm) or the difference multiple (logarithm). The redder color indicates that the gene expression level is higher or the difference factor is larger, and the blue color indicates the opposite.</p>
Full article ">Figure 4
<p>Top 30 GO enrichment pathways of differentially expressed genes in the gonads of Nao-zhou stock large yellow croaker. Note: The horizontal axis shows the gene name, and the vertical axis shows the gene ratio.</p>
Full article ">Figure 5
<p>Top 30 KEGG enrichment pathways of differentially expressed genes in the gonads of Nao-zhou stock large yellow croaker. Note: The horizontal axis shows the gene name, and the vertical axis shows the gene ratio.</p>
Full article ">Figure 6
<p>GO (<b>a</b>) and KEGG (<b>b</b>) enriched pathways of the top 20 differentially expressed genes associated with sex in Nao-zhou stock large yellow croaker.</p>
Full article ">Figure 7
<p>Protein-protein interaction (PPI) network diagram of DEGs in female and male Nao-zhou stock large yellow croaker. Note: Different background colors represent the network degree values of proteins. The inner circle of the PPI network shows hub genes, while the outer two circles are non-hub genes. Number of gene nodes is represented by color depth.</p>
Full article ">Figure 8
<p>Relative expression levels of 15 genes in the testis and ovary of Nao-zhou stock large yellow croaker. Note: Data are presented as mean ± S.E.M. (n = 3). The asterisks indicate that the differences between the mean values are statistically significant between gonads. *: 0.01 &lt; <span class="html-italic">p</span> &lt; 0.05; **: 0.001 &lt; <span class="html-italic">p</span> &lt; 0.01; ***: <span class="html-italic">p</span> &lt; 0.001.</p>
Full article ">Figure 9
<p>qRT-PCR verification of sex-related differentially expressed genes. Note: The horizontal axis shows the gene name, and the vertical axis shows the relative expression level.</p>
Full article ">Figure 10
<p>Chord diagram of the functional classification of twelve candidate genes. Note: The left half represents candidate genes and expression levels, and the right half represents GO enriched pathways related to reproduction.</p>
Full article ">
21 pages, 3828 KiB  
Article
Clinical Phenotypes Associated with the Gut Microbiome in Older Japanese People with Care Needs in a Nursing Home
by Rikako Inoue, Koji Hosomi, Jonguk Park, Haruka Sakaue, Hitomi Yumioka, Hiroko Kamitani, Yoshiharu Kinugasa, Kaori Harano, A. Yasmin Syauki, Miki Doi, Suzumi Kageyama, Kazuhiro Yamamoto, Kenji Mizuguchi, Jun Kunisawa and Yasuyuki Irie
Nutrients 2024, 16(22), 3839; https://doi.org/10.3390/nu16223839 - 8 Nov 2024
Viewed by 465
Abstract
Background: Frailty increases the risk of needing nursing care and significantly affects the life and functional prognosis of older individuals. Early detection and tailored interventions are crucial for maintaining and enhancing their life functions. Recognizing distinct clinical phenotypes is essential for devising appropriate [...] Read more.
Background: Frailty increases the risk of needing nursing care and significantly affects the life and functional prognosis of older individuals. Early detection and tailored interventions are crucial for maintaining and enhancing their life functions. Recognizing distinct clinical phenotypes is essential for devising appropriate interventions. This study aimed to explore diverse frailty phenotypes, focusing on poor nutrition in older Japanese individuals through observational research. Methods: Twenty-one nursing home residents underwent a comprehensive survey covering physical, blood, dietary, cardiac, cognitive, nutritional, nursing care, frailty, agitated behavior, and gut microbiome assessments (high-throughput 16S rRNA gene sequencing). Using clustering analysis with 239 survey items (excluding gut microbiome), participants were classified into subgroups based on clinical phenotypes, and group characteristics were compared through analysis. Results: Individuals with moderate or severe frailty and suspected dementia formed subgroups with distinct clinical phenotypes based on nutritional, defecation, and nursing care statuses. The gut microbiome significantly varied among these groups (p = 0.007), indicating its correlation with changes in clinical phenotype. Nutritional status differences suggested poor nutrition as a differentiating factor in the core clinical phenotype. Conclusions: This study proposes that the gut microbiome differs based on the clinical phenotype of Japanese older individuals with frailty, and targeted interventions addressing the gut microbiome may contribute to preventing frailty in this population. Full article
(This article belongs to the Section Geriatric Nutrition)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>The clinical phenotypes of frail older individuals. (<b>A</b>) Principal component analysis (PCA): The PCA factor map shows four distinct clusters based on participant characteristics. Cluster 1 (green circles) represents similar traits, Cluster 2 (orange triangles) indicates a different subgroup, Cluster 3 (blue squares) highlights another unique group, and Cluster 4 (pink diamonds) depicts a separate cluster. Dim1 accounts for 21.6% of the variance, while Dim2 accounts for 12.7%. Each point corresponds to an individual clustered by analyzed features. (<b>B</b>) Hierarchical clustering on principal components (HCPC) and a comparative analysis: characterization the clinical phenotype in each group; we listed the items with <span class="html-italic">p</span> &lt; 0.01 in a comparative analysis, excluding the questions in each survey and various nutrients.</p>
Full article ">Figure 2
<p>Intergroup comparison of items characterizing the clinical phenotype of each group. (<b>A</b>) Nutritional status: (a) MNA, (b) BMI, (c) eGFR (creatinine), (d) total carnitine. (<b>B</b>) Defecation: (a) toilet defecation, (b) excretions on pad, (c) diarrhea, (d) BSS. (<b>C</b>) Care needs: (a) Clinical Frailty Scale (CFS), (b) care level (point). Data are shown as the mean ± SEM. Comparison between groups: Mann–Whitney U test. ** <span class="html-italic">p</span> &lt; 0.01, * <span class="html-italic">p</span> &lt; 0.05, <sup>†</sup> <span class="html-italic">p</span> &lt; 0.1.</p>
Full article ">Figure 3
<p>Correlation matrix heatmap between items characterizing the clinical phenotype of each group, as determined by a comparative analysis. (<b>a</b>) Correlation analysis (heatmap): Spearman’s correlation: ** <span class="html-italic">p</span> &lt; 0.01, * <span class="html-italic">p</span> &lt; 0.05. (<b>b</b>) Correlation analysis (plot diagram): plot diagrams are shown for items with <span class="html-italic">p</span> &lt; 0.05 according to Spearman’s correlation. (A) eGFRcreat vs. BMI (<span class="html-italic">p</span> = 0.002) (B) BMI vs. MNA (<span class="html-italic">p</span> &lt; 0.001) (C) eGFRcreat vs. total carnitine (<span class="html-italic">p</span> = 0.006) (D) eGFRcreat vs. free carnitine (<span class="html-italic">p</span> = 0.006) (E) eGFRcreat vs. acyl carnitine (<span class="html-italic">p</span> = 0.030) (F) eGFRcreat vs. toilets detected (<span class="html-italic">p</span> = 0.043) (G) eGFRcreat vs. pats excreted (<span class="html-italic">p</span> = 0.035) (H) eGFRcreat vs. BSS (<span class="html-italic">p</span> = 0.039) (I) eGFRcreat vs. frailty (<span class="html-italic">p</span> = 0.047) (J) eGFRcreat vs. MNA (<span class="html-italic">p</span> &lt; 0.001) (K) free carnitine vs. acyl carnitine (<span class="html-italic">p</span> &lt; 0.001) (L) total carnitine vs. acyl carnitine (<span class="html-italic">p</span> &lt; 0.001) (M) free carnitine vs. pats excreted (<span class="html-italic">p</span> &lt; 0.001) (N) toilets defecated vs. pats excreted (<span class="html-italic">p</span> = 0.016) (O) toilets defecated vs. diarrhea (<span class="html-italic">p</span> = 0.004) (P) toilets defecated vs. BSS (<span class="html-italic">p</span> = 0.010) (Q) toilets defecated vs. frailty (<span class="html-italic">p</span> = 0.016) (R) toilets defecated vs. MNA (<span class="html-italic">p</span> = 0.011) (S) Care level vs. toilets defeated (<span class="html-italic">p</span> = 0.042) (T) Diarrhea vs. pats excreted (<span class="html-italic">p</span> = 0.003) (U) BSS vs. pats excreted (<span class="html-italic">p</span> = 0.007) (V) Frailty vs. pats excreted (<span class="html-italic">p</span> = 0.007) (W) MNA vs. pats excreted (<span class="html-italic">p</span> = 0.001) (X) Care level vs. pats excreted (<span class="html-italic">p</span> = 0.024) (Y) BSS vs. diarrhea (<span class="html-italic">p</span> = 0.023) (Z) MNA vs. diarrhea (<span class="html-italic">p</span> = 0.006) (AA) Care level vs. diarrhea (<span class="html-italic">p</span> = 0.010) (AB) MNA vs. BSS (<span class="html-italic">p</span> &lt; 0.001) (AC) MNA vs. frailty (<span class="html-italic">p</span> = 0.025) (AD) Care level vs. frailty (<span class="html-italic">p</span> = 0.006). (<b>c</b>) Correlation matrix heatmap between nutrients and items characterizing the clinical phenotype. Spearman’s correlation: ** <span class="html-italic">p</span> &lt; 0.01, * <span class="html-italic">p</span> &lt; 0.05.</p>
Full article ">Figure 3 Cont.
<p>Correlation matrix heatmap between items characterizing the clinical phenotype of each group, as determined by a comparative analysis. (<b>a</b>) Correlation analysis (heatmap): Spearman’s correlation: ** <span class="html-italic">p</span> &lt; 0.01, * <span class="html-italic">p</span> &lt; 0.05. (<b>b</b>) Correlation analysis (plot diagram): plot diagrams are shown for items with <span class="html-italic">p</span> &lt; 0.05 according to Spearman’s correlation. (A) eGFRcreat vs. BMI (<span class="html-italic">p</span> = 0.002) (B) BMI vs. MNA (<span class="html-italic">p</span> &lt; 0.001) (C) eGFRcreat vs. total carnitine (<span class="html-italic">p</span> = 0.006) (D) eGFRcreat vs. free carnitine (<span class="html-italic">p</span> = 0.006) (E) eGFRcreat vs. acyl carnitine (<span class="html-italic">p</span> = 0.030) (F) eGFRcreat vs. toilets detected (<span class="html-italic">p</span> = 0.043) (G) eGFRcreat vs. pats excreted (<span class="html-italic">p</span> = 0.035) (H) eGFRcreat vs. BSS (<span class="html-italic">p</span> = 0.039) (I) eGFRcreat vs. frailty (<span class="html-italic">p</span> = 0.047) (J) eGFRcreat vs. MNA (<span class="html-italic">p</span> &lt; 0.001) (K) free carnitine vs. acyl carnitine (<span class="html-italic">p</span> &lt; 0.001) (L) total carnitine vs. acyl carnitine (<span class="html-italic">p</span> &lt; 0.001) (M) free carnitine vs. pats excreted (<span class="html-italic">p</span> &lt; 0.001) (N) toilets defecated vs. pats excreted (<span class="html-italic">p</span> = 0.016) (O) toilets defecated vs. diarrhea (<span class="html-italic">p</span> = 0.004) (P) toilets defecated vs. BSS (<span class="html-italic">p</span> = 0.010) (Q) toilets defecated vs. frailty (<span class="html-italic">p</span> = 0.016) (R) toilets defecated vs. MNA (<span class="html-italic">p</span> = 0.011) (S) Care level vs. toilets defeated (<span class="html-italic">p</span> = 0.042) (T) Diarrhea vs. pats excreted (<span class="html-italic">p</span> = 0.003) (U) BSS vs. pats excreted (<span class="html-italic">p</span> = 0.007) (V) Frailty vs. pats excreted (<span class="html-italic">p</span> = 0.007) (W) MNA vs. pats excreted (<span class="html-italic">p</span> = 0.001) (X) Care level vs. pats excreted (<span class="html-italic">p</span> = 0.024) (Y) BSS vs. diarrhea (<span class="html-italic">p</span> = 0.023) (Z) MNA vs. diarrhea (<span class="html-italic">p</span> = 0.006) (AA) Care level vs. diarrhea (<span class="html-italic">p</span> = 0.010) (AB) MNA vs. BSS (<span class="html-italic">p</span> &lt; 0.001) (AC) MNA vs. frailty (<span class="html-italic">p</span> = 0.025) (AD) Care level vs. frailty (<span class="html-italic">p</span> = 0.006). (<b>c</b>) Correlation matrix heatmap between nutrients and items characterizing the clinical phenotype. Spearman’s correlation: ** <span class="html-italic">p</span> &lt; 0.01, * <span class="html-italic">p</span> &lt; 0.05.</p>
Full article ">Figure 4
<p>Characteristics of the intestinal microbiota of each group with different clinical phenotypes. (<b>A</b>) Principal coordinate analysis (PCoA) of the intestinal microbiota at the genus level. (<b>B</b>) Group comparison at the genus (D:5) or family (D:4) level. Data are shown as the mean ± SEM. Comparison between groups was performed using the Mann–Whitney U test; ** <span class="html-italic">p</span> &lt; 0.01, * <span class="html-italic">p</span> &lt; 0.05, <sup>†</sup> <span class="html-italic">p</span> &lt; 0.1.</p>
Full article ">Figure 4 Cont.
<p>Characteristics of the intestinal microbiota of each group with different clinical phenotypes. (<b>A</b>) Principal coordinate analysis (PCoA) of the intestinal microbiota at the genus level. (<b>B</b>) Group comparison at the genus (D:5) or family (D:4) level. Data are shown as the mean ± SEM. Comparison between groups was performed using the Mann–Whitney U test; ** <span class="html-italic">p</span> &lt; 0.01, * <span class="html-italic">p</span> &lt; 0.05, <sup>†</sup> <span class="html-italic">p</span> &lt; 0.1.</p>
Full article ">Figure 5
<p>Correlation analysis between gut microbiome and items characterizing the clinical phenotype. (<b>a</b>) Correlation analysis (heatmap): Spearman’s correlation: ** <span class="html-italic">p</span> &lt; 0.01, * <span class="html-italic">p</span> &lt; 0.05. (<b>b</b>) Correlation analysis (plot diagram): plot diagrams are shown for items with <span class="html-italic">p</span> &lt; 0.05 according to Spearman’s correlation. (A) <span class="html-italic">Acidaminococcaceae</span> vs. pats excreted (<span class="html-italic">p</span> = 0.014) (B) <span class="html-italic">Acidaminococcaceae</span> vs. diarrhea (<span class="html-italic">p</span> = 0.050) (C) <span class="html-italic">Acidaminococcaceae</span> vs. frailty (<span class="html-italic">p</span> = 0.010) (D) <span class="html-italic">Acidaminococcaceae</span> vs. MNA (<span class="html-italic">p</span> = 0.003) (E) <span class="html-italic">Acidaminococcaceae</span> vs. care level (<span class="html-italic">p</span> = 0.005) (F) <span class="html-italic">Ruminococcaceae</span> vs. care level (<span class="html-italic">p</span> = 0.045) (G) <span class="html-italic">Porphyromonadaceae</span> vs. pats excreted (<span class="html-italic">p</span> = 0.047) (H) <span class="html-italic">Porphyromonadaceae</span> vs. diarrhea (<span class="html-italic">p</span> = 0.006) (I) <span class="html-italic">Porphyromonadaceae</span> vs. MNA (<span class="html-italic">p</span> = 0.019) (J) <span class="html-italic">Klebsiella</span> vs. total carnitine (<span class="html-italic">p</span> = 0.040) (K) <span class="html-italic">Klebsiella</span> vs. free carnitine (<span class="html-italic">p</span> = 0.034) (L) <span class="html-italic">Klebsiella</span> vs. BSS (<span class="html-italic">p</span> = 0.019) (M) <span class="html-italic">Phascolarctobacterium</span> vs. pats excreted (<span class="html-italic">p</span> = 0.037) (N) <span class="html-italic">Phascolarctobacterium</span> vs. frailty (<span class="html-italic">p</span> = 0.012) (O) <span class="html-italic">Phascolarctobacterium</span> vs. MNA (<span class="html-italic">p</span> = 0.016). (P) <span class="html-italic">Phascolarctobacterium</span> vs. care level (<span class="html-italic">p</span> = 0.012) (Q) <span class="html-italic">Roseburia</span> vs. eGFR (<span class="html-italic">p</span> = 0.008) (R) <span class="html-italic">Roseburia</span> vs. free carnitine (<span class="html-italic">p</span> = 0.041) (S) <span class="html-italic">Roseburia</span> vs. pats excreted (<span class="html-italic">p</span> = 0.048) (T) <span class="html-italic">Roseburia</span> vs. MNA (<span class="html-italic">p</span> = 0.027) (U) <span class="html-italic">Parabacteriodes</span> vs. toilets detected (<span class="html-italic">p</span> = 0.045) (V) <span class="html-italic">Parabacteriodes</span> vs. pats excreted (<span class="html-italic">p</span> = 0.025) (W) <span class="html-italic">Parabacteriodes</span> vs. diarrhea (<span class="html-italic">p</span> = 0.003) (X) <span class="html-italic">Parabacteriodes</span> vs. BSS (<span class="html-italic">p</span> = 0.019) (Y) <span class="html-italic">Parabacteriodes</span> vs. MNA (<span class="html-italic">p</span> = 0.004) (Z) <span class="html-italic">torques group</span> vs. pats excreted (<span class="html-italic">p</span> = 0.044) (AA) <span class="html-italic">torques group</span> vs. diarrhea (<span class="html-italic">p</span> &lt; 0.001) (AB) <span class="html-italic">torques group</span> vs. MNA (<span class="html-italic">p</span> = 0.022); Spearman’s correlation: ** <span class="html-italic">p</span> &lt; 0.01, * <span class="html-italic">p</span> &lt; 0.05.</p>
Full article ">Figure 5 Cont.
<p>Correlation analysis between gut microbiome and items characterizing the clinical phenotype. (<b>a</b>) Correlation analysis (heatmap): Spearman’s correlation: ** <span class="html-italic">p</span> &lt; 0.01, * <span class="html-italic">p</span> &lt; 0.05. (<b>b</b>) Correlation analysis (plot diagram): plot diagrams are shown for items with <span class="html-italic">p</span> &lt; 0.05 according to Spearman’s correlation. (A) <span class="html-italic">Acidaminococcaceae</span> vs. pats excreted (<span class="html-italic">p</span> = 0.014) (B) <span class="html-italic">Acidaminococcaceae</span> vs. diarrhea (<span class="html-italic">p</span> = 0.050) (C) <span class="html-italic">Acidaminococcaceae</span> vs. frailty (<span class="html-italic">p</span> = 0.010) (D) <span class="html-italic">Acidaminococcaceae</span> vs. MNA (<span class="html-italic">p</span> = 0.003) (E) <span class="html-italic">Acidaminococcaceae</span> vs. care level (<span class="html-italic">p</span> = 0.005) (F) <span class="html-italic">Ruminococcaceae</span> vs. care level (<span class="html-italic">p</span> = 0.045) (G) <span class="html-italic">Porphyromonadaceae</span> vs. pats excreted (<span class="html-italic">p</span> = 0.047) (H) <span class="html-italic">Porphyromonadaceae</span> vs. diarrhea (<span class="html-italic">p</span> = 0.006) (I) <span class="html-italic">Porphyromonadaceae</span> vs. MNA (<span class="html-italic">p</span> = 0.019) (J) <span class="html-italic">Klebsiella</span> vs. total carnitine (<span class="html-italic">p</span> = 0.040) (K) <span class="html-italic">Klebsiella</span> vs. free carnitine (<span class="html-italic">p</span> = 0.034) (L) <span class="html-italic">Klebsiella</span> vs. BSS (<span class="html-italic">p</span> = 0.019) (M) <span class="html-italic">Phascolarctobacterium</span> vs. pats excreted (<span class="html-italic">p</span> = 0.037) (N) <span class="html-italic">Phascolarctobacterium</span> vs. frailty (<span class="html-italic">p</span> = 0.012) (O) <span class="html-italic">Phascolarctobacterium</span> vs. MNA (<span class="html-italic">p</span> = 0.016). (P) <span class="html-italic">Phascolarctobacterium</span> vs. care level (<span class="html-italic">p</span> = 0.012) (Q) <span class="html-italic">Roseburia</span> vs. eGFR (<span class="html-italic">p</span> = 0.008) (R) <span class="html-italic">Roseburia</span> vs. free carnitine (<span class="html-italic">p</span> = 0.041) (S) <span class="html-italic">Roseburia</span> vs. pats excreted (<span class="html-italic">p</span> = 0.048) (T) <span class="html-italic">Roseburia</span> vs. MNA (<span class="html-italic">p</span> = 0.027) (U) <span class="html-italic">Parabacteriodes</span> vs. toilets detected (<span class="html-italic">p</span> = 0.045) (V) <span class="html-italic">Parabacteriodes</span> vs. pats excreted (<span class="html-italic">p</span> = 0.025) (W) <span class="html-italic">Parabacteriodes</span> vs. diarrhea (<span class="html-italic">p</span> = 0.003) (X) <span class="html-italic">Parabacteriodes</span> vs. BSS (<span class="html-italic">p</span> = 0.019) (Y) <span class="html-italic">Parabacteriodes</span> vs. MNA (<span class="html-italic">p</span> = 0.004) (Z) <span class="html-italic">torques group</span> vs. pats excreted (<span class="html-italic">p</span> = 0.044) (AA) <span class="html-italic">torques group</span> vs. diarrhea (<span class="html-italic">p</span> &lt; 0.001) (AB) <span class="html-italic">torques group</span> vs. MNA (<span class="html-italic">p</span> = 0.022); Spearman’s correlation: ** <span class="html-italic">p</span> &lt; 0.01, * <span class="html-italic">p</span> &lt; 0.05.</p>
Full article ">
21 pages, 5851 KiB  
Article
SAM-ResNet50: A Deep Learning Model for the Identification and Classification of Drought Stress in the Seedling Stage of Betula luminifera
by Shiya Gao, Hao Liang, Dong Hu, Xiange Hu, Erpei Lin and Huahong Huang
Remote Sens. 2024, 16(22), 4141; https://doi.org/10.3390/rs16224141 - 6 Nov 2024
Viewed by 417
Abstract
Betula luminifera, an indigenous hardwood tree in South China, possesses significant economic and ecological value. In view of the current severe drought situation, it is urgent to enhance this tree’s drought tolerance. However, traditional artificial methods fall short of meeting the demands [...] Read more.
Betula luminifera, an indigenous hardwood tree in South China, possesses significant economic and ecological value. In view of the current severe drought situation, it is urgent to enhance this tree’s drought tolerance. However, traditional artificial methods fall short of meeting the demands of breeding efforts due to their inefficiency. To monitor drought situations in a high-throughput and automatic approach, a deep learning model based on phenotype characteristics was proposed to identify and classify drought stress in B. luminifera seedlings. Firstly, visible-light images were obtained from a drought stress experiment conducted on B. luminifera shoots. Considering the images’ characteristics, we proposed an SAM-CNN architecture by incorporating spatial attention modules into classical CNN models. Among the four classical CNNs compared, ResNet50 exhibited superior performance and was, thus, selected for the construction of the SAM-CNN. Subsequently, we analyzed the classification performance of the SAM-ResNet50 model in terms of transfer learning, training from scratch, model robustness, and visualization. The results revealed that SAM-ResNet50 achieved an accuracy of 1.48% higher than that of ResNet50, at 99.6%. Furthermore, there was a remarkable improvement of 18.98% in accuracy, reaching 82.31% for the spatial transform images generated from the test set images by applying movement and rotation for robustness testing. In conclusion, the SAM-ResNet50 model achieved outstanding performance, with 99.6% accuracy and realized high-throughput automatic monitoring based on phenotype, providing a new perspective for drought stress classification and technical support for B. luminifera-related breeding work. Full article
(This article belongs to the Special Issue Intelligent Extraction of Phenotypic Traits in Agroforestry)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Image acquisition platform for <span class="html-italic">Betula luminifera</span>.</p>
Full article ">Figure 2
<p>Sample images from the <span class="html-italic">B. luminifera</span> drought stress dataset (<b>a</b>–<b>f</b>) and from the moved and rotated datasets (<b>g</b>–<b>l</b>): (<b>a</b>) CG; (<b>b</b>) DR1; (<b>c</b>) DR2; (<b>d</b>) DR3; (<b>e</b>) DR4; (<b>f</b>) DR5; (<b>g</b>) moved down one-fifth to the left; (<b>h</b>) moved up one-sixth to the left; (<b>i</b>) moved down one-fifth to the right; (<b>j</b>) moved up one-sixth to the right; (<b>k</b>) turned 45° clockwise; (<b>l</b>) turned 45° counterclockwise.</p>
Full article ">Figure 3
<p>Workflow diagram.</p>
Full article ">Figure 4
<p>The SAM-CNN architecture.</p>
Full article ">Figure 5
<p>The change curves of the loss function value and accuracy during training: (<b>a</b>) loss function; (<b>b</b>) accuracy.</p>
Full article ">Figure 6
<p>The confusion matrixes describing the results of different models: (<b>a</b>) AlexNet; (<b>b</b>) VGG16; (<b>c</b>) GoogLeNet; (<b>d</b>) ResNet50; (<b>e</b>) SAM-ResNet50.</p>
Full article ">Figure 7
<p>The change curves of the loss function value and accuracy during training: (<b>a</b>) loss function; (<b>b</b>) accuracy.</p>
Full article ">Figure 8
<p>The change curves of the loss function value and accuracy during training: (<b>a</b>) loss function; (<b>b</b>) accuracy.</p>
Full article ">Figure 9
<p>The confusion matrixes describing the results of different models: (<b>a</b>) ResNet50; (<b>b</b>) SAM-ResNet50.</p>
Full article ">Figure 10
<p>Comparison of the visualization results before and after adding the SAM. The highlighted part of the class activation map represents the network’s attention on the image, with warmer colors indicating higher attention and cooler colors indicating lower attention: (<b>a</b>) the original image; (<b>b</b>) the area of attention before adding the attention mechanism; (<b>c</b>) the area of attention after adding the attention mechanism.</p>
Full article ">
17 pages, 975 KiB  
Review
Drosophila as a Model for Human Disease: Insights into Rare and Ultra-Rare Diseases
by Sergio Casas-Tintó
Insects 2024, 15(11), 870; https://doi.org/10.3390/insects15110870 - 6 Nov 2024
Viewed by 866
Abstract
Rare and ultra-rare diseases constitute a significant medical challenge due to their low prevalence and the limited understanding of their origin and underlying mechanisms. These disorders often exhibit phenotypic diversity and molecular complexity that represent a challenge to biomedical research. There are more [...] Read more.
Rare and ultra-rare diseases constitute a significant medical challenge due to their low prevalence and the limited understanding of their origin and underlying mechanisms. These disorders often exhibit phenotypic diversity and molecular complexity that represent a challenge to biomedical research. There are more than 6000 different rare diseases that affect nearly 300 million people worldwide. However, the prevalence of each rare disease is low, and in consequence, the biomedical resources dedicated to each rare disease are limited and insufficient to effectively achieve progress in the research. The use of animal models to investigate the mechanisms underlying pathogenesis has become an invaluable tool. Among the animal models commonly used in research, Drosophila melanogaster has emerged as an efficient and reliable experimental model for investigating a wide range of genetic disorders, and to develop therapeutic strategies for rare and ultra-rare diseases. It offers several advantages as a research model including short life cycle, ease of laboratory maintenance, rapid life cycle, and fully sequenced genome that make it highly suitable for studying genetic disorders. Additionally, there is a high degree of genetic conservation from Drosophila melanogaster to humans, which allows the extrapolation of findings at the molecular and cellular levels. Here, I examine the role of Drosophila melanogaster as a model for studying rare and ultra-rare diseases and highlight its significant contributions and potential to biomedical research. High-throughput next-generation sequencing (NGS) technologies, such as whole-exome sequencing and whole-genome sequencing (WGS), are providing massive amounts of information on the genomic modifications present in rare diseases and common complex traits. The sequencing of exomes or genomes of individuals affected by rare diseases has enabled human geneticists to identify rare variants and identify potential loci associated with novel gene–disease relationships. Despite these advances, the average rare disease patient still experiences significant delay until receiving a diagnosis. Furthermore, the vast majority (95%) of patients with rare conditions lack effective treatment or a cure. This scenario is enhanced by frequent misdiagnoses leading to inadequate support. In consequence, there is an urgent need to develop model organisms to explore the molecular mechanisms underlying these diseases and to establish the genetic origin of these maladies. The aim of this review is to discuss the advantages and limitations of Drosophila melanogaster, hereafter referred as Drosophila, as an experimental model for biomedical research, and the applications to study human disease. The main question to address is whether Drosophila is a valid research model to study human disease, and in particular, rare and ultra-rare diseases. Full article
(This article belongs to the Section Role of Insects in Human Society)
Show Figures

Figure 1

Figure 1
<p>Schematic representation of the genetic alterations that cause rare diseases. Single gene mutations cause genetic aberrations in the DNA sequence of a specific gene. Copy number variation (CNV): the sequences of the genome are repeated. Mitochondrial mutations, the DNA contained in the mitochondria (mtDNA) is mutated. Chromosomal abnormality, the morphology or the number of chromosomes is altered. Polygenic inheritance, more than one gene is mutated. Image generated in <a href="http://BioRender.com" target="_blank">BioRender.com</a> (10 October 2024).</p>
Full article ">Figure 2
<p>CRISPR/Cas9 system used in <span class="html-italic">Drosophila</span> to generate genetic avatars. Representation of CRISPR/Cas9 system to generate mutants in F1 generation; after the cross of parental lines expressing <span class="html-italic">Cas9</span> in the germ line, embryos are injected with a plasmid containing the required tools to induce the excision of a region of DNA (exonuclease activity of Cas9), and the re-insertion of the mutated form of the same piece of DNA (endonuclease activity of Cas9). The resulting combination produces the substitution of endogenous exons by mutated exons that reproduce the mutations found in patients. In addition, the plasmid carries a GFP under the control of a constitutive promoter (Actin) to identify the flies that undergo CRIPSR/Cas9 substitution. This GFP cDNA is flanked by two FRT siter to be excised if required in an additional cross with flies that express <span class="html-italic">flipase</span>. Image generated in <a href="http://BioRender.com" target="_blank">BioRender.com</a> (accessed on 10 October 2024).</p>
Full article ">
22 pages, 3176 KiB  
Article
Using Multi-Sensor Data Fusion Techniques and Machine Learning Algorithms for Improving UAV-Based Yield Prediction of Oilseed Rape
by Hongyan Zhu, Shikai Liang, Chengzhi Lin, Yong He and Jun-Li Xu
Drones 2024, 8(11), 642; https://doi.org/10.3390/drones8110642 - 5 Nov 2024
Viewed by 561
Abstract
Accurate and timely prediction of oilseed rape yield is crucial in precision agriculture and field remote sensing. We explored the feasibility and potential for predicting oilseed rape yield through the utilization of a UAV-based platform equipped with RGB and multispectral cameras. Genetic algorithm–partial [...] Read more.
Accurate and timely prediction of oilseed rape yield is crucial in precision agriculture and field remote sensing. We explored the feasibility and potential for predicting oilseed rape yield through the utilization of a UAV-based platform equipped with RGB and multispectral cameras. Genetic algorithm–partial least square was employed and evaluated for effective wavelength (EW) or vegetation index (VI) selection. Additionally, different machine learning algorithms, i.e., multiple linear regression (MLR), partial least squares regression (PLSR), least squares support vector machine (LS-SVM), back propagation neural network (BPNN), extreme learning machine (ELM), and radial basis function neural network (RBFNN), were developed and compared. With multi-source data fusion by combining vegetation indices (color and narrow-band VIs), robust prediction models of yield in oilseed rape were built. The performance of prediction models using the combination of VIs (RBFNN: Rpre = 0.8143, RMSEP = 171.9 kg/hm2) from multiple sensors manifested better results than those using only narrow-band VIs (BPNN: Rpre = 0.7655, RMSEP = 188.3 kg/hm2) from a multispectral camera. The best models for yield prediction were found by applying BPNN (Rpre = 0.8114, RMSEP = 172.6 kg/hm2) built from optimal EWs and ELM (Rpre = 0.8118, RMSEP = 170.9 kg/hm2) using optimal VIs. Taken together, the findings conclusively illustrate the potential of UAV-based RGB and multispectral images for the timely and non-invasive prediction of oilseed rape yield. This study also highlights that a lightweight UAV equipped with dual-image-frame snapshot cameras holds promise as a valuable tool for high-throughput plant phenotyping and advanced breeding programs within the realm of precision agriculture. Full article
Show Figures

Figure 1

Figure 1
<p>The general location of the experimental site and an overview of the image obtained by an unmanned aerial vehicle (UAV) remote sensing platform for the oilseed rape field at Zhejiang University.</p>
Full article ">Figure 2
<p>Research technology flowchart.</p>
Full article ">Figure 3
<p>The UAV platform and integrated sensors: (<b>a</b>) POS system; (<b>b</b>) RGB and multispectral imaging sensors; (<b>c</b>) ground control system; (<b>d</b>) flight of the UAV system.</p>
Full article ">Figure 4
<p>The stitching grayscale image of the multispectral image at 889 nm by Agisoft Photoscan.</p>
Full article ">Figure 5
<p>The average reflectance spectra and standard deviation (SD) of the oilseed rape canopy under different nitrogen treatments of ‘0 kg/hm<sup>2</sup>’, ‘75 kg/hm<sup>2</sup>’, ‘150 kg/hm<sup>2</sup>’, and ‘225 kg/hm<sup>2</sup>’ based on the UAV multispectral imaging system.</p>
Full article ">Figure 6
<p>Contour maps of the correlation coefficients between yield and narrow-band vegetation indices using the thorough combinations of two wavebands at λ1 and λ2 nm based on the UAV multispectral imaging system.</p>
Full article ">Figure 6 Cont.
<p>Contour maps of the correlation coefficients between yield and narrow-band vegetation indices using the thorough combinations of two wavebands at λ1 and λ2 nm based on the UAV multispectral imaging system.</p>
Full article ">Figure 7
<p>Selection of effective features by genetic algorithm–partial least squares according to smoothed frequency of selections: (<b>a</b>) wavelengths; (<b>b</b>) vegetation indices.</p>
Full article ">Figure 8
<p>The scatter plot of predicted yield values vs. measured yield values: the optimal models with full wavelengths (<b>a</b>) and selected EWs (<b>c</b>) from the multispectral camera; the optimal models with all the vegetation indices (<b>b</b>) and selected VIs (<b>d</b>) from the dual cameras.</p>
Full article ">
27 pages, 5800 KiB  
Article
Multimodal Deep Learning Integration of Image, Weather, and Phenotypic Data Under Temporal Effects for Early Prediction of Maize Yield
by Danial Shamsuddin, Monica F. Danilevicz, Hawlader A. Al-Mamun, Mohammed Bennamoun and David Edwards
Remote Sens. 2024, 16(21), 4043; https://doi.org/10.3390/rs16214043 - 30 Oct 2024
Viewed by 750
Abstract
Maize (Zea mays L.) has been shown to be sensitive to temperature deviations, influencing its yield potential. The development of new maize hybrids resilient to unfavourable weather is a desirable aim for crop breeders. In this paper, we showcase the development of [...] Read more.
Maize (Zea mays L.) has been shown to be sensitive to temperature deviations, influencing its yield potential. The development of new maize hybrids resilient to unfavourable weather is a desirable aim for crop breeders. In this paper, we showcase the development of a multimodal deep learning model using RGB images, phenotypic, and weather data under temporal effects to predict the yield potential of maize before or during anthesis and silking stages. The main objective of this study was to assess if the inclusion of historical weather data, maize growth captured through imagery, and important phenotypic traits would improve the predictive power of an established multimodal deep learning model. Evaluation of the model performance when training from scratch showed its ability to accurately predict ~89% of hybrids with high-yield potential and demonstrated enhanced explanatory power compared with previously published models. Shapley Additive explanations (SHAP) analysis indicated the top influential features include plant density, hybrid placement in the field, date to anthesis, parental line, temperature, humidity, and solar radiation. Including weather historical data was important for model performance, significantly enhancing the predictive and explanatory power of the model. For future research, the use of the model can move beyond maize yield prediction by fine-tuning the model on other crop data, serving as a potential decision-making tool for crop breeders to determine high-performing individuals from diverse crop types. Full article
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Schematic overview of the multimodal deep learning model development.</p>
Full article ">Figure 2
<p>Overview schematic of the image preprocessing workflow: (<b>A</b>) Orthomosaic generation; Metashape (v2.0.1, Agisoft) was used to combine images taken from UAV platforms. (<b>B</b>) Visualisation of shapefiles on orthomosaics using QGIS (v3.28.7-Firenze, Free Software Foundation); shapefiles were generated using the plotshpcreate package on R software (v2022.12.0 + 353, Posit Software) and imported alongside their respective orthomosaic for each year in QGIS. (<b>C</b>) Cropped images generated from QGIS; individual plot images were cropped from aligned shapefiles using QGIS (v3.28.7-Firenze, Free Software Foundation) and were further rotated horizontally and split into four quadrants using a custom Python script. A1–2 and B1–2 represent replicates of the same hybrid lineage. (<b>D</b>) Final processed image; the quadrants corresponding to one replicate stacked on top of each other for all timepoints. Processed images will serve as inputs to the image-specific model, where the pixel values will used for feature extraction and model learning.</p>
Full article ">Figure 3
<p>Overview of the layers used within each modality. All modules generate a target prediction. Repeated layers are visualised by a dotted line: (<b>A</b>) tab-DNN architecture with embedding layers for each categorical feature. (<b>B</b>) XResNet18 Image-DNN architecture; this framework uses a modified ResNet architecture with 18 layers, including a self-attention layer. (<b>C</b>) DenseNet-121 Image-DNN architecture; This framework uses four dense blocks, each with 6, 12, 24, 16 dense layers, respectively. Each dense block is connected by a transition layer.</p>
Full article ">Figure 4
<p>Overview of the fusion module used within the multimodal architecture. The last layer weights from each module are concatenated, then undergo linear and non-linear transformations before generating a yield prediction.</p>
Full article ">Figure 5
<p>Maize hybrid and yield distribution throughout the years visualised using a boxplot portraying ground truth yields by year (<b>A</b>), boxplot portraying ground truth yields by environment and year (<b>B</b>), and a Venn diagram portraying the number of maize hybrids grown each year.</p>
Full article ">Figure 6
<p>Tab-DNN model holdout results: (<b>A</b>) Comparison between predictions (<span class="html-italic">y</span>-axis) vs. ground truth yields (<span class="html-italic">x</span>-axis) for each sample per year, measured at tonne/ha. (<b>B</b>) Comparison between prediction error (<span class="html-italic">y</span>-axis) vs. environment (<span class="html-italic">x</span>-axis) for each sample per year. Calculated by subtracting the predicted yield from ground truth. (<b>C</b>) Confusion matrix showing the percentage of correctly predicted hybrids that were classified as either high (&gt;10 tonne/ha), medium (6–9.99 tonne/ha), or low (&lt;6 tonne/ha). Ground truth classifications are represented on the <span class="html-italic">y</span>-axis, with predicted performance on the <span class="html-italic">x</span>-axis.</p>
Full article ">Figure 7
<p>The top ten features based on SHAP values: (<b>A</b>) represents the top ten features and their average impact on tab-DNN model holdout predictions by mean absolute SHAP value. (<b>B</b>) represents the top ten features and their directional impact on tab-DNN model holdout predictions by SHAP value per sample. Features are further described in <a href="#app1-remotesensing-16-04043" class="html-app">Supplementary Table S2</a>.</p>
Full article ">Figure 8
<p>Image-DNN model holdout results: (<b>A</b>–<b>C</b>) Comparison between predictions (<span class="html-italic">y</span>-axis) vs. ground truth yields (<span class="html-italic">x</span>-axis) for each sample by environment and year, measured as tonne/ha. (<b>D</b>) Comparison between predictions (<span class="html-italic">y</span>-axis) vs. ground-true yields (<span class="html-italic">x</span>-axis) for each sample per year, measured as tonne/ha. (<b>E</b>) Comparison between hybrid count and yield distribution between ground truth and prediction values. (<b>F</b>) Comparison between prediction error (<span class="html-italic">y</span>-axis) vs. environment (<span class="html-italic">x</span>-axis) for each sample per year. Calculated by subtracting the predicted yield from ground truth. (<b>G</b>) Confusion matrix showing the percentage of correctly predicted hybrids that were classified as either high (&gt;10 tonne/ha), medium (6–9.99 tonne/ha), or low (&lt;6 tonne/ha). Ground truth classifications are represented on the <span class="html-italic">y</span>-axis, with predicted performance on the <span class="html-italic">x</span>-axis.</p>
Full article ">Figure 9
<p>Yield predictions were generated from the multimodal architecture when trained from scratch. Comparison between predictions (<span class="html-italic">y</span>-axis) vs. ground truth yields (<span class="html-italic">x</span>-axis) for each sample by environment and year, measured at tonne/ha using tab-DNN predictions (<b>A</b>), Image-DNN predictions (<b>B</b>), fusion module predictions (<b>C</b>), and weighted predictions (<b>D</b>). (<b>E</b>) Confusion matrix showing the percentage of correctly predicted hybrids that were classified as either high (&gt;10 tonne/ha), medium (6–9.99 tonne/ha), or low (&lt;6 tonne/ha). Ground truth classifications are represented on the <span class="html-italic">y</span>-axis, with predicted performance on the <span class="html-italic">x</span>-axis. (<b>F</b>) Classification error plot showcasing hybrids predicted as medium yield but with a ground truth class of high yield (light purple), and hybrids predicted as medium yield but with a ground truth class of low yield (mustard).</p>
Full article ">Figure 10
<p>Yield predictions generated from the multimodal architecture when using pretrained modules. Comparison between predictions (<span class="html-italic">y</span>-axis) vs. ground truth yields (<span class="html-italic">x</span>-axis) for each sample by environment and year, measured at tonne/ha using tab-DNN predictions (<b>A</b>), Image-DNN predictions (<b>B</b>), fusion module predictions (<b>C</b>), and weighted predictions (<b>D</b>). (<b>E</b>) Confusion matrix showing the percentage of correctly predicted hybrids that were classified as either high (&gt;10 tonne/ha), medium (6–9.99 tonne/ha), or low (&lt;6 tonne/ha). Ground truth classifications represented on the <span class="html-italic">y</span>-axis, with predicted performance on the <span class="html-italic">x</span>-axis. (<b>F</b>) Classification error plot showcasing hybrids predicted as medium yield but with a ground truth class of high yield (light purple), and hybrids predicted as medium yield but with a ground truth class of low yield (mustard).</p>
Full article ">
14 pages, 4130 KiB  
Article
Fermentation Profile, Bacterial Community Structure, Co-Occurrence Networks, and Their Predicted Functionality and Pathogenic Risk in High-Moisture Italian Ryegrass Silage
by Siran Wang, Chenglong Ding, Jipeng Tian, Yunhui Cheng, Nengxiang Xu, Wenjie Zhang, Xin Wang, Mudasir Nazar and Beiyi Liu
Agriculture 2024, 14(11), 1921; https://doi.org/10.3390/agriculture14111921 - 29 Oct 2024
Viewed by 449
Abstract
This study aimed to assess the fermentation characteristics, bacterial community structure, co-occurrence networks, and their predicted functionality and pathogenic risk in high-moisture Italian ryegrass (IR; Lolium multiflorum Lam.) silage. The IR harvested at heading stage (208 g dry matter (DM)/kg fresh weight) was [...] Read more.
This study aimed to assess the fermentation characteristics, bacterial community structure, co-occurrence networks, and their predicted functionality and pathogenic risk in high-moisture Italian ryegrass (IR; Lolium multiflorum Lam.) silage. The IR harvested at heading stage (208 g dry matter (DM)/kg fresh weight) was spontaneously ensiled in plastic silos (10 L scale). Triplicated silos were opened after 1, 3, 7, 15, 30, and 60 days of fermentation, respectively. The bacterial community structure on days 3 and 60 were investigated using high-throughput sequencing technology, and 16S rRNA-gene predicted functionality and phenotypes were determined by PICRUSt2 and BugBase tools, respectively. After 60 days, the IR silage exhibited good ensiling characteristics indicated by large amounts of acetic acid (~58.7 g/kg DM) and lactic acid (~91.5 g/kg DM), relatively low pH (~4.20), acceptable levels of ammonia nitrogen (~87.0 g/kg total nitrogen), and trace amounts of butyric acid (~1.59 g/kg DM). Psychrobacter was prevalent in fresh IR, and Lactobacillus became the most predominant genus after 3 and 60 days. The ensilage process reduced the complexity of the bacterial community networks in IR silage. The bacterial functional pathways in fresh and ensilaged IR are primarily characterized by the metabolism of carbohydrate and amino acid. The pyruvate kinase and 1-phosphofructokinase were critical in promoting lactic acid fermentation. A greater (p < 0.01) abundance of the “potentially pathogenic” label was noticed in the bacterial communities of ensiled IR than fresh IR. Altogether, the findings indicated that the high-moisture IR silage exhibited good ensiling characteristics, but the potential for microbial contamination and pathogens still remained after ensiling. Full article
(This article belongs to the Special Issue Silage Preparation, Processing and Efficient Utilization)
Show Figures

Figure 1

Figure 1
<p>Changes in mono-and di-saccharides in fresh and ensiled Italian ryegrass. DM: dry matter; d: days.</p>
Full article ">Figure 2
<p>Bacterial community diversity of fresh and ensiled Italian ryegrass. (<b>A</b>) Alpha diversity (Ace, Chao1, Shannon, Coverage, Simpson, and Sobs indices) of bacterial community. (<b>B</b>) Rarefaction curves based on the Sobs index. (<b>C</b>) Beta diversity of bacterial community, calculated by principal coordinate analysis (PCoA) based on the Bray−Curtis distance metric. IRFM: fresh Italian ryegrass; IR_3: spontaneously ensiled Italian ryegrass after 3 days; IR_60: spontaneously ensiled Italian ryegrass after 60 days. * 0.01 &lt; <span class="html-italic">p</span> ≤ 0.05; ** 0.001 &lt; <span class="html-italic">p</span> ≤ 0.01.</p>
Full article ">Figure 3
<p>Bacterial community compositions of fresh and ensiled Italian ryegrass. (<b>A</b>) The Venn diagram of the core OTUs in fresh and ensiled Italian ryegrass. (<b>B</b>) The bar plots of bacterial community abundances at the phylum level. (<b>C</b>) The bar plots of bacterial community abundances at the genus level. IRFM: fresh Italian ryegrass; IR_3: spontaneously ensiled Italian ryegrass after 3 days; IR_60: spontaneously ensiled Italian ryegrass after 60 days.</p>
Full article ">Figure 4
<p>Bacterial co-occurrence networks of fresh and ensiled Italian ryegrass. Bacterial co-occurrence networks (Spearman correlation, the most abundant 300 species, <span class="html-italic">p</span> value &lt; 0.05, and correlation &gt; 0.5) of IRFM (<b>A</b>), IR_3 (<b>B</b>), and IR_60 (<b>C</b>). The node represents bacterial species, node color represents bacterial phylum, and node size represents the bacterial abundance. Edges are colored according to negative (green) and positive (red) correlations. (<b>D</b>) Bar plots of node and edge numbers, respectively. (<b>E</b>) Bar plots of negative correlation proportion and correlation number. IRFM: fresh Italian ryegrass; IR_3: spontaneously ensiled Italian ryegrass after 3 days; IR_60: spontaneously ensiled Italian ryegrass after 60 days.</p>
Full article ">Figure 5
<p>The 16S ribosomal RNA gene estimated Kyoto encyclopedia of genes and genomes functional profiles at the second pathway level. IRFM: fresh Italian ryegrass; IR_3: spontaneously ensiled Italian ryegrass after 3 days; IR_60: spontaneously ensiled Italian ryegrass after 60 days.</p>
Full article ">Figure 6
<p>Changes of key enzymes involved in bacterial community metabolism during the ensiling of Italian ryegrass. (<b>A</b>) hexokinase; (<b>B</b>) 1-phosphofructokinase; (<b>C</b>) pyruvate kinase; (<b>D</b>) D-lactate dehydrogenase; (<b>E</b>) L-lactate dehydrogenase; (<b>F</b>) arginine deiminase. EC: reference metabolic pathway highlighting numbers; IRFM: fresh Italian ryegrass; IR_3: spontaneously ensiled Italian ryegrass after 3 days; IR_60: spontaneously ensiled Italian ryegrass after 60 days. Means with different letters differ significantly from each other (<span class="html-italic">p</span> &lt; 0.05).</p>
Full article ">Figure 7
<p>Bacterial phenotypes annotation of fresh and ensiled Italian ryegrass. (<b>A</b>) Phenotypes that reflect bacterial characteristics (Aerobic, Facultatively Anaerobic, Anaerobic, Gram_Negative, Gram_Positive, and Contains_Mobile_Elements) and bacterial resistance (Forms_Biofilms, Potentially_Pathogenic, and Stress_Tolerant). (<b>B</b>) Comparison of the Potentially_Pathogenic category between fresh and ensiled Italian ryegrass. The statistical analysis was conducted using a two-sided <span class="html-italic">t</span>-test. * 0.01 &lt; <span class="html-italic">p</span> ≤ 0.05; ** 0.001 &lt; <span class="html-italic">p</span> ≤ 0.01.</p>
Full article ">
30 pages, 14626 KiB  
Article
Integration of IoT Technologies and High-Performance Phenotyping for Climate Control in Greenhouses and Mitigation of Water Deficit: A Study of High-Andean Oat
by Edwin Villagran, Gabriela Toro-Tobón, Fabián Andrés Velázquez and German A. Estrada-Bonilla
AgriEngineering 2024, 6(4), 4011-4040; https://doi.org/10.3390/agriengineering6040227 - 29 Oct 2024
Viewed by 878
Abstract
Climate change has intensified droughts, severely impacting crops like oats and highlighting the need for effective adaptation strategies. In this context, the implementation of IoT-based climate control systems in greenhouses emerges as a promising solution for optimizing microclimates. These systems allow for the [...] Read more.
Climate change has intensified droughts, severely impacting crops like oats and highlighting the need for effective adaptation strategies. In this context, the implementation of IoT-based climate control systems in greenhouses emerges as a promising solution for optimizing microclimates. These systems allow for the precise monitoring and adjustment of critical variables such as temperature, humidity, vapor pressure deficit (VPD), and photosynthetically active radiation (PAR), ensuring optimal conditions for crop growth. During the experiment, the average daytime temperature was 22.6 °C and the nighttime temperature was 15.7 °C. The average relative humidity was 60%, with a VPD of 0.46 kPa during the day and 1.26 kPa at night, while the PAR reached an average of 267 μmol m−2 s−1. Additionally, the use of high-throughput gravimetric phenotyping platforms enabled precise data collection on the plant–soil–atmosphere relationship, providing exhaustive control over water balance and irrigation. This facilitated the evaluation of the physiological response of plants to abiotic stress. Inoculation with microbial consortia (PGPB) was used as a tool to mitigate water stress. In this 69-day study, irrigation was suspended in specific treatments to simulate drought, and it was observed that inoculated plants maintained chlorophyll b and carotenoid levels akin to those of irrigated plants, indicating greater tolerance to water deficit. These plants also exhibited greater efficiency in dissipating light energy and rapid recovery after rehydration. The results underscore the potential of combining IoT monitoring technologies, advanced phenotyping platforms, and microbial consortia to enhance crop resilience to climate change. Full article
(This article belongs to the Section Sensors Technology and Precision Agriculture)
Show Figures

Figure 1

Figure 1
<p>Geometric scheme and some characteristics of the experimental greenhouse.</p>
Full article ">Figure 2
<p>Description of the operation of the remote climate control system.</p>
Full article ">Figure 3
<p>Active climate control strategies installed in the experimental greenhouse. (<b>A</b>) automated natural ventilation; (<b>B</b>) evaporative cooling equipment; and (<b>C</b>) fog cooling.</p>
Full article ">Figure 4
<p>Schematic of the high-throughput phenotyping and big data system.</p>
Full article ">Figure 5
<p>Setup of the experimental trial on the phenotyping platform.</p>
Full article ">Figure 6
<p>The distribution of the experimental trial in the 32 pots of the physiological phenotyping platform.</p>
Full article ">Figure 7
<p>Temporal distribution of outdoor temperature data and hourly mean values during the experimental period.</p>
Full article ">Figure 8
<p>The temporal distribution of outdoor relative humidity data and hourly mean values during the experimental period.</p>
Full article ">Figure 9
<p>The temporal distribution of solar radiation data and hourly mean values during the experimental period.</p>
Full article ">Figure 10
<p>The temporal distribution of wind speed and hourly mean values during the experimental period.</p>
Full article ">Figure 11
<p>The temporal distribution of temperature and mean hourly values during the experimental period inside the greenhouse.</p>
Full article ">Figure 12
<p>The temporal distribution of relative humidity and mean hourly values during the experimental period inside the greenhouse.</p>
Full article ">Figure 13
<p>The temporal distribution of the vapor pressure deficit (VPD) and hourly average values during the experimental period inside the greenhouse.</p>
Full article ">Figure 14
<p>The temporal distribution of PAR and mean hourly values during the experimental period inside the greenhouse.</p>
Full article ">Figure 15
<p>Hourly average values of the thermal differential between the inside and outside environment of the greenhouse.</p>
Full article ">Figure 16
<p>The water balance parameters of oat plants with and without inoculation with a microbial consortium composed of <span class="html-italic">A. brasilense, Herbaspirillum</span> sp. and <span class="html-italic">R. leguminosarum</span> (PGPB), after 17 days of irrigation suspension (Stress), with their respective irrigated controls (Irrigation). (<b>A</b>) The volumetric water content and (<b>B</b>) the daily plant transpiration.</p>
Full article ">Figure 17
<p>The stomatal conductance of high-Andean oat plants without and with PGPB inoculation under moderate (Ms) and severe (Ss) water deficit stress and five days after hydration (H), with their respective irrigated controls. Data were analyzed using Tukey’s test for parametric data and the Kruskal–Wallis test for non-parametric data. Different letters indicate significant differences between treatments (α = 0.05).</p>
Full article ">Figure 18
<p>The relative water content (RWC) of high-Andean oat plants without and with PGPB inoculation (<b>A</b>) under severe water deficit stress and (<b>B</b>) after 5 days following hydration, with their respective irrigated controls. Asterisks indicate significant differences between treatments after a two-sided Student’s <span class="html-italic">t</span>-test. **** <span class="html-italic">p</span> ≤ 0.0001.</p>
Full article ">Figure 19
<p>Dry biomass produced by oat plants with and without inoculation with PGPB after 17 days of irrigation suspension, with their respective irrigated controls. (<b>A</b>) Shoot dry biomass, (<b>B</b>) root dry biomass, and (<b>C</b>) total dry biomass. Asterisks indicate significant differences between treatments after a two-sided Student’s <span class="html-italic">t</span>-test. * <span class="html-italic">p</span> ≤ 0.05.</p>
Full article ">Figure 20
<p>The photochemical and non-photochemical quenching parameters of chlorophyll a fluorescence of high-Andean oat plants with and without PGPB inoculation, under moderate stress (Ms), severe stress (Ss), and after five days of hydration (H) water deficit stress, with their irrigated controls. (<b>A</b>) The relative electron transfer rate (ETR); (<b>B</b>) the photochemical quantum yield of PSII—Y (II); (<b>C</b>) the quantum yield of unregulated heat dissipation—Y (NO); (<b>D</b>) the non-photochemical quantum yield—Y (NPQ). Data were analyzed using Tukey’s test for parametric data and the Kruskal–Wallis test for non-parametric data. Different letters indicate significant differences between treatments (α = 0.05).</p>
Full article ">Figure 21
<p>The content of (<b>A</b>) chlorophyll a, (<b>B</b>) chlorophyll b, and (<b>C</b>) carotenoids in high-Andean oat plants inoculated and non-inoculated with the species <span class="html-italic">A. brasilense</span>, <span class="html-italic">Herbaspirillum</span> sp. and <span class="html-italic">R. leguminosarum</span> (PGPB) under moderate water deficit stress (Ms), severe water deficit stress (Ss), and five days after hydration (H). Compared with their respective irrigated controls. Data were analyzed using Tukey’s test for parametric data and the Kruskal–Wallis test for non-parametric data. Different letters indicate significant differences between treatments (α = 0.05).</p>
Full article ">
17 pages, 11054 KiB  
Article
Advanced Plant Phenotyping: Unmanned Aerial Vehicle Remote Sensing and CimageA Software Technology for Precision Crop Growth Monitoring
by Hongyu Fu, Jianning Lu, Guoxian Cui, Jihao Nie, Wei Wang, Wei She and Jinwei Li
Agronomy 2024, 14(11), 2534; https://doi.org/10.3390/agronomy14112534 - 28 Oct 2024
Viewed by 652
Abstract
In production activities and breeding programs, large-scale investigations of crop high-throughput phenotype information are needed to help improve management and decision-making. The development of UAV (unmanned aerial vehicle) remote sensing technology provides a new means for the large-scale, efficient, and accurate acquisition of [...] Read more.
In production activities and breeding programs, large-scale investigations of crop high-throughput phenotype information are needed to help improve management and decision-making. The development of UAV (unmanned aerial vehicle) remote sensing technology provides a new means for the large-scale, efficient, and accurate acquisition of crop phenotypes, but its practical application and popularization are hindered due to the complicated data processing required. To date, there is no automated system that can utilize the canopy images acquired through UAV to conduct a phenotypic character analysis. To address this bottleneck, we developed a new scalable software called CimageA. CimageA uses crop canopy images obtained by UAV as materials. It can combine machine vision technology and machine learning technology to conduct the high-throughput processing and phenotyping of crop remote sensing data. First, zoning tools are applied to draw an area-of-interest (AOI). Then, CimageA can rapidly extract vital remote sensing information such as the color, texture, and spectrum of the crop canopy in the plots. In addition, we developed data analysis modules that estimate and quantify related phenotypes (such as leaf area index, canopy coverage, and plant height) by analyzing the association between measured crop phenotypes and CimageA-derived remote sensing eigenvalues. Through a series of experiments, we confirmed that CimageA performs well in extracting high-throughput remote sensing information regarding crops, and verified the reliability of retrieving LAI (R2 = 0.796) and estimating plant height (R2 = 0.989) and planting area using CimageA. In short, CimageA is an efficient and non-destructive tool for crop phenotype analysis, which is of great value for monitoring crop growth and guiding breeding decisions. Full article
(This article belongs to the Section Precision and Digital Agriculture)
Show Figures

Figure 1

Figure 1
<p>Plots’ distribution and multi-channel remote sensing images of ramie test area. Result_RGB is the highest-resolution visible image of the test area, Result_Blue shows the blue channel spectral reflectance image, Result_Green shows the green channel spectral reflectance image, Result_Red shows the red channel spectral reflectance image, Result_RedEdge shows the red edge channel spectral reflectance image and Result_Nir is the spectral reflectance image of the near infrared channel.</p>
Full article ">Figure 2
<p>User interface and instructions of CimageA. Numbers in the figure provide guidance: (1) File processing. (2) Image processing and parameter settings. (3) Data processing. (4) Data analysis. (5) Phenotypic visualization.</p>
Full article ">Figure 3
<p>The operation process of CimageA.</p>
Full article ">Figure 4
<p>Tilt correction of CimageA. (<b>A</b>) Image before tilt correction. (<b>B</b>) Image after tilt correction.</p>
Full article ">Figure 5
<p>Intelligent drawing of the AOI. (<b>A</b>) Plots divided by diagonal zoning. (<b>B</b>) Plots divided by hypotenuse zoning.</p>
Full article ">Figure 6
<p>The result of ramie plant segmentation based on ExR. (<b>A</b>) Original image. (<b>B</b>) Grayscale image in ExR channel. (<b>C</b>) The segmented image.</p>
Full article ">Figure 7
<p>Three ramie varieties with different leaf color. (<b>A</b>) Dazhu ramie. (<b>B</b>) Xiangzhu 7. (<b>C</b>) Changshun ramie.</p>
Full article ">Figure 8
<p>Extraction of ramie plant height based on UAV remote sensing images.</p>
Full article ">Figure 9
<p>Crop planting area measurement. Area 1 (red area) is the ramie planting area, area 2 (green area) is the jute planting area, area 3 (blue area) is the cabbage planting area, and area 4 includes all planting areas.</p>
Full article ">Figure 10
<p>Quantitative color features of the ramie with different leaf color classes.</p>
Full article ">Figure 11
<p>Quantitative leaf color evaluation of three ramie varieties.</p>
Full article ">Figure 12
<p>Correlation analysis between ramie LAI and remote sensing eigenvalues.</p>
Full article ">Figure 13
<p>Relationship between the measured plant height and the estimated plant height.</p>
Full article ">Figure 14
<p>Temporal changes of the estimated plant height.</p>
Full article ">Figure 15
<p>Precision of crop area extracted by CimageA. ** indicates a significant correlation.</p>
Full article ">
20 pages, 4450 KiB  
Article
Comprehensive Analysis of Phenotypic Traits in Chinese Cabbage Using 3D Point Cloud Technology
by Chongchong Yang, Lei Sun, Jun Zhang, Xiaofei Fan, Dongfang Zhang, Tianyi Ren, Minggeng Liu, Zhiming Zhang and Wei Ma
Agronomy 2024, 14(11), 2506; https://doi.org/10.3390/agronomy14112506 - 25 Oct 2024
Viewed by 533
Abstract
Studies on the phenotypic traits and their associations in Chinese cabbage lack precise and objective digital evaluation metrics. Traditional assessment methods often rely on subjective evaluations and experience, compromising accuracy and reliability. This study develops an innovative, comprehensive trait evaluation method based on [...] Read more.
Studies on the phenotypic traits and their associations in Chinese cabbage lack precise and objective digital evaluation metrics. Traditional assessment methods often rely on subjective evaluations and experience, compromising accuracy and reliability. This study develops an innovative, comprehensive trait evaluation method based on 3D point cloud technology, with the aim of enhancing the precision, reliability, and standardization of the comprehensive phenotypic traits of Chinese cabbage. By using multi-view image sequences and structure-from-motion algorithms, 3D point clouds of 50 plants from each of the 17 Chinese cabbage varieties were reconstructed. Color-based region growing and 3D convex hull techniques were employed to measure 30 agronomic traits. Comparisons between 3D point cloud-based measurements of the plant spread, plant height, leaf area, and leaf ball volume and traditional methods yielded R2 values greater than 0.97, with root mean square errors of 1.27 cm, 1.16 cm, 839.77 cm3, and 59.15 cm2, respectively. Based on the plant spread and plant height, a linear regression prediction of Chinese cabbage weights was conducted, yielding an R2 value of 0.76. Integrated optimization algorithms were used to test the parameters, reducing the measurement time from 55 min when using traditional methods to 3.2 min. Furthermore, in-depth analyses including variation, correlation, principal component analysis, and clustering analyses were conducted. Variation analysis revealed significant trait variability, with correlation analysis indicating 21 pairs of traits with highly significant positive correlations and 2 pairs with highly significant negative correlations. The top six principal components accounted for 90% of the total variance. Using the elbow method, k-means clustering determined that the optimal number of clusters was four, thus classifying the 17 cabbage varieties into four distinct groups. This study provides new theoretical and methodological insights for exploring phenotypic trait associations in Chinese cabbage and facilitates the breeding and identification of high-quality varieties. Compared with traditional methods, this system provides significant advantages in terms of accuracy, speed, and comprehensiveness, with its low cost and ease of use making it an ideal replacement for manual methods, being particularly suited for large-scale monitoring and high-throughput phenotyping. Full article
Show Figures

Figure 1

Figure 1
<p>Sequence photo collection process of Chinese cabbage. (<b>A</b>) Planting site at Sanfen Farm at Hebei Agricultural University, where 17 cabbage varieties were grown using a randomized block design. (<b>B</b>) Image acquisition equipment used to capture images of the Chinese cabbage, ensuring high-quality data collection. (<b>C</b>) Image capture process, with the camera rotating around each cabbage sample on a transparent platform, collecting 60–70 images from top and bottom angles for comprehensive phenotypic analysis.</p>
Full article ">Figure 2
<p>Flowchart of 3D point cloud reconstruction of Chinese cabbage. (<b>A</b>) Image acquisition of Chinese cabbage from multiple angles. (<b>B</b>) Initial 3D point cloud reconstruction using structure from motion (SfM) algorithms and Agisoft Metashape 1.7.0. (<b>C</b>) Preprocessed point cloud models showing top and bottom views after noise removal and background segmentation using pass-through filtering, color threshold-based segmentation, and statistical filtering algorithms. (<b>D</b>) Registered and aligned point cloud models from the top and bottom views achieved through coarse registration with a reference cube and fine-tuning using the iterative closest point (ICP) algorithm.</p>
Full article ">Figure 3
<p>Measurement flow chart. (<b>A</b>) The formula for the leaf shape index, calculated using the leaf length and width. (<b>B</b>) The selection of positions for the leaf base and tip. (<b>C</b>) The convex hull of the cabbage leaf point cloud. (<b>D</b>) The segmentation of the midrib into three parts for width and thickness analysis. (<b>E</b>) The bounding box surrounding the cabbage leaf point cloud. (<b>F</b>) The point cloud of the cabbage leaf. (<b>G</b>) The segmentation of the leaf using a color-based region-growing algorithm. (<b>H</b>) The point cloud of the midrib in the cabbage leaf. (<b>I</b>) The convex hull of the midrib point cloud. (<b>J</b>) The manual counting of the total leaf and bulb leaf numbers. (<b>K</b>) The weight measurement of the cabbage using a high-precision electronic scale. (<b>L</b>) A 3D point cloud model of the cabbage plant. (<b>M</b>) The convex hull of the cabbage plant point cloud. (<b>N</b>) The bounding box of the midrib point cloud. (<b>O</b>) The formula for calculating the compactness of the leaf ball. (<b>P</b>) The bounding box of the cabbage plant point cloud. (<b>Q</b>) The point cloud model of the central pillar of the cabbage. (<b>R</b>) The bounding box of the central pillar point cloud. (<b>S</b>) The convex hull calculation for the central pillar. (<b>T</b>) The measurements of the upper and lower widths of the center column, with models showing the original, PCA-aligned, and separated point clouds.</p>
Full article ">Figure 4
<p>Comparative analysis of data. Note: Comparison of 3D point cloud and conventional measurements for Chinese cabbage using regression lines and statistical metrics.</p>
Full article ">Figure 5
<p>Cabbage weight regression analysis. Note: This 3D scatter plot illustrates the regression relationship between the cabbage weight and the variables of plant height and spread. The red plane represents the model’s fit, with annotations for the regression equation, R<sup>2</sup>, and RMSE.</p>
Full article ">Figure 6
<p>Correlation analysis heatmap and hierarchical clustering dendrogram. Note: This figure shows two visualizations: a heatmap (<b>A</b>) which highlights the strength of correlations between these variables and a dendrogram (<b>B</b>) representing hierarchical clustering of variables. The heatmap uses color intensity to depict the strength and direction of these correlations, where warm colors show positive relationships and cool colors indicate negative ones. The dendrogram groups variables based on similarity, with shorter branches indicating higher correlation.</p>
Full article ">Figure 7
<p>Visualization of principal component analysis results. (<b>A</b>) A scree plot displaying the variance explained by each principal component. (<b>B</b>) A cumulative variance explained plot indicating that the first six components together accounted for 90% of the data variability.</p>
Full article ">Figure 8
<p>Cluster analysis diagram. (<b>A</b>) The elbow method was used to determine the optimal number of clusters by plotting the sum of squared distances against the number of clusters. The point of inflection at four clusters indicates the optimal k value. (<b>B</b>) Cluster analysis was conducted using the k-means method after applying principal component analysis (PCA) for dimensionality reduction. The plot shows 17 cabbage varieties divided into four clusters based on the first two principal components. Each cluster represents distinct morphological and physiological characteristics, providing insights into the diversity and potential for variety improvement.</p>
Full article ">
13 pages, 4060 KiB  
Article
Monitoring of Broccoli Flower Head Development in Fields Using Drone Imagery and Deep Learning Methods
by Chenzi Zhang, Xiaoxue Sun, Shuxin Xuan, Jun Zhang, Dongfang Zhang, Xiangyang Yuan, Xiaofei Fan and Xuesong Suo
Agronomy 2024, 14(11), 2496; https://doi.org/10.3390/agronomy14112496 - 25 Oct 2024
Viewed by 481
Abstract
For different broccoli materials, it used to be necessary to manually plant in a large area for the investigation of flower ball information, and this method is susceptible to subjective influence, which is not only time-consuming and laborious but may also cause some [...] Read more.
For different broccoli materials, it used to be necessary to manually plant in a large area for the investigation of flower ball information, and this method is susceptible to subjective influence, which is not only time-consuming and laborious but may also cause some damage to the broccoli in the process of investigation. Therefore, the rapid and nondestructive monitoring of flower heads is key to acquiring high-throughput phenotypic information on broccoli crops. In this study, we used an unmanned aerial vehicle (UAV) to acquire hundreds of images of field-grown broccoli to evaluate their flower head development rate and sizes during growth. First, YOLOv5 and YOLOv8 were used to complete the position detection and counting statistics at the seedling and heading stages. Then, UNet, PSPNet, DeepLabv3+, and SC-DeepLabv3+ were used to segment the flower heads in the images. The improved SC-DeepLabv3+ model excelled in segmenting flower heads, showing Precision, reconciled mean F1-score, mean intersection over union, and mean pixel accuracy values of 93.66%, 95.24%, 91.47%, and 97.24%, respectively, which were 0.57, 1.12, 1.16, and 1.70 percentage points higher than the respective values achieved with the DeepLabv3+ model. Flower head sizes were predicted on the basis of the pixel value of individual flower heads and ground sampling distance, yielding predictions with an R2 value of 0.67 and root-mean-squared error of 1.81 cm. Therefore, the development rate and sizes of broccoli flower heads during growth were successively estimated and calculated. Compared with the existing technology, it greatly improves work efficiency and can help to obtain timely information on crop growth in the field. Our methodology provides a convenient, fast, and reliable way for investigating field traits in broccoli breeding. Full article
Show Figures

Figure 1

Figure 1
<p>Network structure diagrams. (<b>a</b>) YOLOv5 network structure. (<b>b</b>) YOLOv8 network structure.</p>
Full article ">Figure 2
<p>Network architecture of SC-DeepLabv3+, an improved DeepLabv3+ model.</p>
Full article ">Figure 3
<p>Broccoli target detection results: (<b>a</b>) two seedlings were not successfully balled; (<b>b</b>) all of the seedlings successfully balled.</p>
Full article ">Figure 4
<p>Assessment of the segmentation analysis. (<b>a</b>) Flower head segmentation results from the different models. (<b>b</b>) Segmentation accuracy of the different models. The value inside the center of the radar plot is 90% and the value at the outer boundary is 98%. (<b>c</b>) Evaluation results of different segmentation models.</p>
Full article ">Figure 5
<p>Broccoli flower head size estimation. (<b>a</b>) Results of flower head segmentation. (<b>b</b>) Calculated flower head sizes in the field area. (<b>c</b>) Comparison of predicted flower head sizes with actual measurements.</p>
Full article ">
16 pages, 7713 KiB  
Article
Digital Magnetic Sorting for Fractionating Cell Populations with Variable Antigen Expression in Cell Therapy Process Development
by Savannah Bshara-Corson, Andrew Burwell, Timothy Tiemann and Coleman Murray
Magnetochemistry 2024, 10(11), 81; https://doi.org/10.3390/magnetochemistry10110081 - 23 Oct 2024
Viewed by 794
Abstract
Cellular therapies exhibit immense potential in treating complex diseases with sustained responses. The manufacture of cell therapies involves the purification and engineering of specific cells from a donor or patient to achieve a therapeutic response upon injection. Magnetic cell sorting targeting the presence [...] Read more.
Cellular therapies exhibit immense potential in treating complex diseases with sustained responses. The manufacture of cell therapies involves the purification and engineering of specific cells from a donor or patient to achieve a therapeutic response upon injection. Magnetic cell sorting targeting the presence or absence of surface markers is commonly used for upfront purification. However, emerging research shows that optimal therapeutic phenotypes are characterized not only by the presence or absence of specific antigens but also by antigen density. Unfortunately, current cell purification tools like magnetic or fluorescence-activated cell sorting (FACS) lack the resolution to differentiate populations based on antigen density while maintaining scalability. Utilizing a technique known as digital magnetic sorting (DMS), we demonstrate proof of concept for a scalable, magnetic-based approach to fractionate cell populations based on antigen density level. Targeting CD4 on human leukocytes, DMS demonstrated fractionation into CD4Hi T cells and CD4Low monocytes and neutrophils as quantified by flow cytometry and single-cell RNA seq. DMS also demonstrated high throughput processing at throughputs 3–10× faster than FACS. We believe DMS can be leveraged and scaled to enable antigen density-based sorting in cell therapy manufacturing, leading to the production of more potent and sustainable cellular therapies. Full article
(This article belongs to the Section Applications of Magnetism and Magnetic Materials)
Show Figures

Figure 1

Figure 1
<p>(<b>A</b>) The DMS instrument houses a closed loop controlled internal magnetic wheel of rare earth magnets arranged in a Halbach array orientation to generate a cycling magnetic field and a peristaltic pump to perform automated sample processing through the DMS cartridge. The DMS cartridge is loaded into the benchtop DMS instrument and connected to the peristaltic pump. A tablet with a simple graphic user interface walks the user through a standardized protocol. (<b>B</b>) The DMS cartridge workflow consists of: (1) Buffer priming to fill the cartridge with working buffer. (2) Sample injection, where the magnetically tagged cell sample is circulated through the base of the cartridge chamber, where magnetized cells are captured and pulled onto the micromagnetic substrate. (3) Digital magnetic sorting, where cells are transported vertically across the micromagnetic substrate and separate into High Magnetic and Low Magnetic Fractions. (4) Extraction after DMS separation. A flexible “extraction lock” valve is actuated to separate the High Magnetic Fraction from the Low Magnetic Fraction. Each fraction can be eluted from the cartridge via Luer Lock ports.</p>
Full article ">Figure 2
<p>Flow cytometry analysis of leukocytes processed through the DMS high/low workflow. (<b>A</b>) Before separation, the sample was observed to have CD4<sup>Hi</sup> and CD4<sup>Low</sup> populations, with &gt;99% of the cells being CD45+. The peak-to-peak difference between CD4<sup>Hi</sup> and CD4<sup>Low</sup> populations was observed to be approximately 0.6 log. (<b>B</b>) The HMF demonstrated enrichment of CD4<sup>Hi</sup> cells with 81.8% ± 4.1 (N = 3 trials) of cells occupying the CD4<sup>Hi</sup> gate, while (<b>C</b>) the LMF showed enrichment of CD4<sup>Low</sup> cells with 78.8% ± 4.3% (N = 3 trials). (<b>D</b>) Histogram plots of CD4 expression between HMF and LMF show significant peak expression with a peak-to-peak difference that correlates with the original sample. ANOVA analysis of the HMFs and LMFs showed statistically significant differences in CD4 expression (<span class="html-italic">p</span> value &lt; 0.0001).</p>
Full article ">Figure 3
<p>Multi-donor fractionation performance of the DMS system of CD4<sup>Hi</sup> and CD4<sup>Low</sup> populations into the High Magnetic Fraction (HMF) and Low Magnetic Fraction (LMF), respectively. Donor samples had varying ratios of CD4(−), CD4<sup>Hi</sup>, and CD4<sup>Low</sup> phenotypes typically split between 40 and 50% CD4(−), between 40 and 50% CD4<sup>Low</sup>, and from 5 to 10% CD4<sup>Hi</sup>. Phenotypic gating based on CD4 expression was carried out using flow cytometry of samples before separation and in the LMFs and HMFs. Separations across multiple donors show high enrichment of CD4<sup>Hi</sup> cells in the HMF (77.4 ± 4.9% donor average) and high enrichment of CD4<sup>Low</sup> cells in the LMF (79.9 ± 7.5% donor average). All DMS fractionations show &gt;90% purity of CD4(+) cells in both HMFs and LMFs.</p>
Full article ">Figure 4
<p>Single-cell RNA seq using a 10× GEM chip was performed with leukocyte samples before separation and after fractionation into LMF and HMF. (<b>A</b>) Before separation, leukocytes exhibited expected populations of B cells; monocytes and neutrophils; T cells and their subtypes; and other various cell subpopulations. (<b>B</b>) tSNE plot of LMF shows enrichment of monocytes (78%) and some T cells (11%). (<b>C</b>) The HMF shows enrichment of T cells (52%) and subpopulations (~14%) with some monocytes and neutrophils (31%). These data confirm that the DMS system is fractionating populations, which have been shown to correlate with CD4 high/low expression and align with published data and scientific publications.</p>
Full article ">Figure 5
<p>Single-cell RNA sequencing of LMF and DMS fractions. (<b>A</b>) Cells in HMF show significantly higher expression of multiple T cell markers compared to the LMF (<b>B</b>), which showed significantly higher expression of monocyte and neutrophil markers. These data show that the HMF preferentially enriches T cell populations, while the LMF is preferentially enriching for monocyte and neutrophil populations. Note * <span class="html-italic">p</span> &lt; 0.05, ** <span class="html-italic">p</span> &lt; 0.01, and **** <span class="html-italic">p</span> &lt; 0.0001.</p>
Full article ">Figure 6
<p>Throughput comparison between DMS, FACS, magnetic enrichment/depletion, and FACS. (<b>A</b>) Due to its parallelized nature, DMS can achieve higher throughput processing of antigen density subpopulations. DMS is centrifugation-free and requires 30 min of magnetic labeling, rapid magnetic debulking, or preselection followed by antigen density fractionation via DMS. (<b>B</b>) Contrastingly, FACS workflows require all cells in a sample to be analyzed and sorted, which leads to long sort times that negatively impact viability and throughput. A FACS workflow requires antibody labeling and multiple centrifugation and wash steps before proceeding to FACS. Assuming a 1 × 10<sup>4</sup> cell per second throughput, a sample of 1 × 10<sup>8</sup> cells would take about 2.7 h to complete. (<b>C</b>) Oftentimes magnetic enrichment or depletion is performed to debulk a sample before FACS. The first magnetic enrichment/depletion is performed by magnetic labeling and pulldown via bulk magnet or column. Then, the sample is processed for FACS isolation. While this FACS + MACS method increases throughput, it requires the addition of multiple sample processing steps, which reduces target cell yield and increases consumable and reagent costs significantly. Typically, a post-sort process evaluation step is performed by sampling cell viability and cell counts with hemocytometers or automated cell counting techniques. Commonly, post-sort cell count and viability checks are performed regardless of sorting platform (DMS, FACS, or MACS + FACS) and typically take ~5 min to complete.</p>
Full article ">Figure 7
<p>Comparison of DMS and other cell sorting footprints. MACS GMP scale systems for cell therapy production can accommodate large cell quantities but cannot fractionate subpopulations based on antigen density. FACS sorters have quantitative gating capabilities but are too low throughput to accommodate cell therapy production scale. DMS offers a scalable approach to sorting based on antigen density by cartridge parallelization. A single DMS benchtop sorter can achieve a 10<sup>8</sup> cell/h throughput in a single run. The scaled DMS GMP system has eight parallelized cartridges that operate in semicontinuous batch mode to achieve a total throughput of 10<sup>9</sup> cells/h.</p>
Full article ">Figure A1
<p>(<b>A</b>) Initial gating on lymphocytes was performed on the native sample, utilizing forward scatter (FSC) and side scatter (SSC) parameters to segregate the lymphocyte population based on their intrinsic physical properties. (<b>B</b>) The ‘Single Cells’ gate was established to ensure analysis of individual cell events, thereby excluding multi-cell aggregates and enhancing the accuracy of the subsequent gating steps. (<b>C</b>) Viability assessment was conducted using Thermo Fisher Scientific’s SYTOX™ Advanced™ Dead Cell Stain Kit (Catalog No. S10349). This step involved gating live cells to exclude non-viable cells, ensuring the analysis focused solely on live cellular constituents. (<b>D</b>) Discrimination of CD4<sup>Hi</sup> and CD4<sup>Low</sup> cell populations was then performed. CD4<sup>Hi</sup> cells, typically representative of T cell populations, were gated based on their characteristic expression patterns. In contrast, CD4<sup>Low</sup> cells, correlating with expected monocyte populations, were gated accordingly. This distinction was pivotal in reflecting the anticipated frequencies of these cell types within peripheral blood mononuclear cells (PBMCs).</p>
Full article ">Figure A2
<p>Primer on digital magnetic sorting (DMS), also known as ratcheting cytometry. (<b>A</b>) The DMS benchtop instrument accepts a DMS cartridge into a “saddle”, which is then hooked to a tubing and peristaltic pump to draw fluid through the cartridge. Internal and centered with the cartridge is a rotating magnetic wheel in radial Halbach array orientation. This rotating wheel drives the digital magnetic sorting behavior across the magnetic chip within the cartridge. (<b>B</b>) The cartridge itself consists of a buffer/sample input reservoir, a flexible extraction lock valve to sequester LMF and HMF regions in the chip, a three-way Luer Lock valve to connect to the peristaltic pump, and Luer Lock extraction ports to remove cells. (<b>C</b>) The DMS cartridge is a clamshell assembly containing the DMS magnetic chip, which is secured to the cartridge chamber via a backplate and die-cut double-sided adhesive strip. (<b>D</b>) The DMS chip itself is a 25 × 75 mm glass substrate with regions of high-density nickel/iron micro-bars with a 5 µm width, 25 mm length, and a thickness between 3 and 5 µm. The first region is the LMF and is 25 × 25 mm with a microbar pitch of 10 µm. This region is directly connected into the HMF region, which consists of multiple 2 mm long × 25 mm wide sections of increasing pitch. The first region has a 12 µm pitch and then increments by 2 µm with each 2 mm section (12 µm, 14 µm, 16 µm, …, 50 µm). (<b>E</b>) In response to a rotating magnetic field generated by the DMS benchtop instrument, the microbars magnetize in alignment to the bulk field, modifying the magnetic potential energy landscape and introducing potential wells in which superparamagnetic particles/cells migrate. As the wheel is cycled, particles/cells will follow the potential energy wells and jump from bar to bar against gravity based on size and bound iron oxide content. Using a chip consisting of microstructure arrays with gradient horizontal pitch, cells with different levels of magnetic binding will traverse the array until reaching their critical pitch, where they will collect and oscillate. Magnetized cells with increasing magnetic content will have correspondingly higher critical pitches and can therefore be separated. This results in temporally stable distributions of cells based on levels of magnetic content per cell or particles per cell. In this way, cell populations can be differentiated based on their level of magnetic binding (Denoted by vertical gradient arrow). Another analogy to conceptualize the process of digital magnetic sorting is to imagine the cells are climbing a “magnetic ladder” whose rungs are moving further apart the higher they climb. Magnetized cells will “jump” to a higher rung (or Ni/Fe microbar) with each cycle of the magnetic wheel as they chase the dynamic potential energy wells up the chip. At some point the magnetically tagged cells will have insufficient strength to jump to the next rung of the ladder. This “critical pitch” is determined by the amount of magnetic content bound to the cell, the pitch of the Ni/Fe microbars, and the characteristics of the magnetic field at that point on the chip (flux density [T] and rotational frequency of the wheel). This enables highly magnetic cells to climb the chip further than weakly magnetic cells, enabling multi-population sorting based on bound magnetic content.</p>
Full article ">Figure A3
<p>Numerical modeling of magnetic content per cell. We have already demonstrated in previous work [<a href="#B19-magnetochemistry-10-00081" class="html-bibr">19</a>] that cells transport and concentrate onto magnetic pitch ranges depending on the microelement pitch, cell-bound magnetic content, and the magnetic field frequency. (<b>A</b>) Using the numerical model developed previously, we can determine what level of iron oxide content per cell corresponds to a specific magnetic microstructure pitch under a given magnetic frequency. (<b>B</b>) The predictive model matches well with multiple empirical datasets where individual cells were fluorescently imaged through DMS chips to quantify the number of bound magnetic beads per cell with known levels of magnetic content per bead.</p>
Full article ">Figure A4
<p>T cell and monocyte distributions on the DMS chip. Two key system and cartridge design parameters needed to be considered for fractionating the high and low antigen density populations: Primarily the frequency of the magnetic field that the magnetic cartridge was subjected to and the micromagnetic pitch range that the cells magnetically equilibrate to. From numerical modeling, we know the ranges that cells separate to with varying magnetic contents per cell. Empirical separations needed to be carried out to correlate CD4 surface expression on human cell populations with iron oxide content and ultimately, micropillar pitch ranges and frequencies on the DMS cartridge. To achieve this, we procured pre-fractionated PBMC samples, which were split into T cells and monocytes (custom order from AllCells, Inc., Alameda, CA, USA). Using established assays from previous work [<a href="#B20-magnetochemistry-10-00081" class="html-bibr">20</a>], we stained the T cells with calcein AM red and the monocytes with calcein AM green to create a two-population mixture. After magnetic labeling and separation on a DMS chip, we fluorescently imaged the chip (a) to determine where the populations equilibrated under various magnetic labeling conditions and magnetic frequencies. (b) Image analysis enabled us to develop distributions of T cells vs. monocytes and determine that a 20 Hz separation frequency and a cutoff off pitch of 34 µm enabled maximum separability between the CD4<sup>Hi</sup> T cells and CD4<sup>Low</sup> monocytes/neutrophils. T cells (red) separated above the 34 µm pitch at a 20 Hz frequency, while the monocytes were mainly retained in the bottom, below the 34 µm pitch. This experiment set allowed us to optimize DMS separation parameters and cartridge designs to determine physical cartridge cutoff locations (extraction lock location and separation frequency).</p>
Full article ">Figure A5
<p>Quantification of cell viability and activation from digital magnetic sorting. (<b>A</b>) Cell viability after DMS sorting was ascertained via flow cytometry (SYTOX™ Advanced™ Dead Cell Stain Kit, Thermo Fisher Scientific, Carlsbad, CA, USA) and was found to be 97.3 ± 0.05% across three independent donors. (<b>B</b>) Analysis of cell activation on sorted T cells was also performed using flow 1 h after DMS processing, observing CD69 (early activation) and CD25 (late activation). We observed the population of CD69+ cells being slightly elevated (from 0.1% to 4.6%) compared to CD25, which showed no observable increase in CD25+ T cells. Some early-stage activation of cells is expected in any cell processing step but still shows a minor impact on cell activation.</p>
Full article ">
23 pages, 7974 KiB  
Article
Maize Phenotypic Parameters Based on the Constrained Region Point Cloud Phenotyping Algorithm as a Developed Method
by Qinzhe Zhu, Miaoyuan Bai and Ming Yu
Agronomy 2024, 14(10), 2446; https://doi.org/10.3390/agronomy14102446 - 21 Oct 2024
Viewed by 546
Abstract
As one of the world’s most crucial food crops, maize plays a pivotal role in ensuring food security and driving economic growth. The diversification of maize variety breeding is significantly enhancing the cumulative benefits in these areas. Precise measurement of phenotypic data is [...] Read more.
As one of the world’s most crucial food crops, maize plays a pivotal role in ensuring food security and driving economic growth. The diversification of maize variety breeding is significantly enhancing the cumulative benefits in these areas. Precise measurement of phenotypic data is pivotal for the selection and breeding of maize varieties in cultivation and production. However, in outdoor environments, conventional phenotyping methods, including point cloud processing techniques based on region growing algorithms and clustering segmentation, encounter significant challenges due to the low density and frequent loss of point cloud data. These issues substantially compromise measurement accuracy and computational efficiency. Consequently, this paper introduces a Constrained Region Point Cloud Phenotyping (CRPCP) algorithm that proficiently detects the phenotypic traits of multiple maize plants in sparse outdoor point cloud data. The CRPCP algorithm consists primarily of three core components: (1) a constrained region growth algorithm for effective segmentation of maize stem point clouds in complex backgrounds; (2) a radial basis interpolation technique to bridge gaps in point cloud data caused by environmental factors; and (3) a multi-level parallel decomposition strategy leveraging scene blocking and plant instances to enable high-throughput real-time computation. The results demonstrate that the CRPCP algorithm achieves a segmentation accuracy of 96.2%. When assessing maize plant height, the algorithm demonstrated a strong correlation with manual measurements, evidenced by a coefficient of determination R2 of 0.9534, a root mean square error (RMSE) of 0.4835 cm, and a mean absolute error (MAE) of 0.383 cm. In evaluating the diameter at breast height (DBH) of the plants, the algorithm yielded an R2 of 0.9407, an RMSE of 0.0368 cm, and an MAE of 0.031 cm. Compared to the PointNet point cloud segmentation method, the CRPCP algorithm reduced segmentation time by more than 44.7%. The CRPCP algorithm proposed in this paper enables efficient segmentation and precise phenotypic measurement of low-density maize multi-plant point cloud data in outdoor environments. This algorithm offers an automated, high-precision, and highly efficient solution for large-scale field phenotypic analysis, with broad applicability in precision breeding, agronomic management, and yield prediction. Full article
(This article belongs to the Section Precision and Digital Agriculture)
Show Figures

Figure 1

Figure 1
<p>(<b>a</b>) The three-dimensional point cloud data of corn plants and (<b>b</b>) low-density single corn plant point cloud, with the red box indicating regions of sparse data due to occlusion and missing points.</p>
Full article ">Figure 2
<p>A schematic diagram of the principle of the voxel denoising algorithm.</p>
Full article ">Figure 3
<p>Ground point separation schematic diagram.</p>
Full article ">Figure 4
<p>(<b>a</b>) Single corn plant point cloud data and (<b>b</b>) Point cloud data of three corn plants.</p>
Full article ">Figure 5
<p>CRPCP algorithm for segmentation of maize stem point cloud segmentation process, with arrows illustrating the segmentation workflow.</p>
Full article ">Figure 6
<p>CRPCP algorithm segmentation of corn stalk point cloud completion procedure, with arrows illustrating the completion process.</p>
Full article ">Figure 7
<p>Schematic diagram of the synchronous incremental segmentation completion process, with arrows indicating the key steps in the process.</p>
Full article ">Figure 8
<p>Parallel processing architecture diagram for the CRPCP algorithm, with arrows indicating the flow of execution steps.</p>
Full article ">Figure 9
<p>Schematic diagram of plant height calculation.</p>
Full article ">Figure 10
<p>Schematic diagram of breast diameter calculation.</p>
Full article ">Figure 11
<p>Detailed segmentation schematic diagram.</p>
Full article ">Figure 12
<p>Effect of segmented segmentation strategy on a corn plant.</p>
Full article ">Figure 13
<p>Effect of segmented segmentation strategy on multiple corn plants.</p>
Full article ">Figure 14
<p>Comparison of multiple corn plant point clouds before and after completion.</p>
Full article ">Figure 15
<p>(<b>a</b>) Scatterplot of plant height measurements obtained by the algorithm and manual measurements. and (<b>b</b>) Scatterplot of stem-breast diameter measurements obtained by the algorithm and manual measurements.</p>
Full article ">
22 pages, 6617 KiB  
Article
Contrasting Alleles of OsNRT1.1b Fostering Potential in Improving Nitrogen Use Efficiency in Rice
by Jonaliza L. Siangliw, Mathurada Ruangsiri, Cattarin Theerawitaya, Suriyan Cha-um, Wasin Poncheewin, Decha Songtoasesakul, Burin Thunnom, Vinitchan Ruanjaichon and Theerayut Toojinda
Plants 2024, 13(20), 2932; https://doi.org/10.3390/plants13202932 - 19 Oct 2024
Viewed by 777
Abstract
Nitrogen use efficiency (NUE) is important for the growth and development of rice and is significant in reducing the costs of rice production. OsNRT1.1b is involved in nitrate assimilation, and the alleles at position 21,759,092 on chromosome 10 clearly separate indica (Pathum Thani [...] Read more.
Nitrogen use efficiency (NUE) is important for the growth and development of rice and is significant in reducing the costs of rice production. OsNRT1.1b is involved in nitrate assimilation, and the alleles at position 21,759,092 on chromosome 10 clearly separate indica (Pathum Thani 1 (PTT1) and Homcholasit (HCS)) and japonica (Azucena and Leum Pua (LP)) rice varieties. Rice morphological and physiological traits were collected at three nitrogen levels (N0 = 0 kg ha−1, N7 = 43.75 kg ha−1, and N14 = 87.5 kg ha−1). Leaf and tiller numbers in PTT1 and HCS at N7 and N14 were two to three times higher than those at N0. At harvest, the biomass yield in PTT1 was the highest, while the total grain number in HCS was the maximum. The leaf widths and total chlorophyll contents (SPAD units) of Azucena and LP increased with nitrogen application as well as photosynthetic pigment parameters; for example, plant senescence reflectance indices (PSRIs), structure-insensitive pigment indices (SIPIs), and modified chlorophyll absorption ratio indices (MCARIs) were highly related in the japonica varieties. PTT1 and HCS, both carrying the A allele at OsNRT1.1b, had better NUE than Azucena and LP with the G allele. HCS, overall, had better NUE than PTT1. The translation to grain yield of assimilates was remarkable in PTT1 and HCS compared with Azucena and LP. In addition, HCS converted biomass for a 75% higher yield than PTT1. The ability of HCS to produce high yields was achieved even at N7 nitrogen fertilization, manifesting efficient use of nitrogen. Full article
(This article belongs to the Section Plant Nutrition)
Show Figures

Graphical abstract

Graphical abstract
Full article ">Figure 1
<p>Agronomic response of different rice varieties (PTT1, HCS, Azucena, and LP) at different times (30 DAS (seedling) (T1), 45 DAS (T2), 58 DAS (T3), 72 DAS (T4), 86 DAS (T5), and 100 DAS (grain filling) (T6)) and nitrogen levels (N0 (0 kg/ha), N7 (43.8 kg/ha), and N14 (87.5 kg/ha)). *, **, and *** represent significant differences at 0.05, 0.01, and 0.001 levels, respectively.</p>
Full article ">Figure 2
<p>Relative chlorophyll content (SPAD) and chlorophyll fluorescence parameter response of four rice varieties (PTT1, HCS, Azucena, and LP) at different growth stages (30 DAS (seedling) (T1), 45 DAS (T2), 58 DAS (T3), 72 DAS (T4), 86 DAS (T5), and 100 DAS (grain filling) (T6)) and nitrogen levels (N0 (0 kg/ha), N7 (43.8 kg/ha), and N14 (87.5 kg/ha)). *, **, and *** represent significant differences at 0.05, 0.01, and 0.001 levels, respectively.</p>
Full article ">Figure 3
<p>Yield and yield components of different rice varieties (PTT1, HCS, Azucena, and LP) at different nitrogen levels (N0 (0 kg/ha), N7 (43.8 kg/ha), and N14 (87.5 kg/ha)) at harvest. Small letters (abc) represent significant differences of the traits at different nitrogen levels for each variety.</p>
Full article ">Figure 4
<p>Agronomic response of rice varieties carrying G and A alleles at <span class="html-italic">OsNRT1.1b</span> at different nitrogen levels (N0 (0 kg/ha), N7 (43.8 kg/ha), and N14 (87.5 kg/ha)).</p>
Full article ">Figure 5
<p>Physiological responses measured by LI-COR 6400 XT, SPAD, and leaf chlorophyl fluorescence of rice varieties carrying the G and A alleles at <span class="html-italic">OsNRT1.1b</span> at different nitrogen levels (N0 (0 kg/ha), N7 (43.8 kg/ha), and N14 (87.5 kg/ha)).</p>
Full article ">Figure 6
<p>Hyperspectral reflectance parameters of rice varieties carrying the G and A alleles at <span class="html-italic">OsNRT1.1b</span> at different nitrogen levels (N0 (0 kg/ha), N7 (43.8 kg/ha), and N14 (87.5 kg/ha)).</p>
Full article ">Figure 7
<p>Yields and yield components of rice varieties carrying the G and A alleles at <span class="html-italic">OsNRT1.1b</span> at different nitrogen levels (N0 (0 kg/ha), N7 (43.8 kg/ha), and N14 (87.5 kg/ha)). Small letters represent significant differences of the trait at different nitrogen levels with the same allele, while capital letters represent significant differences of the trait at different nitrogen levels between two alleles.</p>
Full article ">Figure 8
<p>Correlation matrix of morphological and physiological traits at N7 (43.8 kg/ha) nitrogen level among the lines tested. Positive and negative correlations are shown with blue and red squares, respectively. Color shading is proportional to the correlation coefficients, with their values corresponding to the color intensity bar. Significance levels are indicated with asterisks (<span class="html-italic">p</span> &lt; 0.05 *; <span class="html-italic">p</span> &lt; 0.01 **; <span class="html-italic">p</span> &lt; 0.001 ***).</p>
Full article ">Figure 9
<p>Correlation matrix of morphological and physiological traits at the N14 (87.5 kg/ha) nitrogen level among the lines tested. Positive and negative correlations are shown with blue and red squares, respectively. Color shading is proportional to the correlation coefficients, with their values corresponding to the color intensity bar. Significance levels are indicated with asterisks (<span class="html-italic">p</span> &lt; 0.05 *; <span class="html-italic">p</span> &lt; 0.01 **; <span class="html-italic">p</span> &lt; 0.001 ***).</p>
Full article ">Figure 10
<p>Trait loading scores of the morphological and physiological traits for each principal component under two nitrogen levels: (<b>A</b>) N7 (43.8 kg/ha) and (<b>B</b>) N14 (87.5 kg/ha).</p>
Full article ">Figure 11
<p>Nitrogen use efficiency based on biomass (BIO), total grain number (TGN), and yield per plant (YLDPLT) of (<b>A</b>) varieties carrying the A allele (PTT1 and HCS) and the G allele (Azucena and LP) of <span class="html-italic">OsNRT1.1b</span> and (<b>B</b>) allele groups (A and G). Small letters represent significant differences in NUE at N7 and N14 in each rice variety (<b>A</b>,<b>B</b>) differences in NUE between alleles at the same nitrogen level.</p>
Full article ">
Back to TopTop