Gene Structure Prediction Research Articles

Genome-wide identification of gene families related to miRNA biogenesis in Mangifera indica L. and their possible role during heat stress.

Mango is a popular tropical fruit that requires quarantine hot water treatment (QHWT) for postharvest sanitation, which can cause abiotic stress. Plants have various defense mechanisms to cope with stress; miRNAs mainly regulate the expression of these defense responses. Proteins involved in the biogenesis of miRNAs include DICER-like (DCL), ARGONAUTE (AGO), HYPONASTIC LEAVES 1 (HYL1), SERRATE (SE), HUA ENHANCER1 (HEN1), HASTY (HST), and HEAT-SHOCK PROTEIN 90 (HSP90), among others. According to our analysis, the mango genome contains five DCL, thirteen AGO, six HYL, two SE, one HEN1, one HST, and five putative HSP90 genes. Gene structure prediction and domain identification indicate that sequences contain key domains for their respective gene families, including the RNase III domain in DCL and PAZ and PIWI domains for AGOs. In addition, phylogenetic analysis indicates the formation of clades that include the mango sequences and their respective orthologs in other flowering plant species, supporting the idea these are functional orthologs. The analysis of cis-regulatory elements of these genes allowed the identification of MYB, ABRE, GARE, MYC, and MeJA-responsive elements involved in stress responses. Gene expression analysis showed that most genes are induced between 3 to 6 h after QHWT, supporting the early role of miRNAs in stress response. Interestingly, our results suggest that mango rapidly induces the production of miRNAs after heat stress. This research will enable us to investigate further the regulation of gene expression and its effects on commercially cultivated fruits, such as mango, while maintaining sanitary standards.

Chromosome-level genome assembly of Hippophae rhamnoides variety

Fructus hippophae (Hippophae rhamnoides spp. mongolica×Hippophae rhamnoides sinensis), a hybrid variety of sea buckthorn that Hippophae rhamnoides spp. mongolica serves as the female parent and Hippophae rhamnoidessinensis serves as the male parent, is a traditional plant with great potentials of economic and medical values. Herein, we gained a chromosome-level genome of Fructus hippophae about 918.59 Mb, with the scaffolds N50 reaching 83.65 Mb. Then, we anchored 440 contigs with 97.17% of the total genome sequences onto 12 pseudochromosomes. Next, de-novo, homology and transcriptome assembly strategies were adopted for gene structure prediction. This predicted 36475 protein-coding genes, of which 36226 genes could be functionally annotated. Simultaneously, various strategies were used for quality assessment, both the complete BUSCO value (98.80%) and the mapping rate indicated the high assembly quality. Repetitive elements, which occupied 63.68% of the genome, and 1483600 bp of non-coding RNA were annotated. Here, we provide genomic information on female plants of a popular variety, which can provide data for pan-genomic construction of sea buckthorn and for the resolution of the mechanism of sex differentiation.

Genome‐wide study of OsWRKY gene in Oryza sativa subsp. japonica

AbstractWRKY transcription factor (TF) family is well known to govern essential physiological functioning as well as regulate plant response to biotic and abiotic stress. In this study, we have identified 108 OsWRKY genes in the genome of Oryza sativa subsp. japonica, using the updated genomic data from the Rice Annotation Project Database and Oryzabase, which were further used to conduct the phylogenetic study, motif analysis, gene structure analysis, chromosomal mapping, and prediction of sub‐cellular localization. The multiple sequence alignment OsWRKY proteins revealed the presence of nine different types of alterations in the conserved heptapeptide sequence WRKYGQK associated with 19 OsWRKY genes. Physiochemical analysis discloses the hydrophobic amino acid‐rich, thermally stable, and polar nature of OsWRKY proteins. These genes were noted as highly conserved between the two cultivated sub‐species of Oryza sativa, that is Indica and japonica type. Additionally, from motif analysis, we have found a new motif, which was categorized as hAT family C‐terminal dimerization region associated with four members of group IIc. We have identified 21 stress‐responsive OsWRKY genes, and their significance to the different biotic and abiotic stress‐mediated cascades was further evaluated by analysing 1500 kb upstream sequences and this disclosed the presence of important phytohormone‐responsive cis‐elements in the OsWRKY gene, suggesting its direct involvement in defence against a wide range of external stressors and these 21 OsWRKY genes are tentatively listed as possible candidates for more study.

How to overcome constraints imposed by microsporidian genome features to ensure gene prediction?

Since the advent of sequencing techniques and due to their continuous evolution, it has become easier and less expensive to obtain the complete genome sequence of any organism. Nevertheless, to elucidate all biological processes governing organism development, quality annotation is essential. In genome annotation, predicting gene structure is one of the most important and captivating challenges for computational biology. This aspect of annotation requires continual optimization, particularly for genomes as unusual as those of microsporidia. Indeed, this group of fungal-related parasites exhibits specific features (highly reduced gene sizes, sequences with high rate of evolution) linked to their evolution as intracellular parasites, requiring the implementation of specific annotation approaches to consider all these features. This review aimed to outline these characteristics and to assess the increasingly efficient approaches and tools that have enhanced the accuracy of gene prediction for microsporidia, both in terms of sensitivity and specificity. Subsequently, a final part will be dedicated to postgenomic approaches aimed at reinforcing the annotation data generated by prediction software. These approaches include the characterization of other understudied genes, such as those encoding regulatory noncoding RNAs or very small proteins, which also play crucial roles in the life cycle of these microorganisms.

An ancient bacterial zinc acquisition system identified from a cyanobacterial exoproteome.

Bacteria have developed fine-tuned responses to cope with potential zinc limitation. The Zur protein is a key player in coordinating this response in most species. Comparative proteomics conducted on the cyanobacterium Anabaena highlighted the more abundant proteins in a zur mutant compared to the wild type. Experimental evidence showed that the exoprotein ZepA mediates zinc uptake. Genomic context of the zepA gene and protein structure prediction provided additional insights on the regulation and putative function of ZepA homologs. Phylogenetic analysis suggests that ZepA represents a primordial system for zinc acquisition that has been conserved for billions of years in a handful of species from distant bacterial lineages. Furthermore, these results show that Zur may have been one of the first regulators of the FUR family to evolve, consistent with the scarcity of zinc in the ecosystems of the Archean eon.

Genome-wide investigation and expression analysis of OSCA gene family in response to abiotic stress in alfalfa.

Alfalfa is an excellent leguminous forage crop that is widely cultivated worldwide, but its yield and quality are often affected by drought and soil salinization. Hyperosmolality-gated calcium-permeable channel (OSCA) proteins are hyperosmotic calcium ion (Ca2+) receptors that play an essential role in regulating plant growth, development, and abiotic stress responses. However, no systematic analysis of the OSCA gene family has been conducted in alfalfa. In this study, a total of 14 OSCA genes were identified from the alfalfa genome and classified into three groups based on their sequence composition and phylogenetic relationships. Gene structure, conserved motifs and functional domain prediction showed that all MsOSCA genes had the same functional domain DUF221. Cis-acting element analysis showed that MsOSCA genes had many cis-regulatory elements in response to abiotic or biotic stresses and hormones. Tissue expression pattern analysis demonstrated that the MsOSCA genes had tissue-specific expression; for example, MsOSCA12 was only expressed in roots and leaves but not in stem and petiole tissues. Furthermore, RT-qPCR results indicated that the expression of MsOSCA genes was induced by abiotic stress (drought and salt) and hormones (JA, SA, and ABA). In particular, the expression levels of MsOSCA3, MsOSCA5, MsOSCA12 and MsOSCA13 were significantly increased under drought and salt stress, and MsOSCA7, MsOSCA10, MsOSCA12 and MsOSCA13 genes exhibited significant upregulation under plant hormone treatments, indicating that these genes play a positive role in drought, salt and hormone responses. Subcellular localization results showed that the MsOSCA3 protein was localized on the plasma membrane. This study provides a basis for understanding the biological information and further functional analysis of the MsOSCA gene family and provides candidate genes for stress resistance breeding in alfalfa.

GINGER: an integrated method for high-accuracy prediction of gene structure in higher eukaryotes at the gene and exon level.

The prediction of gene structure within the genome sequence is the starting point of genome analysis, and its accuracy has a significant impact on the quality of subsequent analyses. Gene structure prediction is roughly divided into RNA-Seq-based methods, ab initio-based methods, homology-based methods, and the integration of individual prediction methods. Integrated methods are mainstream in recent genome projects because they improve prediction accuracy by combining or taking the best individual prediction findings; however, adequate prediction accuracy for eukaryotic species has not yet been achieved. Therefore, we developed an integrated tool, GINGER, that solves various issues related to gene structure prediction in higher eukaryotes. By handling artefacts in alignments of RNA and protein sequences, reconstructing gene structures via dynamic programming with appropriately weighted and scored exon/intron/intergenic regions, and applying different prediction processes and filtering criteria to multi-exon and single-exon genes, we achieved a significant improvement in accuracy compared to the existing integration methods. The feature of GINGER is its high prediction accuracy at the gene and exon levels, which is pronounced for species with more complex gene architectures. GINGER is implemented using Nextflow, which allows for the efficient and effective use of computing resources.

Flagellin and mannitol modulate callose biosynthesis and deposition in soybean seedlings.

Callose is a polymer deposited on the cell wall and is necessary for plant growth and development. Callose is synthesized by genes from the glucan synthase-like family (GSL) and dynamically responds to various types of stress. Callose can inhibit pathogenic infection, in the case of biotic stresses, and maintain cell turgor and stiffen the plant cell wall in abiotic stresses. Here, we report the identification of 23 GSL genes (GmGSL) in the soybean genome. We performed phylogenetic analyses, gene structure prediction, duplication patterns, and expression profiles on several RNA-Seq libraries. Our analyses show that WGD/Segmental duplication contributed to expanding this gene family in soybean. Next, we analyzed the callose responses in soybean under abiotic and biotic stresses. The data show that callose is induced by both osmotic stress and flagellin 22 (flg22) and is related to the activity of β-1,3-glucanases. By using RT-qPCR, we evaluated the expression of GSL genes during the treatment of soybean roots with mannitol and flg22. The GmGSL23 gene was upregulated in seedlings treated with osmotic stress or flg22, showing the essential role of this gene in the soybean defense response to pathogenic organisms and osmotic stress. Our results provide an important understanding of the role of callose deposition and regulation of GSL genes in response to osmotic stress and flg22 infection in soybean seedlings.

Comparative genomics studies on the stk gene family in vertebrates: From the bighead carp (Hypophthalmichthys nobilis) genome

The mammalian sterile 20-like (MST) family belongs to the serine/threonine protein kinase (STK) superfamily and participates in a variety of biological processes, such as cell apoptosis, polarity, migration, immune regulation, inflammatory responses, and cancer. In the economically important bighead carp (Hypophthalmichthys nobilis), the STK gene family and immune-related biological functions may be helpful in increasing its economic yield. However, the comprehensive role of STKs in the bighead carp remains unclear. In this study, the five stk sequences from the bighead carp were divided into two classes: stk3/4 and stk24/25/26. Gene structure and motif prediction analyses confirmed that stk is conserved in the bighead carp. Compared to 26 other vertebrate species, teleosts (including bighead carp) possess more stk members because of teleost-specific whole-genome duplication. Synteny analysis revealed that stk3, stk24, stk25, and stk26 have been relatively conserved in bighead carp during evolution. Meanwhile, stk4 was lost in most Cyprinid species, including bighead carp, during evolution. RNA-seq data revealed that STK expression was associated with various pathogens, and the expression of these STKs (Hnstk3, Hnstk24a, Hnstk24b, Hnstk25, and Hnstk26) was different in seven tissues of bighead carp. In addition, we showed that STK expression levels were dramatically altered in the head kidney and that stk24 was involved in defense against Aeromonas hydrophila. This study provides a molecular basis for the analysis of stk function in bighead carp, and can be used as a reference for further phylogenomics.

Genome-wide identification and analysis of the evolution and expression pattern of the HVA22 gene family in three wild species of tomatoes.

Wild tomato germplasm is a valuable resource for improving biotic and abiotic stresses in tomato breeding. The HVA22 is widely present in eukaryotes and involved in growth and development as well as stress response, such as cold, salt, drought, and biotic stress. In the present study, we identified 45 HVA22 genes in three wild species of tomatoes. The phylogenetic relationships, gene localization to chromosomes, gene structure, gene collinearity, protein interactions, and cis-acting element prediction of all 45 HVA22 genes (14 in Solanum pennellii, 15 in S. pimpinellifolium, and 16 in S. lycopersicoides) were analyzed. The phylogenetic analysis showed that the all HVA22 proteins from the family Solanaceae were divided into three branches. The identified 45 HVA22 genes were grouped into four subfamilies, which displayed similar number of exons and expanded in a fragmentary replication manner. The distribution of HVA22 genes on the chromosomes of the three wild tomato species was also highly similar. RNA-seq and qRT-PCR revealed that HVA22 genes were expressed in different tissues and induced by drought, salt, and phytohormone treatments. These results might be useful for explaining the evolution, expression patterns, and functional divergence of HVA22 genes in Lycopersicon.

Transcriptome analysis and gene expression analysis related to salinity-alkalinity and low temperature adaptation of Triplophysa yarkandensis.

T. yarkandensis is a common species of Triplophysa, and it is distributed in Shule river of Hexi Corridor, of Gansu province in China. In order to enrich gene database resources and explore the environment adaptation of T. yarkandensis, fifteen tissues were collected from three adult T. yarkandensis for transcriptome sequencing and de novo assembly. Nine major international gene annotation databases (NR, COG, egg_NOG, TrEMBL, Pfam, KOG, Swiss prot, KEGG and Gene Ontology) were utilized to annotate unigenes. A detailed study was conducted to explore the gene expression and the differentially expressed genes among five tissues (brain, heart, kidney, liver and spleen). In addition, the current study showed that candidate genes involved in salinity-alkalinity and low temperature adaptation were differentially expressed in tissues of T. yarkandensis. Precisely, mapk1, abcc1, gpx1, gpx4, cat and aqp1 genes participated in the regulation process of salinity-alkalinity adaptation, and elovl4, acaca, fasn, acaa2, acox1 and acox3 genes were involved in fatty acid metabolism and closely associated with low temperature adaptation. On the one hand, it was found that the expression of these genes varied among different tissues, and the important pathways involved in these genes were mapped. Furthermore, we analyzed mapk1 and acox1 genes in depth to obtain the predicted gene structure and important amino acid sites. The transcriptome information in this study will be conducive to provide further understanding for the molecular level research and exploration of the environmental adaptation of T. yarkandensis.

Whole Genome Analysis of SLs Pathway Genes and Functional Characterization of DlSMXL6 in Longan Early Somatic Embryo Development.

Strigolactones (SLs), a new class of plant hormones, are implicated in the regulation of various biological processes. However, the related family members and functions are not identified in longan (Dimocarpus longan Lour.). In this study, 23 genes in the CCD, D27, and SMXL family were identified in the longan genome. The phylogenetic relationships, gene structure, conserved motifs, promoter elements, and transcription factor-binding site predictions were comprehensively analysed. The expression profiles indicated that these genes may play important roles in longan organ development and abiotic stress responses, especially during early somatic embryogenesis (SE). Furthermore, GR24 (synthetic SL analogue) and Tis108 (SL biosynthesis inhibitor) could affect longan early SE by regulating the levels of endogenous IAA (indole-3-acetic acid), JA (jasmonic acid), GA (gibberellin), and ABA (abscisic acid). Overexpression of SMXL6 resulted in inhibition of longan SE by regulating the synthesis of SLs, carotenoids, and IAA levels. This study establishes a foundation for further investigation of SL genes and provides novel insights into their biological functions.

Genome-wide identification and expression analysis of the AhTrx family genes in peanut

Thioredoxins (Trx) are small multifunctional redox proteins that contain thioredoxin conserved domain and active site WCXXC. The Trx family has an important role in multiple processes, including electron transport, seed germination, redox regulation, biotic and abiotic stresses resistance, etc. Although Trx genes have been extensively characterized in some plants, they have not been reported in peanut until now. The identification of AhTrx genes provides potential candidate genes for studying their effects and regulatory mechanisms in peanut (Arachis hypogaea L.) growth and development, especially under aluminium (Al) stress. It is also helpful to further analyze the Al resistance pathway in plants. Seventy AhTrx genes were identified using a genome-wide search method and conservative domain analysis. Then the basic physicochemical properties, phylogenetic relationship, gene structure, chromosomal localization, and promoter prediction were studied by the bioinformatic methods. Furthermore, the expressions of AhTrx genes under different Al treatment times in two peanut cultivars were tested using a real-time quantitative polymerase chain reaction. Seventy AhTrx genes were identified and characterized. Phylogenetic tree analysis showed that all AhTrx members could be classified into 9 groups with different conserved domains. Motif 1 was found to exist in every sequence, with an active site. Furthermore, the gene structures showed that the AhTrx family was complicated and changeable during evolution. The chromosomal localization indicated that the distribution and density of the Trx family on 20 peanut chromosomes were uneven. Predictive promoter analysis indicated that AhTrx proteins might play a role in phytohormones synthesis and stress response. Finally, the expression patterns of the AhTrx genes showed that every gene was differently expressed under Al treatment in different peanut cultivars, some were obvious, others had no significant difference, some were at a high level, while others were at a low level. This study systematically identifies the Trx gene family in peanut, providing some candidates for further study on its effects and regulatory mechanism under Al stress in peanut.

Genome-Wide Identification and Expression Analysis of XTH Gene Family during Flower-Opening Stages in Osmanthus fragrans.

Osmanthus fragrans is an aromatic plant which is widely used in landscaping and garden greening in China. However, the process of flower opening is significantly affected by ambient temperature changes. Cell expansion in petals is the primary factor responsible for flower opening. Xyloglucan endoglycolase/hydrolase (XTH) is a cell-wall-loosening protein involved in cell expansion or cell-wall weakening. Through whole-genome analysis, 38 OfXTH genes were identified in O. fragrans which belong to the four main phylogenetic groups. The gene structure, chromosomal location, synteny relationship, and cis-acting elements prediction and expression patterns were analyzed on a genome-wide scale. The expression patterns showed that most OfXTHs were closely associated with the flower-opening period of O. fragrans. At the early flower-opening stage (S1 and S2), transcriptome and qRT-PCR analysis revealed the expression of OfXTH24, 27, 32, 35, and 36 significantly increased under low ambient temperature (19 °C). It is speculated that the five genes might be involved in the regulation of flower opening by responding to ambient temperature changes. Our results provide solid foundation for the functional analysis of OfXTH genes and help to explore the mechanism of flower opening responding to ambient temperature in O. fragrans.

Genome-Wide Identification of APETALA2/ETHYLENE RESPONSIVE FACTOR Transcription Factors in Cucurbita moschata and Their Involvement in Ethylene Response.

APETALA2/ETHYLENE RESPONSIVE FACTOR (AP2/ERF), a plant-specific transcription factor (TF) family, plays an essential role in the growth and development of plants, and in their response to biotic and abiotic stresses. However, information on AP2/ERF in Cucurbita moschata (pumpkin), an edible and medicinal vegetable used worldwide, is scarce. A total of 212 AP2/ERF genes were identified in the C. moschata genome (CmoAP2/ERFs). Based on phylogenetic analysis, they were divided into four groups—28 AP2s, 92 ERFs, 86 dehydration-responsive element-binding (DREB) factors, and 6 ABI3/VPs (RAV). The 212 AP2/ERF genes were unevenly distributed on the 20 chromosomes of C. moschata. The results of structural analysis showed the absence of introns on 132 CmoAP2/ERFs. Four pairs of tandem duplication and 155 pairs of segmental duplication events were identified, which indicated that segmental duplications might be the main reason for the expansion of the CmoAP2/ERF family. The analysis of cis-regulatory elements (CREs) showed that most of the CmoAP2/ERFs contained hormone response elements (ABREs, EREs) in their promoters, suggesting that AP2/ERFs could contribute to the processes regulated by ethylene and abscisic acid. By comparing the transcriptome of ethephon-treated and control plants, we found that 16 CmoAP2/ERFs were significantly upregulated after ethephon treatment. Furthermore, we determined the expression patterns of these genes at different developmental stages of female and male flowers. This study provides insights into the identification, classification, physicochemical property, phylogenetic analysis, chromosomal location, gene structure, motif identification, and CRE prediction of the AP2/ERF superfamily in C. moschata. Sixteen CmoAP2/ERF genes were identified as ethylene-inducible genes. The results of this study will be valuable for understanding the roles of CmoAP2/ERFs in ethylene response and should provide a foundation for elucidating the function of AP2/ERF TFs in C. moschata.

Enhancing the prediction of protein coding regions in biological sequence via a deep learning framework with hybrid encoding

Protein coding regions prediction is a very important but overlooked subtask for tasks such as prediction of complete gene structure, coding/noncoding RNA. Many machine learning methods have been proposed for this problem, they first encode a biological sequence into numerical values and then feed them into a classifier for final prediction. However, encoding schemes directly influence the classifier's capability to capture coding features and how to choose a proper encoding scheme remains uncertain. Recently, we proposed a protein coding region prediction method in transcript sequences based on a bidirectional recurrent neural network with non-overlapping 3-mer feature, and achieved considerable improvement over existing methods, but there is still much room to improve the performance. First, 3-mer feature that counts the occurrence frequency of trinucleotides in a biological sequence only reflects local sequence order information between the most contiguous nucleotides, which loses almost all the global sequence order information. Second, kmer features of length k larger than three (e.g., hexamer) may also contain useful information. Based on the two points, we here present a deep learning framework with hybrid encoding for protein coding regions prediction in biological sequences, which effectively exploit global sequence order information, non-overlapping gapped kmer (gkm) features and statistical dependencies among coding labels. 3-fold cross-validation tests on human and mouse biological sequence demonstrate that our proposed method significantly outperforms existing state-of-the-art methods.

Foster thy young: enhanced prediction of orphan genes in assembled genomes.

Proteins encoded by newly-emerged genes (‘orphan genes’) share no sequence similarity with proteins in any other species. They provide organisms with a reservoir of genetic elements to quickly respond to changing selection pressures. Here, we systematically assess the ability of five gene prediction pipelines to accurately predict genes in genomes according to phylostratal origin. BRAKER and MAKER are existing, popular ab initio tools that infer gene structures by machine learning. Direct Inference is an evidence-based pipeline we developed to predict gene structures from alignments of RNA-Seq data. The BIND pipeline integrates ab initio predictions of BRAKER and Direct inference; MIND combines Direct Inference and MAKER predictions. We use highly-curated Arabidopsis and yeast annotations as gold-standard benchmarks, and cross-validate in rice. Each pipeline under-predicts orphan genes (as few as 11 percent, under one prediction scenario). Increasing RNA-Seq diversity greatly improves prediction efficacy. The combined methods (BIND and MIND) yield best predictions overall, BIND identifying 68% of annotated orphan genes, 99% of ancient genes, and give the highest sensitivity score regardless dataset in Arabidopsis. We provide a light weight, flexible, reproducible, and well-documented solution to improve gene prediction.

Prediction of the sea buckthorn AQP gene structure and its spatiotemporal expression pattern under drought stress

Prediction of the sea buckthorn AQP gene structure and its spatiotemporal expression pattern under drought stress

Improved recognition of splice sites in A. thaliana by incorporating secondary structure information into sequence-derived features: a computational study.

The online version contains supplementary material available at 10.1007/s13205-021-03036-8.

High-quality reference genome of Fasciola gigantica: Insights into the genomic signatures of transposon-mediated evolution and specific parasitic adaption in tropical regions.

Fasciola gigantica and Fasciola hepatica are causative pathogens of fascioliasis, with the widest latitudinal, longitudinal, and altitudinal distribution; however, among parasites, they have the largest sequenced genomes, hindering genomic research. In the present study, we used various sequencing and assembly technologies to generate a new high-quality Fasciola gigantica reference genome. We improved the integration of gene structure prediction, and identified two independent transposable element expansion events contributing to (1) the speciation between Fasciola and Fasciolopsis during the Cretaceous-Paleogene boundary mass extinction, and (2) the habitat switch to the liver during the Paleocene-Eocene Thermal Maximum, accompanied by gene length increment. Long interspersed element (LINE) duplication contributed to the second transposon-mediated alteration, showing an obvious trend of insertion into gene regions, regardless of strong purifying effect. Gene ontology analysis of genes with long LINE insertions identified membrane-associated and vesicle secretion process proteins, further implicating the functional alteration of the gene network. We identified 852 predicted excretory/secretory proteins and 3300 protein-protein interactions between Fasciola gigantica and its host. Among them, copper/zinc superoxide dismutase genes, with specific gene copy number variations, might play a central role in the phase I detoxification process. Analysis of 559 single-copy orthologs suggested that Fasciola gigantica and Fasciola hepatica diverged at 11.8 Ma near the Middle and Late Miocene Epoch boundary. We identified 98 rapidly evolving gene families, including actin and aquaporin, which might explain the large body size and the parasitic adaptive character resulting in these liver flukes becoming epidemic in tropical and subtropical regions.