Incongruence in the phylogenomics era.

Genome-scale data and the development of novel statistical phylogenetic approaches have greatly aided the reconstruction of a broad sketch of the tree of life and resolved many of its branches. However, incongruence - the inference of conflicting evolutionary histories - remains pervasive in phylogenomic data, hampering our ability to reconstruct and interpret the tree of life. Biological factors, such as incomplete lineage sorting, horizontal gene transfer, hybridization, introgression, recombination and convergent molecular evolution, can lead to gene phylogenies that differ from the species tree. In addition, analytical factors, including stochastic, systematic and treatment errors, can drive incongruence. Here, we review these factors, discuss methodological advances to identify and handle incongruence, and highlight avenues for future research.

Multiple Origins of Bioluminescence in Beetles and Evolution of Luciferase Function.

Bioluminescence in beetles has long fascinated biologists, with diverse applications in biotechnology. To date, however, our understanding of its evolutionary origin and functional variation mechanisms remains poor. To address these questions, we obtained high-quality reference genomes of luminous and nonluminous beetles in 6 Elateroidea families. We then reconstructed a robust phylogenetic relationship for all luminous families and related nonluminous families. Comparative genomic analyses and biochemical functional experiments suggested that gene evolution within Elateroidea played a crucial role in the origin of bioluminescence, with multiple parallel origins observed in the luminous beetle families. While most luciferase-like proteins exhibited a conserved nonluminous amino acid pattern (TLA346 to 348) in the luciferin-binding sites, luciferases in the different luminous beetle families showed divergent luminous patterns at these sites (TSA/CCA/CSA/LVA). Comparisons of the structural and enzymatic properties of ancestral, extant, and site-directed mutant luciferases further reinforced the important role of these sites in the trade-off between acyl-CoA synthetase and luciferase activities. Furthermore, the evolution of bioluminescent color demonstrated a tendency toward hypsochromic shifts and variations among the luminous families. Taken together, our results revealed multiple parallel origins of bioluminescence and functional divergence within the beetle bioluminescent system.

The influence of the number of tree searches on maximum likelihood inference in phylogenomics.

Maximum likelihood (ML) phylogenetic inference is widely used in phylogenomics. As heuristic searches most likely find suboptimal trees, it is recommended to conduct multiple (e.g., ten) tree searches in phylogenetic analyses. However, beyond its positive role, how and to what extent multiple tree searches aid ML phylogenetic inference remains poorly explored. Here, we found that a random starting tree was not as effective as the BioNJ and parsimony starting trees in inferring ML gene tree and that RAxML-NG and PhyML were less sensitive to different starting trees than IQ-TREE. We then examined the effect of the number of tree searches on ML tree inference with IQ-TREE and RAxML-NG, by running 100 tree searches on 19,414 gene alignments from 15 animal, plant, and fungal phylogenomic datasets. We found that the number of tree searches substantially impacted the recovery of the best-of-100 ML gene tree topology among 100 searches for a given ML program. In addition, all of the concatenation-based trees were topologically identical if the number of tree searches was ≥ 10. Quartet-based ASTRAL trees inferred from 1 to 80 tree searches differed topologically from those inferred from 100 tree searches for 6 /15 phylogenomic datasets. Lastly, our simulations showed that gene alignments with lower difficulty scores had a higher chance of finding the best-of-100 gene tree topology and were more likely to yield the correct trees.

OrthoSNAP: A tree splitting and pruning algorithm for retrieving single-copy orthologs from gene family trees.

Molecular evolution studies, such as phylogenomic studies and genome-wide surveys of selection, often rely on gene families of single-copy orthologs (SC-OGs). Large gene families with multiple homologs in 1 or more species-a phenomenon observed among several important families of genes such as transporters and transcription factors-are often ignored because identifying and retrieving SC-OGs nested within them is challenging. To address this issue and increase the number of markers used in molecular evolution studies, we developed OrthoSNAP, a software that uses a phylogenetic framework to simultaneously split gene families into SC-OGs and prune species-specific inparalogs. We term SC-OGs identified by OrthoSNAP as SNAP-OGs because they are identified using a splitting and pruning procedure analogous to snapping branches on a tree. From 415,129 orthologous groups of genes inferred across 7 eukaryotic phylogenomic datasets, we identified 9,821 SC-OGs; using OrthoSNAP on the remaining 405,308 orthologous groups of genes, we identified an additional 10,704 SNAP-OGs. Comparison of SNAP-OGs and SC-OGs revealed that their phylogenetic information content was similar, even in complex datasets that contain a whole-genome duplication, complex patterns of duplication and loss, transcriptome data where each gene typically has multiple transcripts, and contentious branches in the tree of life. OrthoSNAP is useful for increasing the number of markers used in molecular evolution data matrices, a critical step for robustly inferring and exploring the tree of life.

ClipKIT: A multiple sequence alignment trimming software for accurate phylogenomic inference.

Highly divergent sites in multiple sequence alignments (MSAs), which can stem from erroneous inference of homology and saturation of substitutions, are thought to negatively impact phylogenetic inference. Thus, several different trimming strategies have been developed for identifying and removing these sites prior to phylogenetic inference. However, a recent study reported that doing so can worsen inference, underscoring the need for alternative alignment trimming strategies. Here, we introduce ClipKIT, an alignment trimming software that, rather than identifying and removing putatively phylogenetically uninformative sites, instead aims to identify and retain parsimony-informative sites, which are known to be phylogenetically informative. To test the efficacy of ClipKIT, we examined the accuracy and support of phylogenies inferred from 14 different alignment trimming strategies, including those implemented in ClipKIT, across nearly 140,000 alignments from a broad sampling of evolutionary histories. Phylogenies inferred from ClipKIT-trimmed alignments are accurate, robust, and time saving. Furthermore, ClipKIT consistently outperformed other trimming methods across diverse datasets, suggesting that strategies based on identifying and retaining parsimony-informative sites provide a robust framework for alignment trimming.

A comparative analysis of mitochondrial ORFs provides new insights on expansion of mitochondrial genome size in Arcidae.

BACKGROUND: Arcidae, comprising about 260 species of ark shells, is an ecologically and economically important lineage of bivalve mollusks. Interestingly, mitochondrial genomes of several Arcidae species are 2-3 times larger than those of most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. The large mitochondrial genome size is mainly due to expansion of unassigned regions (regions that are functionally unassigned). Previous work on unassigned regions of Arcidae mtDNA genomes has focused on nucleotide-level analyses to observe sequence characteristics, however the origin of expansion remains unclear. RESULTS: We assembled six new mitogenomes and sequenced six transcriptomes of Scapharca broughtonii to identify conserved functional ORFs that are transcribed in unassigned regions. Sixteen lineage-specific ORFs with different copy numbers were identified from seven Arcidae species, and 11 of 16 ORFs were expressed and likely biologically active. Unassigned regions of 32 Arcidae mitogenomes were compared to verify the presence of these novel mitochondrial ORFs and their distribution. Strikingly, multiple structural analyses and functional prediction suggested that these additional mtDNA-encoded proteins have potential functional significance. In addition, our results also revealed that the ORFs have a strong connection to the expansion of Arcidae mitochondrial genomes and their large-scale duplication play an important role in multiple expansion events. We discussed the possible origin of ORFs and hypothesized that these ORFs may originate from duplication of mitochondrial genes. CONCLUSIONS: The presence of lineage-specific mitochondrial ORFs with transcriptional activity and potential functional significance supports novel features for Arcidae mitochondrial genomes. Given our observation and analyses, these ORFs may be products of mitochondrial gene duplication. These findings shed light on the origin and function of novel mitochondrial genes in bivalves and provide new insights into evolution of mitochondrial genome size in metazoans.

Rooting the Animal Tree of Life.

Identifying our most distant animal relatives has emerged as one of the most challenging problems in phylogenetics. This debate has major implications for our understanding of the origin of multicellular animals and of the earliest events in animal evolution, including the origin of the nervous system. Some analyses identify sponges as our most distant animal relatives (Porifera-sister hypothesis), and others identify comb jellies (Ctenophora-sister hypothesis). These analyses vary in many respects, making it difficult to interpret previous tests of these hypotheses. To gain insight into why different studies yield different results, an important next step in the ongoing debate, we systematically test these hypotheses by synthesizing 15 previous phylogenomic studies and performing new standardized analyses under consistent conditions with additional models. We find that Ctenophora-sister is recovered across the full range of examined conditions, and Porifera-sister is recovered in some analyses under narrow conditions when most outgroups are excluded and site-heterogeneous CAT models are used. We additionally find that the number of categories in site-heterogeneous models is sufficient to explain the Porifera-sister results. Furthermore, our cross-validation analyses show CAT models that recover Porifera-sister have hundreds of additional categories and fail to fit significantly better than site-heterogenuous models with far fewer categories. Systematic and standardized testing of diverse phylogenetic models suggests that we should be skeptical of Porifera-sister results both because they are recovered under such narrow conditions and because the models in these conditions fit the data no better than other models that recover Ctenophora-sister.

PhyKIT: a broadly applicable UNIX shell toolkit for processing and analyzing phylogenomic data.

MOTIVATION: Diverse disciplines in biology process and analyze multiple sequence alignments (MSAs) and phylogenetic trees to evaluate their information content, infer evolutionary events and processes and predict gene function. However, automated processing of MSAs and trees remains a challenge due to the lack of a unified toolkit. To fill this gap, we introduce PhyKIT, a toolkit for the UNIX shell environment with 30 functions that process MSAs and trees, including but not limited to estimation of mutation rate, evaluation of sequence composition biases, calculation of the degree of violation of a molecular clock and collapsing bipartitions (internal branches) with low support. RESULTS: To demonstrate the utility of PhyKIT, we detail three use cases: (1) summarizing information content in MSAs and phylogenetic trees for diagnosing potential biases in sequence or tree data; (2) evaluating gene-gene covariation of evolutionary rates to identify functional relationships, including novel ones, among genes and (3) identify lack of resolution events or polytomies in phylogenetic trees, which are suggestive of rapid radiation events or lack of data. We anticipate PhyKIT will be useful for processing, examining and deriving biological meaning from increasingly large phylogenomic datasets. AVAILABILITY AND IMPLEMENTATION: PhyKIT is freely available on GitHub (https://github.com/JLSteenwyk/PhyKIT), PyPi (https://pypi.org/project/phykit/) and the Anaconda Cloud (https://anaconda.org/JLSteenwyk/phykit) under the MIT license with extensive documentation and user tutorials (https://jlsteenwyk.com/PhyKIT). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Mitogenomics reveals phylogenetic relationships of Arcoida (Mollusca, Bivalvia) and multiple independent expansions and contractions in mitochondrial genome size.

Arcoida, comprising about 570 species of blood cockles, is an ecologically and economically important lineage of bivalve molluscs. Current classification of arcoids is largely based on morphology, which shows widespread homoplasy. Despite two recent studies employing multi-locus analyses with broad sampling of Arcoida, evolutionary relationships among major lineages remain controversial. Interestingly, mitochondrial genomes of several ark shell species are 2-3 times larger than those found in most bilaterians, and are among the largest bilaterian mitochondrial genomes reported to date. These results highlight the need of detailed phylogenetic study to explore evolutionary relationships within Arcoida so that the evolution of mitochondrial genome size can be understood. To this end, we sequenced 17 mitochondrial genomes and compared them with publicly available data, including those from other lineages of Arcoida with emphasis on the subclade Arcoidea species. Our phylogenetic analyses indicate that Noetiidae, Cucullaeidae and Glycymerididae are nested within a polyphyletic Arcidae. Moreover, we find multiple independent expansions and potential contractions of mitochondrial genome size, suggesting that the large mitochondrial genome is not a shared ancestral feature in Arcoida. We also examined tandem repeats and inverted repeats in non-coding regions and investigated the presence of such repeats with relation to genome size variation. Our results suggest that tandem repeats might facilitate intraspecific mitochondrial genome size variation, and that inverted repeats, which could be derived from transposons, might be responsible for mitochondrial genome expansions and contractions. We show that mitochondrial genome size in Arcoida is more dynamic than previously understood and provide insights into evolution of mitochondrial genome size variation in metazoans.

Posterior Fusiform and Midfusiform Contribute to Distinct Stages of Facial Expression Processing.

Though the fusiform is well-established as a key node in the face perception network, its role in facial expression processing remains unclear, due to competing models and discrepant findings. To help resolve this debate, we recorded from 17 subjects with intracranial electrodes implanted in face sensitive patches of the fusiform. Multivariate classification analysis showed that facial expression information is represented in fusiform activity and in the same regions that represent identity, though with a smaller effect size. Examination of the spatiotemporal dynamics revealed a functional distinction between posterior fusiform and midfusiform expression coding, with posterior fusiform showing an early peak of facial expression sensitivity at around 180 ms after subjects viewed a face and midfusiform showing a later and extended peak between 230 and 460 ms. These results support the hypothesis that the fusiform plays a role in facial expression perception and highlight a qualitative functional distinction between processing in posterior fusiform and midfusiform, with each contributing to temporally segregated stages of expression perception.

Genomic adaptations to chemosymbiosis in the deep-sea seep-dwelling tubeworm Lamellibrachia luymesi.

BACKGROUND: Symbiotic relationships between microbes and their hosts are widespread and diverse, often providing protection or nutrients, and may be either obligate or facultative. However, the genetic mechanisms allowing organisms to maintain host-symbiont associations at the molecular level are still mostly unknown, and in the case of bacterial-animal associations, most genetic studies have focused on adaptations and mechanisms of the bacterial partner. The gutless tubeworms (Siboglinidae, Annelida) are obligate hosts of chemoautotrophic endosymbionts (except for Osedax which houses heterotrophic Oceanospirillales), which rely on the sulfide-oxidizing symbionts for nutrition and growth. Whereas several siboglinid endosymbiont genomes have been characterized, genomes of hosts and their adaptations to this symbiosis remain unexplored. RESULTS: Here, we present and characterize adaptations of the cold seep-dwelling tubeworm Lamellibrachia luymesi, one of the longest-lived solitary invertebrates. We sequenced the worm's ~ 688-Mb haploid genome with an overall completeness of ~ 95% and discovered that L. luymesi lacks many genes essential in amino acid biosynthesis, obligating them to products provided by symbionts. Interestingly, the host is known to carry hydrogen sulfide to thiotrophic endosymbionts using hemoglobin. We also found an expansion of hemoglobin B1 genes, many of which possess a free cysteine residue which is hypothesized to function in sulfide binding. Contrary to previous analyses, the sulfide binding mediated by zinc ions is not conserved across tubeworms. Thus, the sulfide-binding mechanisms in sibgolinids need to be further explored, and B1 globins might play a more important role than previously thought. Our comparative analyses also suggest the Toll-like receptor pathway may be essential for tolerance/sensitivity to symbionts and pathogens. Several genes related to the worm's unique life history which are known to play important roles in apoptosis, cell proliferation, and aging were also identified. Last, molecular clock analyses based on phylogenomic data suggest modern siboglinid diversity originated in 267 mya (± 70 my) support previous hypotheses indicating a Late Mesozoic or Cenozoic origins of approximately 50-126 mya for vestimentiferans. CONCLUSIONS: Here, we elucidate several specific adaptations along various molecular pathways that link phenome to genome to improve understanding of holobiont evolution. Our findings of adaptation in genomic mechanisms to reducing environments likely extend to other chemosynthetic symbiotic systems.

Conservation of mitochondrial genome arrangements in brittle stars (Echinodermata, Ophiuroidea).

Brittle stars are conspicuous members of benthic ecosystems, fill many ecological niches and are the most speciose of all classes of echinoderms. With high levels of biodiversity, elucidating the evolutionary history of this group is important. Understanding of higher-level relationships within Ophiuroidea has been aided by multilocus nuclear data and DNA barcoding. However, the degree of consistency between mitochondrial and nuclear data within ophiuroids remains unclear and deserves further assessment. In this study, 17 mitochondrial genomes spanning the taxonomic breadth of Ophiuroidea were utilized to explore evolutionary relationships through maximum likelihood analyses, Bayesian inference and comparative assessment of gene order. Our phylogenetic analyses, based on both nucleotide and amino acid residues, support recent findings based on multilocus nuclear data and morphology, in that the brittle star clades Ophintegrida and Euryophiurida were recovered as monophyletic with the latter comprising Euyalida, Ophiuridae and Ophiopyrgidae. Only three different arrangements of the 13 protein coding and 2 ribosomal RNA genes were observed. As expected, tRNA genes were more likely to have undergone rearrangement but the order of all 37 genes was found to be conserved in all sampled Euryalida and Ophiuridae. Both Euryalida and the clade comprised of Ophiuridae and Ophiopyrgidae, each had their own conserved rearrangement of protein coding genes and ribosomal genes, after divergence from their last common ancestor. Euryalida has a rearrangement of the two ribosomal RNA genes, rrnS and rrnL, in contrast to Ophiuridae and Ophiopyrgidae, which had an inversion of the genes nad1, nad2, and cob relative to Ophintegrida. Further, our data support the gene order found in all sampled Euryalida as the most likely ancestral order for all Ophiuroidea.

Decoding and disrupting left midfusiform gyrus activity during word reading.

The nature of the visual representation for words has been fiercely debated for over 150 y. We used direct brain stimulation, pre- and postsurgical behavioral measures, and intracranial electroencephalography to provide support for, and elaborate upon, the visual word form hypothesis. This hypothesis states that activity in the left midfusiform gyrus (lmFG) reflects visually organized information about words and word parts. In patients with electrodes placed directly in their lmFG, we found that disrupting lmFG activity through stimulation, and later surgical resection in one of the patients, led to impaired perception of whole words and letters. Furthermore, using machine-learning methods to analyze the electrophysiological data from these electrodes, we found that information contained in early lmFG activity was consistent with an orthographic similarity space. Finally, the lmFG contributed to at least two distinguishable stages of word processing, an early stage that reflects gist-level visual representation sensitive to orthographic statistics, and a later stage that reflects more precise representation sufficient for the individuation of orthographic word forms. These results provide strong support for the visual word form hypothesis and demonstrate that across time the lmFG is involved in multiple stages of orthographic representation.

Genome sequence of walking catfish (Clarias batrachus) provides insights into terrestrial adaptation.

BACKGROUND: Walking catfish (Clarias batrachus) is a freshwater fish capable of air-breathing and locomotion on land. It usually inhabits various low-oxygen habitats, burrows inside the mudflat, and sometimes "walks" to search for suitable environments during summer. It has evolved accessory air-breathing organs for respiring air and corresponding mechanisms to survive in such challenging environments. Thereby, it serves as a great model for understanding adaptations to terrestrial life. RESULTS: Comparative genomics with channel catfish (Ictalurus punctatus) revealed specific adaptations of C. batrachus in DNA repair, enzyme activator activity, and small GTPase regulator activity. Comparative analysis with 11 non-air-breathing fish species suggested adaptive evolution in gene expression and nitrogenous waste metabolic processes. Further, myoglobin, olfactory receptor related to class A G protein-coupled receptor 1, and sulfotransferase 6b1 genes were found to be expanded in the air-breathing walking catfish genome, with 15, 15, and 12 copies, respectively, compared to non-air-breathing fishes that possess only 1-2 copies of these genes. Additionally, we sequenced and compared the transcriptomes of the gill and the air-breathing organ to characterize the mechanism of aerial respiration involved in elastic fiber formation, oxygen binding and transport, angiogenesis, ion homeostasis and acid-base balance. The hemoglobin genes were expressed dramatically higher in the air-breathing organ than in the gill of walking catfish. CONCLUSIONS: This study provides an important genomic resource for understanding the adaptive mechanisms of walking catfish to terrestrial environments. It is possible that the coupling of enhanced abilities for oxygen storage and oxygen transport through genomic expansion of myoglobin genes and transcriptomic up-regulation of hemoglobin and angiogenesis-related genes are important components of the molecular basis for adaptation of this aquatic species to terrestrial life.

Multi-Connection Pattern Analysis: Decoding the representational content of neural communication.

The lack of multivariate methods for decoding the representational content of interregional neural communication has left it difficult to know what information is represented in distributed brain circuit interactions. Here we present Multi-Connection Pattern Analysis (MCPA), which works by learning mappings between the activity patterns of the populations as a factor of the information being processed. These maps are used to predict the activity from one neural population based on the activity from the other population. Successful MCPA-based decoding indicates the involvement of distributed computational processing and provides a framework for probing the representational structure of the interaction. Simulations demonstrate the efficacy of MCPA in realistic circumstances. In addition, we demonstrate that MCPA can be applied to different signal modalities to evaluate a variety of hypothesis associated with information coding in neural communications. We apply MCPA to fMRI and human intracranial electrophysiological data to provide a proof-of-concept of the utility of this method for decoding individual natural images and faces in functional connectivity data. We further use a MCPA-based representational similarity analysis to illustrate how MCPA may be used to test computational models of information transfer among regions of the visual processing stream. Thus, MCPA can be used to assess the information represented in the coupled activity of interacting neural circuits and probe the underlying principles of information transformation between regions.

Evolution of Sulfur Binding by Hemoglobin in Siboglinidae (Annelida) with Special Reference to Bone-Eating Worms, Osedax.

Most members of Siboglinidae (Annelida) harbor endosymbiotic bacteria that allow them to thrive in extreme environments such as hydrothermal vents, methane seeps, and whale bones. These symbioses are enabled by specialized hemoglobins (Hbs) that are able to bind hydrogen sulfide for transportation to their chemosynthetic endosymbionts. Sulfur-binding capabilities are hypothesized to be due to cysteine residues at key positions in both vascular and coelomic Hbs, especially in the A2 and B2 chains. Members of the genus Osedax, which live on whale bones, do not have chemosynthetic endosymbionts, but instead harbor heterotrophic bacteria capable of breaking down complex organic compounds. Although sulfur-binding capabilities are important in other siboglinids, we questioned whether Osedax retained these cysteine residues and the potential ability to bind hydrogen sulfide. To answer these questions, we used high-throughput DNA sequencing to isolate and analyze Hb sequences from 8 siboglinid lineages. For Osedax mucofloris, we recovered three (A1, A2, and B1) Hb chains, but the B2 chain was not identified. Hb sequences from gene subfamilies A2 and B2 were translated and aligned to determine conservation of cysteine residues at previously identified key positions. Hb linker sequences were also compared to determine similarity between Osedax and siboglinids/sulfur-tolerant annelids. For O. mucofloris, our results found conserved cysteines within the Hb A2 chain. This finding suggests that Hb in O. mucofloris has retained some capacity to bind hydrogen sulfide, likely due to the need to detoxify this chemical compound that is abundantly produced within whale bones.

Mitogenomics reveals phylogeny and repeated motifs in control regions of the deep-sea family Siboglinidae (Annelida).

Deep-sea tubeworms in the annelid family Siboglinidae have drawn considerable interest regarding their ecology and evolutionary biology. As adults, they lack a digestive tract and rely on endosymbionts for nutrition. Moreover, they are important members of chemosynthetic environments including hydrothermal vents, cold seeps, muddy sediments, and whale bones. Evolution and diversification of siboglinids has been associated with host-symbiont relationships and reducing habitats. Despite their importance, the taxonomy and phylogenetics of this clade are debated due to conflicting results. In this study, 10 complete and 2 partial mitochondrial genomes and one transcriptome were sequenced and analyzed to address siboglinid evolution. Notably, repeated nucleotide motifs were found in control regions of these mt genomes, which may explain previous challenges of sequencing siboglinid mt genomes. Phylogenetic analyses of amino acid and nucleotide datasets were conducted in order to infer evolutionary history. Both analyses generally had strong nodal support and suggest Osedax is most closely related to the Vestimentifera+Sclerolinum clade, rather than Frenulata, as recently reported. These results imply Osedax, the only siboglinid lineage with heterotrophic endosymbionts, evolved from a lineage utilizing chemoautotrophic symbionts.

Enhancing neural encoding models for naturalistic perception with a multi-level integration of deep neural networks and cortical networks.

Cognitive neuroscience aims to develop computational models that can accurately predict and explain neural responses to sensory inputs in the cortex. Recent studies attempt to leverage the representation power of deep neural networks (DNNs) to predict the brain response and suggest a correspondence between artificial and biological neural networks in their feature representations. However, typical voxel-wise encoding models tend to rely on specific networks designed for computer vision tasks, leading to suboptimal brain-wide correspondence during cognitive tasks. To address this challenge, this work proposes a novel approach that upgrades voxel-wise encoding models through multi-level integration of features from DNNs and information from brain networks. Our approach combines DNN feature-level ensemble learning and brain atlas-level model integration, resulting in significant improvements in predicting whole-brain neural activity during naturalistic video perception. Furthermore, this multi-level integration framework enables a deeper understanding of the brain's neural representation mechanism, accurately predicting the neural response to complex visual concepts. We demonstrate that neural encoding models can be optimized by leveraging a framework that integrates both data-driven approaches and theoretical insights into the functional structure of the cortical networks.

Quantification of Hypsarrhythmia in Infantile Spasmatic EEG: A Large Cohort Study.

Infantile spasms (IS) is a neurological disorder causing mental and/or developmental retardation in many infants. Hypsarrhythmia is a typical symptom in the electroencephalography (EEG) signals with IS. Long-term EEG/video monitoring is most frequently employed in clinical practice for IS diagnosis, from which manual screening of hypsarrhythmia is time consuming and lack of sufficient reliability. This study aims to identify potential biomarkers for automatic IS diagnosis by quantitative analysis of the EEG signals. A large cohort of 101 IS patients and 155 healthy controls (HC) were involved. Typical hypsarrhythmia and non-hypsarrhythmia EEG signals were annotated, and normal EEG were randomly picked from the HC. Root mean square (RMS), teager energy (TE), mean frequency, sample entropy (SamEn), multi-channel SamEn, multi-scale SamEn, and nonlinear correlation coefficient were computed in each sub-band of the three EEG signals, and then compared using either a one-way ANOVA or a Kruskal-Wallis test (based on their distribution) and the receiver operating characteristic (ROC) curves. The effects of infant age on these features were also investigated. For most of the employed features, significant ( ) differences were observed between hypsarrhythmia EEG and non-hypsarrhythmia EEG or HC, which seem to increase with increased infant age. RMS and TE produce the best classification in the delta and theta bands, while entropy features yields the best performance in the gamma band. Our study suggests RMS and TE (delta and theta bands) and entropy features (gamma band) to be promising biomarkers for automatic detection of hypsarrhythmia in long-term EEG monitoring. The findings of our study indicate the feasibility of automated IS diagnosis using artificial intelligence.