Abstract

Gene families contain genes that come from the same ancestor and have similar sequences and structures. They perform certain specific functions within and among different species. Currently, there is no complete process or platform for the rapid analysis of plant gene families. In this study, a comprehensive query and analysis platform of plant gene families, the Plant Gene Family Platform (PlantGF), was constructed. The platform is composed of four main parts: Search, Tools, Statistics and Auxiliary. A total of 2 909 580 gene family members were identified from 138 plant species in PlantGF. The data can be queried in the Search section through a user-friendly interface. A general process for gene family analysis, having nine steps, is provided. The platform also includes four online tools (HMM-Search, BLAST, MAFFT and HMMER) in the Tools section for useful additional analyses. The statistical analysis of the relevant gene families is shown on the Statistics page. Auxiliary pages are provided for data downloading. The datasets for all 138 plant species’ protein sequences and their gene families can be acquired on the Download page. A user’s manual and some useful links are displayed on the Manual and Links pages, respectively. To the best of our knowledge, PlantGF is the first comprehensive platform for studying plant gene families, and it will make important contributions to plant gene family-related research.

Database URL: http://biodb.sdau.edu.cn/PGF/index.html

Key points

  • Genome-level annotated protein sequences from 138 plant species have been used and analyzed in this research.

  • We searched all potential families and domains from the 138 species.

  • A total of seven useful online bioinformatic tools have been provided for the diverse analysis of the families.

  • Complete data are provided for users to start their own comprehensive analysis.

  • Based on all the above data, a web-friendly database has been constructed for all the researchers even without any family basic knowledge to make a rapid and effective analysis.

Introduction

A gene family is formed through the duplication and mutation of the same ancestor. Additionally, family members are defined as containing the same domains. For instance, WRKY genes, which all contain the W-R-K-Y domain, are important components of plant defense response-related signal transduction (1). Generally, these domains have conserved sequences that easily form stable three-dimensional structures, which then determine their particular function.

With improved sequencing technologies, the genome-wide sequencing of a series of important plant species has been completed, promoting research on plant genomics at the molecular level. Evolutionary biologists are now exploring the evolutionary laws of genomes using whole-genome data. To date, more than 500 plant species have been sequenced and released on public platforms (2). Currently, a huge number of gene families in one or more plants have been studied, but most of this research has been focused on families involved in specific plant characteristics, such as the SWEET family in pineapple (3), which is involved in the sugar transport process, the PPR family in tomato (4), which is involved in growth and development (4), and the WRKY family in tobacco (5), which is involved in stress resistance processes. However, there is no comprehensive platform to display and analyze all the available gene families of a plant species. This prompted the construction of a gene family database that would provide convenient access to data for all plant gene families.

Two comprehensive databases, Pfam (6) (http://pfam.xfam.org/) and InterPro (7) (http://www.ebi.ac.uk/interpro/), focus on gene family research. The Pfam database is a comprehensive platform for gene family processing, and it is dedicated to collecting specific domains and then identifying conserved domains using a Hidden Markov Model (HMM) algorithm. Interpro provides functional analyses of proteins through classification and domain and important site predictions. However, these databases mainly focus on gene families and their functions. There is currently no complete plant gene family study and self-analysis database. Therefore, it is necessary to establish an omnibus platform of plant gene families. In this study, the plant gene family database was developed to query all the gene families of sequenced plants and their functional annotations, and it provides some analytical tools and useful links. The database is a resource and analysis platform through which plant gene families can be well studied, as well as evolutionary relationships.

Platform content and web interface

Plant Gene Family (PlantGF) is a database committed to the gene families of sequenced plants, and it provides family analysis-related online tools. The database contains four main components: Search, Tools, Statistics and Auxiliary. Search and Tools are the two core components of the platform, and there are some necessary related modules (i.e. Manual, About us and Links). The result is an open and user-friendly web interface (Figure 1).

Main PlantGF web page. (A) PlantGF homepage: provides quick entry paths to all main parts. (B) Search: contains 2 909 580 gene family members and their specific detail annotations. (C) Tools: consists of nine steps in gene family analysis. Among that, four online tools (HMM-Search, BLAST, MAFFT and HMMER) also exist. (D) Statistics: statistics among these gene families. (E) Datasets Download: provides download links of 138 plant species’ protein sequences and their gene families.
Figure 1.

Main PlantGF web page. (A) PlantGF homepage: provides quick entry paths to all main parts. (B) Search: contains 2 909 580 gene family members and their specific detail annotations. (C) Tools: consists of nine steps in gene family analysis. Among that, four online tools (HMM-Search, BLAST, MAFFT and HMMER) also exist. (D) Statistics: statistics among these gene families. (E) Datasets Download: provides download links of 138 plant species’ protein sequences and their gene families.

PlantGF’s main components are as follows: (i) Search; this section contains detailed data on 138 plant gene families. These gene families can be queried using species name, family type, species family ID, family name, Pfam accession and Description. In addition, some exhaustive gene annotations, including Pfam (6), Prosite (8), EMBL (9), KEGG (10) and GO (11) are displayed on secondary webpages. (ii) Tools; a series of popular and convenient tools for gene family analyses are shown. There are nine common steps: Data Acquisition, Family Identification, Physicochemical Property Analysis, Structural Analysis, Phylogenetic Analysis, Collinearity Analysis, Annotation Analysis, Gene Location and Expression Patterns. Among them, four specific online tools applicable to this database, HMM (12), BLAST (13), MAFFT (14) and HMMER (12), are marked in red on this page. (iii) Statistics; this section contains the numbers of family members and their specific distributions. (iv) Auxiliary; this section provides the protein sequences and gene families found in 138 plant species. Some useful links and the user’s manual are also available. In addition, gene family types and member numbers in different species are displayed and can be downloaded in file tree form on the Statistics page.

Gene family

All the gene Family, Domain, Coiled-coil, Disordered, Motif and Repeat have been identified in the 138 plant species’ gene sequences attained from public platforms (Figure 2A). In total, 2 909 580 gene family members were identified, and some necessary annotations were developed for each gene. Approximately 80% of plant species contain the PPR gene family, followed by the Mito_carr gene family (Figure 2B). Triticum aestivum has the largest number of gene families at 121 667, whereas Ostreococcus lucimarinus has the smallest number at 3686 (Figure 2C).

Statistics of plant gene families. (A) The composition structure of 138 plant species’ genes. (B) Statistics of the largest number of gene families in each species. (C) The number of gene families of each species.
Figure 2.

Statistics of plant gene families. (A) The composition structure of 138 plant species’ genes. (B) Statistics of the largest number of gene families in each species. (C) The number of gene families of each species.

All the gene family datasets and their annotations are stored on the Search page. Users can search these data using checkboxes: Species (138 plant species), Type (Family, Domain, Coiled-coil, Disordered, Motif and Repeat), PGF/short name-gene id (PGF-ID), Family name, Pfam accession and Description. The resulting specified data will appear in a dynamic table after users click ‘Search’ in the current page. The annotation page contains three main parts: Gene sequences, Family/Domain information and Family Gene annotations. (i) Gene sequences; the detailed sequence information and download link are shown. (ii) Family/Domain information; using the annotation data from Pfam (6) and Prosite (8), basic information and a detailed description of the family are provided. In addition, the HMM can be downloaded. (iii) Family Gene annotations; detailed annotations for the queried gene, including SwissProt, EMBL (9), KEGG (10) and GO (11) results are shown. Furthermore, a table of species names, including their short and common names, is provided on the Search page to facilitate user queries.

Tools

Nine analyses-related steps powered by 26 software programs are shown on the Tools page. The nine steps are Data Acquisition, Family Identification, Physicochemical Property Analysis, Structural Analysis, Phylogenetic Analysis, Collinearity Analysis, Annotation Analysis, Gene Location and Expression Patterns (Table 1).

Table 1.

Nine steps of gene family analysis

StepsNameTools
1Data AcquisitionExpression Patterns; Homologous Families (Genera) Database; Single species Database
2Family IdentificationHMM-Search (12); BLAST (13); MAFFT (14); HMMER (12); CDD-NCBI (24); DNAman; MEME (25); SMART (26)
3Physicochemical PropertiesCompute pI/Mw tool (27); TMHMM (28); SingnalP (29); CELLO (30)
4Structural AnalysisGSDS (31); PDB (32); PredictProtein (33)
5Phylogenetic AnalysisMEGA (34); Evolview (35)
6Collinearity AnalysisCricos (36); MCscanX (37)
7Annotation AnalysisGene Ontology (11); Wego (38); KEGG (10); KAAS (39)
8Gene locationMapChart (40)
9Expression PatternsArrayExpress (41); NCBI-GEO (42)
StepsNameTools
1Data AcquisitionExpression Patterns; Homologous Families (Genera) Database; Single species Database
2Family IdentificationHMM-Search (12); BLAST (13); MAFFT (14); HMMER (12); CDD-NCBI (24); DNAman; MEME (25); SMART (26)
3Physicochemical PropertiesCompute pI/Mw tool (27); TMHMM (28); SingnalP (29); CELLO (30)
4Structural AnalysisGSDS (31); PDB (32); PredictProtein (33)
5Phylogenetic AnalysisMEGA (34); Evolview (35)
6Collinearity AnalysisCricos (36); MCscanX (37)
7Annotation AnalysisGene Ontology (11); Wego (38); KEGG (10); KAAS (39)
8Gene locationMapChart (40)
9Expression PatternsArrayExpress (41); NCBI-GEO (42)
Table 1.

Nine steps of gene family analysis

StepsNameTools
1Data AcquisitionExpression Patterns; Homologous Families (Genera) Database; Single species Database
2Family IdentificationHMM-Search (12); BLAST (13); MAFFT (14); HMMER (12); CDD-NCBI (24); DNAman; MEME (25); SMART (26)
3Physicochemical PropertiesCompute pI/Mw tool (27); TMHMM (28); SingnalP (29); CELLO (30)
4Structural AnalysisGSDS (31); PDB (32); PredictProtein (33)
5Phylogenetic AnalysisMEGA (34); Evolview (35)
6Collinearity AnalysisCricos (36); MCscanX (37)
7Annotation AnalysisGene Ontology (11); Wego (38); KEGG (10); KAAS (39)
8Gene locationMapChart (40)
9Expression PatternsArrayExpress (41); NCBI-GEO (42)
StepsNameTools
1Data AcquisitionExpression Patterns; Homologous Families (Genera) Database; Single species Database
2Family IdentificationHMM-Search (12); BLAST (13); MAFFT (14); HMMER (12); CDD-NCBI (24); DNAman; MEME (25); SMART (26)
3Physicochemical PropertiesCompute pI/Mw tool (27); TMHMM (28); SingnalP (29); CELLO (30)
4Structural AnalysisGSDS (31); PDB (32); PredictProtein (33)
5Phylogenetic AnalysisMEGA (34); Evolview (35)
6Collinearity AnalysisCricos (36); MCscanX (37)
7Annotation AnalysisGene Ontology (11); Wego (38); KEGG (10); KAAS (39)
8Gene locationMapChart (40)
9Expression PatternsArrayExpress (41); NCBI-GEO (42)

In addition, four online tools are provided in the Tools section to promote analyses among gene families. HMM-Search provides a simple way to search the HMMs Users just need to input keywords in the textbox and click ‘Submit’. The results will be displayed in a new window. The BLAST software was mainly developed using PlantGF, and the BLAST library consists of 138 plants. Users can input query sequences and select appropriate parameters to inform their results. MAFFT is a multiple sequence alignment program for Unix-like operating systems that can be used for preparing multiple sequence alignment files to develop phylogenetic trees and for HMMER Build. Two main HMMER modules, HMMER Build and HMMER Search, are also provided. HMMs can be developed using HMMER Build and the results of MAFFT. HMMER Search requires users to upload the target HMM and choose plant species to identify their target genes. All the results can be visualized online or can be download to a local computer through the results page.

Statistics

Statistical analyses are displayed using the different charts available in this section. The construction of a gene family in 138 species and the gene families with the most members in each species are shown in two pie charts. Their detailed information appears when the mouse is passed over the target area. Statistics on gene family types and numbers of different species are shown in the form of tree files. The species are classified using the initials of their scientific names. The number of gene families in each species is shown using a dynamic histogram at the bottom of the page. Users can also turn the mouse wheel to see specific species.

Materials and methods

Dataset collection

Identified protein sequences were obtained from species-specific databases and public comprehensive platforms. In total, 13 Brassicaceae family species, including Aethionema arabicum, Arabidopsis thaliana and Brassica napus, were downloaded from the Brassica Database (http://brassicadb.org/brad/index.php) (15); 12 Solanaceae family species, including Nicotiana tabacum, Solanum tuberosum and Solanum lycopersicum, were downloaded from the Sol Genomics Network (https://solgenomics.net/) (16); 4 Rosaceae family species, including Prunus mume, Pyrus bretschneideri and Malus × domestica, were downloaded from Genome Database for Rosaceae (https://www.rosaceae.org/) (17); 16 Grass family species, including Aegilops tauschii, Hordeum vulgare and Oryza sativa, were downloaded from Gramene (http://gramene.org/) (18) and 3 Cucurbitaceae family species, Citrullus lanatus, Cucumis melo and Cucumis sativus, were downloaded from Cucurbit Genomics Database (http://cucurbitgenomics.org/) (19). The remaining species came from comprehensive plant databases, such as Phytozome (http://www.phytozome.net) (20), PlantGDB (http://www.plantgdb.org/) (21), NCBI (https://www.ncbi.nlm.nih.gov/), EnsemblPlants (http://plants.ensembl.org/index.html) and PMDbase (http://www.sesame-bioinfo.org/PMDBase/) (22).

Gene family identity

In this study, Perl script-based Pfamscan software was used to input each protein sequence in a Linux environment with the default parameters. All the data processing and statistics were performed in Perl script, R script and Echarts.

Gene family annotation

Currently, different annotation databases, like Pfam (6), GO (11), KEGG (10), Uniprot (23) and Prosite (8), contain massive amounts of accurate annotation data, and data can be exchanged among these databases. A one-to-one correspondence for these gene families was annotated using the Pfam accession number of each gene family. Furthermore, it was not possible to obtain accurate annotations for every sequence in every plant. To help the user explore an unknown gene with our platform, we chose genes of well-studied model plants, such as A. thaliana and O. sativa, as the annotation sources for the gene family. Using this information, researchers may infer the function and origin of their research targets.

Database implementation

The Python Web framework is popular for constructing databases. First, the detailed annotations of gene families were stored in MySQL, in which data manipulation and maintenance were also performed. Then, uWSGI was used together with HTML and Bootstrap to construct the users’ access interface. Additionally, Flask, BioPython, Perl Scripts, Echarts and Javascript were all required to connect MySQL and uWSGI.

Conclusions

PlantGF is a comprehensive platform that was developed for the study of plant gene families. It provides 2 909 580 gene family members and their specific detail annotations from 138 plants. Furthermore, the incorporation of several useful tools makes it easy for users unfamiliar with bioinformatics to perform plant family-related scientific research. The platform will be updated continuously as new plant sequences are generated and new bioinformatics tools emerge.

Funding

This research work was supported by the Foundation of Innovation Team Project for Modern Agricultural Industrious Technology System of Shandong Province (SDAIT-25-01). And we thank Supercomputing Center in Shandong Agricultural University for technical support.

Conflict of interest

None declared.

References

1.

Phukan
U.J.
,
Jeena
G.S.
and
Shukla
R.K.
(
2016
)
WRKY transcription factors: molecular regulation and stress responses in plants
.
Front. Plant Sci.
,
7
, 760.doi: .

2.

Song
X.
,
Liu
Z.
,
Wan
H.
 et al.  (
2021
)
Editorial: comparative genomics and functional genomics analyses in plants
.
Front. Genet.
,
12
, 687966.doi: .

3.

Guo
C.
,
Li
H.
,
Xia
X.
 et al.  (
2018
)
Functional and evolution characterization of SWEET sugar transporters in Ananas comosus
.
Biochem. Biophys. Res. Commun.
,
496
,
407
414
.

4.

Ding
A.
,
Li
L.
,
Qu
X.
 et al.  (
2014
)
Genome-wide identification and bioinformatic analysis of PPR gene family in tomato
.
Yi Chuan
,
36
,
77
84
.

5.

Schluttenhofer
C.
and
Yuan
L.
(
2015
)
Regulation of specialized metabolism by WRKY transcription factors
.
Plant Physiol.
,
167
,
295
306
.

6.

El-Gebali
S.
,
Mistry
J.
,
Bateman
A.
 et al.  (
2019
)
The Pfam protein families database in 2019
.
Nucleic Acids Res.
,
47
,
D427
D432
.

7.

Mitchell
A.L.
,
Attwood
T.K.
,
Babbitt
P.C.
 et al.  (
2019
)
InterPro in 2019: improving coverage, classification and access to protein sequence annotations
.
Nucleic Acids Res.
,
47
,
D351
D360
.

8.

Hulo
N.
,
Bairoch
A.
,
Bulliard
V.
 et al.  (
2006
)
The PROSITE database
.
Nucleic Acids Res.
,
34
,
D227
D230
.

9.

Madeira
F.
,
Park
Y.M.
,
Lee
J.
 et al.  (
2019
)
The EMBL-EBI search and sequence analysis tools APIs in 2019
.
Nucleic Acids Res.
,
47
,
W636
W641
.

10.

Kanehisa
M.
,
Furumichi
M.
,
Tanabe
M.
 et al.  (
2017
)
KEGG: new perspectives on genomes, pathways, diseases and drugs
.
Nucleic Acids Res.
,
45
,
D353
D361
.

11.

The Gene Ontology Consortium
. (
2019
)
The gene ontology resource: 20 years and still GOing strong
.
Nucleic Acids Res.
,
47
,
D330
D338
.

12.

Eddy
S.R.
(
2011
)
Accelerated profile HMM searches
.
PLoS Comput. Biol.
,
7
, e1002195.

13.

Zhang
Z.
,
Schwartz
S.
,
Wagner
L.
 et al.  (
2000
)
A greedy algorithm for aligning DNA sequences
.
J. Comput. Biol.
,
7
,
203
214
.

14.

Katoh
K.
and
Standley
D.M.
(
2013
)
MAFFT multiple sequence alignment software version 7: improvements in performance and usability
.
Mol. Biol. Evol.
,
30
,
772
780
.

15.

Wang
X.
,
Wu
J.
,
Liang
J.
 et al.  (
2015
)
Brassica database (BRAD) version 2.0: integrating and mining Brassicaceae species genomic resources
.
Database: the journal of biological databases and curation
,
2015
, bav093.doi: .

16.

Fernandez-Pozo
N.
,
Menda
N.
,
Edwards
J.D.
 et al.  (
2015
)
The Sol Genomics Network (SGN)—from genotype to phenotype to breeding
.
Nucleic Acids Res.
,
43
,
D1036
D1041
.

17.

Jung
S.
,
Lee
T.
,
Cheng
C.-H.
 et al.  (
2019
)
15 years of GDR: new data and functionality in the genome database for rosaceae
.
Nucleic Acids Res.
,
47
,
D1137
D1145
.

18.

Tello-Ruiz
M.K.
,
Naithani
S.
,
Stein
J.C.
 et al.  (
2018
)
Gramene 2018: unifying comparative genomics and pathway resources for plant research
.
Nucleic Acids Res.
,
46
,
D1181
D1189
.

19.

Zheng
Y.
,
Wu
S.
,
Bai
Y.
 et al.  (
2019
)
Cucurbit Genomics Database (CuGenDB): a central portal for comparative and functional genomics of cucurbit crops
.
Nucleic Acids Res.
,
47
,
D1128
D1136
.

20.

Goodstein
D.M.
,
Shu
S.
,
Howson
R.
 et al.  (
2012
)
Phytozome: a comparative platform for green plant genomics
.
Nucleic Acids Res.
,
40
,
D1178
D1186
.

21.

Duvick
J.
,
Fu
A.
,
Muppirala
U.
 et al.  (
2008
)
PlantGDB: a resource for comparative plant genomics
.
Nucleic Acids Res.
,
36
,
D959
D965
.

22.

Yu
J.
,
Dossa
K.
,
Wang
L.
 et al.  (
2017
)
PMDBase: a database for studying microsatellite DNA and marker development in plants
.
Nucleic Acids Res.
,
45
,
D1046
D1053
.

23.

The UniProt Consortium
. (
2017
)
UniProt: the universal protein knowledgebase
.
Nucleic Acids Res.
,
45
,
D158
D169
.

24.

Marchler-Bauer
A.
,
Derbyshire
M.K.
,
Gonzales
N.R.
 et al.  (
2015
)
CDD: NCBI’s conserved domain database
.
Nucleic Acids Res.
,
43
,
D222
D226
.

25.

Bailey
T.L.
,
Johnson
J.
,
Grant
C.E.
 et al.  (
2015
)
The MEME suite
.
Nucleic Acids Res.
,
43
,
W39
W49
.

26.

Letunic
I.
,
Khedkar
S.
and
Bork
P.
(
2021
)
SMART: recent updates, new developments and status in 2020
.
Nucleic Acids Res.
,
49
,
D458
D460
.

27.

Wilkins
M.R.
,
Gasteiger
E.
,
Bairoch
A.
 et al.  (
1999
)
Protein identification and analysis tools in the ExPASy server
.
Methods Mol. Biol.
,
112
,
531
552
.

28.

Krogh
A.
,
Larsson
B.
,
von Heijne
G.
 et al.  (
2001
)
Predicting transmembrane protein topology with a hidden Markov model: application to complete genomes
.
J. Mol. Biol.
,
305
,
567
580
.

29.

Almagro Armenteros
J.J.
,
Tsirigos
K.D.
,
Sønderby
C.K.
 et al.  (
2019
)
SignalP 5.0 improves signal peptide predictions using deep neural networks
.
Nat. Biotechnol.
,
37
,
420
423
.

30.

Bernstein
M.N.
,
Ma
Z.
,
Gleicher
M.
 et al.  (
2021
)
CellO: comprehensive and hierarchical cell type classification of human cells with the Cell Ontology
.
iScience
,
24
, 101913.

31.

Hu
B.
,
Jin
J.
,
Guo
A.-Y.
 et al.  (
2015
)
GSDS 2.0: an upgraded gene feature visualization server
.
Bioinformatics (Oxford, England)
,
31
,
1296
1297
.

32.

Berman
H.M.
,
Westbrook
J.
,
Feng
Z.
 et al.  (
2000
)
The protein data bank
.
Nucleic Acids Res.
,
28
,
235
242
.

33.

Bernhofer
M.
,
Dallago
C.
,
Karl
T.
 et al.  (
2021
)
PredictProtein - predicting protein structure and function for 29 years
.
Nucleic Acids Res.
,
49
,
W535
W540
.

34.

Kumar
S.
,
Stecher
G.
,
Li
M.
 et al.  (
2018
)
MEGA X: molecular evolutionary genetics analysis across computing platforms
.
Mol. Biol. Evol.
,
35
,
1547
1549
.

35.

Subramanian
B.
,
Gao
S.
,
Lercher
M.J.
 et al.  (
2019
)
Evolview v3: a webserver for visualization, annotation, and management of phylogenetic trees
.
Nucleic Acids Res.
,
47
,
W270
W275
.

36.

Krzywinski
M.
,
Schein
J.
,
Birol
I.
 et al.  (
2009
)
Circos: an information aesthetic for comparative genomics
.
Genome Res.
,
19
,
1639
1645
.

37.

Wang
Y.
,
Tang
H.
,
Debarry
J.D.
 et al.  (
2012
)
MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity
.
Nucleic Acids Res.
,
40
, e49.

38.

Ye
J.
,
Zhang
Y.
,
Cui
H.
 et al.  (
2018
)
WEGO 2.0: a web tool for analyzing and plotting GO annotations, 2018 update
.
Nucleic Acids Res.
,
46
,
W71
W75
.

39.

Moriya
Y.
,
Itoh
M.
,
Okuda
S.
 et al.  (
2007
)
KAAS: an automatic genome annotation and pathway reconstruction server
.
Nucleic Acids Res.
,
35
,
W182
W185
.

40.

Voorrips
R.E.
(
2002
)
MapChart: software for the graphical presentation of linkage maps and QTLs
.
J. Hered.
,
93
,
77
78
.

41.

Athar
A.
,
Füllgrabe
A.
,
George
N.
 et al.  (
2019
)
ArrayExpress update - from bulk to single-cell expression data
.
Nucleic Acids Res.
,
47
,
D711
D715
.

42.

Barrett
T.
,
Wilhite
S.E.
,
Ledoux
P.
 et al.  (
2013
)
NCBI GEO: archive for functional genomics data sets—update
.
Nucleic Acids Res.
,
41
,
D991
D995
.

Author notes

Jiaxuan Li and Shuai Yang contributed equally to this work and share first authorship.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0/), which permits unrestricted reuse, distribution, and reproduction in any medium, provided the original work is properly cited.