Finding Needles in the Haystack: Strategies for Uncovering Noncoding Regulatory Variants.

Despite accumulating evidence implicating noncoding variants in human diseases, unraveling their functionality remains a significant challenge. Systematic annotations of the regulatory landscape and the growth of sequence variant data sets have fueled the development of tools and methods to identify causal noncoding variants and evaluate their regulatory effects. Here, we review the latest advances in the field and discuss potential future research avenues to gain a more in-depth understanding of noncoding regulatory variants.

Author Correction: The Apostasia genome and the evolution of orchids.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Single-molecule optofluidic microsensor with interface whispering gallery modes.

Label-free sensors are highly desirable for biological analysis and early-stage disease diagnosis. Optical evanescent sensors have shown extraordinary ability in label-free detection, but their potentials have not been fully exploited because of the weak evanescent field tails at the sensing surfaces. Here, we report an ultrasensitive optofluidic biosensor with interface whispering gallery modes in a microbubble cavity. The interface modes feature both the peak of electromagnetic-field intensity at the sensing surface and high-Q factors even in a small-sized cavity, enabling a detection limit as low as 0.3 pg/cm2 The sample consumption can be pushed down to 10 pL due to the intrinsically integrated microfluidic channel. Furthermore, detection of single DNA with 8 kDa molecular weight is realized by the plasmonic-enhanced interface mode.

PbABCG1 and PbABCG2 transporters are required for the emission of floral monoterpenes in Phalaenopsis bellina.

Terrestrial plants emit volatiles into the atmosphere to attract both pollinators and the enemies of herbivores, for defense. Phalaenopsis bellina is a scented orchid species in which the main scent components are monoterpenes, including linalool and geraniol, and their derivatives. Here, we investigated whether ABC transporters are involved in floral scent emission. We carried out whole-genome identification of ABC transporter-related genes using four floral transcriptomics libraries of P. bellina. We identified 86 ABC subfamily G genes related to terpenoid transport. After comparing the gene expression patterns of P. bellina with that of Phalaenopsis aphrodite subsp. formosana, a scentless species, followed by gene-to-gene correlation analysis, PbABCG1 and PbABCG2 were selected. The temporal expression of both PbABCG1 and PbABCG2 was highly correlated with that of the key enzyme PbGDPS and the major transcription factor PbbHLH4 in monoterpene biosynthesis, with optimal expression on day 5 post-anthesis. Spatial gene expression analysis showed that PbABCG1 was highly expressed in sepals, whereas PbABCG2 was expressed in the lip. Subcellular localization with a GFP fusion protein revealed that both PbABCG1 and PbABCG2 are cytoplasmic membrane proteins. Co-downregulation of PbABCG1 and PbABCG2 using both double-strand RNA interference and tobacco rattle virus-based gene silencing led to a significant decrease in monoterpene emission, accompanied by an increase in the internal monoterpene pools. Furthermore, ectopic expression of PbABCG1 and PbABCG2 in an ABC16- mutant yeast strain rescued its tolerance to geraniol. Altogether, our results indicate that PbABCG1 and PbABCG2 play substantial roles in monoterpene transport/emission in P. bellina floral scent.

Investigating Brain-Gut Microbiota Dynamics and Inflammatory Processes in an Autistic-like Rat Model Using MRI Biomarkers During Childhood and Adolescence.

Autism spectrum disorder (ASD) is characterized by social interaction deficits and repetitive behaviors. Recent research has linked that gut dysbiosis may contribute to ASD-like behaviors. However, the exact developmental time point at which gut microbiota alterations affect brain function and behavior in patients with ASD remains unclear. We hypothesized that ASD-related brain microstructural changes and gut dysbiosis induce metabolic dysregulation and proinflammatory responses, which collectively contribute to the social behavioral deficits observed in early childhood. We used an autistic-like rat model that was generated via prenatal valproic acid exposure. We analyzed brain microstructural changes using diffusion tensor imaging (DTI) and examined microbiota, blood, and fecal samples for inflammation biomarkers. The ASD model rats exhibited significant brain microstructural changes in the anterior cingulate cortex, hippocampus, striatum, and thalamus; reduced microbiota diversity (Prevotellaceae and Peptostreptococcaceae); and altered metabolic signatures. The shift in microbiota diversity and density observed at postnatal day (PND) 35, which is a critical developmental period, underscored the importance of early ASD interventions. We identified a unique metabolic signature in the ASD model, with elevated formate and reduced acetate and butyrate levels, indicating a dysregulation in short-chain fatty acid (SCFA) metabolism. Furthermore, increased astrocytic and microglial activation and elevated proinflammatory cytokines-interleukin-1 beta (IL-1ß), interleukin-6 (IL-6), interferon-gamma (IFN-γ), and tumor necrosis factor-alpha (TNF-α)-were observed, indicating immune dysregulation. This study provided insights into the complex interplay between the brain and the gut, and indicated DTI metrics as potential imaging-based biomarkers in ASD, thus emphasizing the need for early childhood interventions.

Testing plastomes and nuclear ribosomal DNA sequences as the next-generation DNA barcodes for species identification and phylogenetic analysis in Acer.

BACKGROUND: Acer is a taxonomically intractable and speciose genus that contains over 150 species. It is challenging to distinguish Acer species only by morphological method due to their abundant variations. Plastome and nuclear ribosomal DNA (nrDNA) sequences are recommended as powerful next-generation DNA barcodes for species discrimination. However, their efficacies were still poorly studied. The current study will evaluate the application of plastome and nrDNA in species identification and perform phylogenetic analyses for Acer. RESULT: Based on a collection of 83 individuals representing 55 species (c. 55% of Chinese species) from 13 sections, our barcoding analyses demonstrated that plastomes exhibited the highest (90.47%) species discriminatory power among all plastid DNA markers, such as the standard plastid barcodes matK + rbcL + trnH-psbA (61.90%) and ycf1 (76.19%). And the nrDNA (80.95%) revealed higher species resolution than ITS (71.43%). Acer plastomes show abundant interspecific variations, however, species identification failure may be due to the incomplete lineage sorting (ILS) and chloroplast capture resulting from hybridization. We found that the usage of nrDNA contributed to identifying those species that were unidentified by plastomes, implying its capability to some extent to mitigate the impact of hybridization and ILS on species discrimination. However, combining plastome and nrDNA is not recommended given the cytonuclear conflict caused by potential hybridization. Our phylogenetic analysis covering 19 sections (95% sections of Acer) and 128 species (over 80% species of this genus) revealed pervasive inter- and intra-section cytonuclear discordances, hinting that hybridization has played an important role in the evolution of Acer. CONCLUSION: Plastomes and nrDNA can significantly improve the species resolution in Acer. Our phylogenetic analysis uncovered the scope and depth of cytonuclear conflict in Acer, providing important insights into its evolution.

A 3D structural SARS-CoV-2-human interactome to explore genetic and drug perturbations.

Emergence of new viral agents is driven by evolution of interactions between viral proteins and host targets. For instance, increased infectivity of SARS-CoV-2 compared to SARS-CoV-1 arose in part through rapid evolution along the interface between the spike protein and its human receptor ACE2, leading to increased binding affinity. To facilitate broader exploration of how pathogen-host interactions might impact transmission and virulence in the ongoing COVID-19 pandemic, we performed state-of-the-art interface prediction followed by molecular docking to construct a three-dimensional structural interactome between SARS-CoV-2 and human. We additionally carried out downstream meta-analyses to investigate enrichment of sequence divergence between SARS-CoV-1 and SARS-CoV-2 or human population variants along viral-human protein-interaction interfaces, predict changes in binding affinity by these mutations/variants and further prioritize drug repurposing candidates predicted to competitively bind human targets. We believe this resource ( http://3D-SARS2.yulab.org ) will aid in development and testing of informed hypotheses for SARS-CoV-2 etiology and treatments.

miR-194 Up-Regulates Cholesterol 7 Alpha-Hydroxylase Expression via ß-Catenin Signaling and Aggravates Cholestatic Liver Diseases.

miR-194 is abundantly expressed in hepatocytes, and its depletion increases hepatic resistance to acetaminophen-induced acute injuries. In this study, the biological role of miR-194 in cholestatic liver injury was investigated by using miR-194/miR-192 cluster liver-specific knockout (LKO) mice, in which no liver injuries or metabolic disorders were predisposed. Bile duct ligation (BDL) and 1-naphthyl isothiocyanate (ANIT) were applied to LKO and matched control wild-type (WT) mice to induce hepatic cholestasis. Periportal liver damage, mortality rate, and liver injury biomarkers in LKO mice were significantly less than in WT mice after BDL and ANIT injection. Intrahepatic bile acid level was significantly lower in the LKO liver within 48 hours of BDL- and ANIT-induced cholestasis compared with WT. Western blot analysis showed that ß-catenin (CTNNB1) signaling and genes involved in cellular proliferation were activated in BDL- and ANIT-treated mice. The expression levels of cholesterol 7 alpha-hydroxylase (CYP7A1), pivotal in bile synthesis, and its upstream regulator hepatocyte nuclear factor 4α were reduced in primary LKO hepatocytes and liver tissues compared with WT. The knockdown of miR-194 using miRNA inhibitors reduced CYP7A1 expression in WT hepatocytes. In contrast, the knockdown of CTNNB1 and overexpression of miR-194, but not miR-192, in LKO hepatocytes and AML12 cells increased CYP7A1 expression. In conclusion, the results suggest that the loss of miR-194 ameliorates cholestatic liver injury and may suppress CYP7A1 expression via activation of CTNNB1 signaling.

Feedback insensitivity in a self-chaotic microcavity laser.

Insensitivity to external optical feedback is experimentally demonstrated in a self-chaotic deformed square microcavity laser for the first time, to the best of our knowledge. Both the optical and radio frequency (RF) spectra of the microlaser remain unaffected for external optical feedback with feedback strength as high as 9.9 dB. In addition, the autocorrelation function curve exhibits no time-delayed peaks. The insensitivity makes the self-chaotic microcavity laser promising for applications in feedback-insensitive optical sources.

Cigarette smoke-induced dysbiosis: comparative analysis of lung and intestinal microbiomes in COPD mice and patients.

BACKGROUND: The impact of cigarette smoke (CS) on lung diseases and the role of microbiome dysbiosis in chronic obstructive pulmonary disease (COPD) have been previously reported; however, the relationships remain unclear. METHODS: Our research examined the effects of 20-week cigarette smoke (CS) exposure on the lung and intestinal microbiomes in C57BL/6JNarl mice, alongside a comparison with COPD patients' intestinal microbiome data from a public dataset. RESULTS: The study found that CS exposure significantly decreased forced vital capacity (FVC), thickened airway walls, and induced emphysema. Increased lung damage was observed along with higher lung keratinocyte chemoattractant (KC) levels by CS exposure. Lung microbiome analysis revealed a rise in Actinobacteriota, while intestinal microbiome showed significant diversity changes, indicating dysbiosis. Principal coordinate analysis highlighted distinct intestinal microbiome compositions between control and CS-exposed groups. In the intestinal microbiome, notable decreases in Patescibacteria, Campilobacterota, Defferibacterota, Actinobacteriota, and Desulfobacterota were observed. We also identified correlations between lung function and dysbiosis in both lung and intestinal microbiomes. Lung interleukins, interferon-É£, KC, and 8-isoprostane levels were linked to lung microbiome dysbiosis. Notably, dysbiosis patterns in CS-exposed mice were similar to those in COPD patients, particularly of Global Initiative for Chronic Obstructive Lung Disease (GOLD) stage 4 patients. This suggests a systemic impact of CS exposure. CONCLUSION: In summary, CS exposure induces significant dysbiosis in lung and intestinal microbiomes, correlating with lung function decline and injury. These results align with changes in COPD patients, underscoring the important role of microbiome in smoke-related lung diseases.

In Situ-Generated Formamidine as a Carbon/Nitrogen Source for Enaminone Formation: One-Pot Synthesis of Functionalized 4-Acyl-1,2,3-triazoles.

N,N-Dimethylformamide was reacted with hexamethyldisilazane to generate an N,N-dimethylformimidamide intermediate; thereafter, a reaction with acetophenones/ß-diketones was induced to form enaminones. The one-pot synthetic protocol described in this paper can be applied to synthesize 1,4-disubstituted 1,2,3-triazoles and 1,4,5-trisubstituted 1,2,3-triazoles, in which organic azides are used as substrates under optimized conditions. Furthermore, this protocol uses readily available materials, is nearly free of solvent, can be applied to gram-scale operations, and leads to the formation of structurally diverse products with favorable yields.

Monitoring the Biochemical Activity of Single Anammox Granules with Microbarometers.

Granule-based anaerobic ammonium oxidation (Anammox) is a promising biotechnology for wastewater treatments with extraordinary performance in nitrogen removal. However, traditional analytical methods often delivered an average activity of a bulk sample consisting of millions and even billions of Anammox granules with distinct sizes and components. Here, we developed a novel technique to monitor the biochemical activity of individual Anammox granules in real-time by recording the production rate of nitrogen gas with a microbarometer in a sealed chamber containing only one granule. It was found that the specific activity of a single Anammox granule not only varied by tens of folds among different individuals with similar sizes (activity heterogeneity) but also revealed significant breath-like dynamics over time (temporal fluctuation). Statistical analysis on tens of individuals further revealed two subpopulations with distinct color and specific activity, which were subsequently attributed to the different expression levels of heme c content and hydrazine dehydrogenase activity. This study not only provides a general methodology for various kinds of gas-producing microbial processes but also establishes a bottom-up strategy for exploring the structural-activity relationship at a single sludge granule level, with implications for developing a better Anammox process.

Potential Role of the Anammoxosome in the Adaptation of Anammox Bacteria to Salinity Stress.

The underlying adaptative mechanisms of anammox bacteria to salt stress are still unclear. The potential role of the anammoxosome in modulating material and energy metabolism in response to salinity stress was investigated in this study. The results showed that anammox bacteria increased membrane fluidity and decreased mechanical properties by shortening the ladderane fatty acid chain length of anammoxosome in response to salinity shock, which led to the breakdown of the proton motive force driving ATP synthesis and retarded energy metabolism activity. Afterward, the fatty acid chain length and membrane properties were recovered to enhance the energy metabolic activity. The relative transmission electron microscopy (TEM) area proportion of anammoxosome decreased from 55.9 to 38.9% under salinity stress. The 3D imaging of the anammox bacteria based on Synchrotron soft X-ray tomography showed that the reduction in the relative volume proportion of the anammoxosome and the concave surfaces was induced by salinity stress, which led to the lower energy expenditure of the material transportation and provided more binding sites for enzymes. Therefore, anammox bacteria can modulate nitrogen and energy metabolism by changing the membrane properties and morphology of the anammoxosome in response to salinity stress. This study broadens the response mechanism of anammox bacteria to salinity stress.

Exploring potential relationships between acoustic indices and ecosystem functions: a test on insect herbivory.

Biodiversity loss is a global concern. Current technological advances allow the development of novel tools that can monitor biodiversity remotely with minimal disturbance. One example is passive acoustic monitoring (PAM), which involves recording the soundscape of an area using autonomous recording units, and processing these data using acoustic indices, for example, to estimate the diversity of various vocal animal groups. We explored the hypothesis that data obtained through PAM could also be used to study ecosystem functions. Specifically, we investigated the potential relationship between seven commonly used acoustic indices and insect leaf herbivory, measured as total leaf damage and as the damage from three major insect feeding guilds. Herbivory was quantified on seedlings in 13 plots in four subtropical forests in south China, and acoustic data, representing insect acoustic complexity, were obtained by recording the evening soundscapes in those same locations. Herbivory levels correlated positively with the acoustic entropy index, commonly reported as one of the best-performing indices, whose high values indicate higher acoustic complexity, likely due to greater insect diversity. Relationships for specific feeding guilds were moderately stronger for chewers, indicating that the acoustic indices capture some insect groups more than others (e.g., chewers include soniferous taxa such as crickets, whereas miners are mostly silent). Our findings suggest that the use of PAM to monitor ecosystem functions deserves to be explored further, as this is a research field with unexplored potential. Well-designed targeted studies could help us better understand how to best use novel technologies to monitor ecosystem functions.

Discovering clinical drug-drug interactions with known pharmacokinetics mechanisms using spontaneous reporting systems and electronic health records.

OBJECTIVE: Although the mechanisms behind pharmacokinetic (PK) drug-drug interactions (DDIs) are well-documented, bridging the gap between this knowledge and clinical evidence of DDIs, especially for serious adverse drug reactions (SADRs), remains challenging. While leveraging the FDA Adverse Event Reporting System (FAERS) database along with disproportionality analysis tends to detect a vast number of DDI signals, this abundance complicates further investigation, such as validation through clinical trials. Our study proposed a framework to efficiently prioritize these signals and assessed their reliability using multi-source Electronic Health Records (EHR) to identify top candidates for further investigation. METHODS: We analyzed FAERS data spanning from January 2004 to March 2023, employing four established disproportionality methods: Proportional Reporting Ratio (PRR), Reporting Odds Ratio (ROR), Multi-item Gamma Poisson Shrinker (MGPS), and Bayesian Confidence Propagating Neural Network (BCPNN). Building upon these models, we developed four ranking models to prioritize DDI-SADR signals and cross-referenced signals with DrugBank. To validate the top-ranked signals, we employed longitudinal EHRs from Vanderbilt University Medical Center and the All of Us research program. The performance of each model was assessed by counting how many of the top-ranked signals were confirmed by EHRs and calculating the average ranking of these confirmed signals. RESULTS: Out of 189 DDI-SADR signals identified by all four disproportionality methods, only two were documented in the DrugBank database. By prioritizing the top 20 signals as determined by each of the four disproportionality methods and our four ranking models, 58 unique DDI-SADR signals were selected for EHR validations. Of these, five signals were confirmed. The ranking model, which integrated the MGPS and BCPNN, demonstrated superior performance by assigning the highest priority to those five EHR-confirmed signals. CONCLUSION: The fusion of disproportionality analysis with ranking models, validated through multi-source EHRs, presents a groundbreaking approach to pharmacovigilance. Our study's confirmation of five significant DDI-SADRs, previously unrecorded in the DrugBank database, highlights the essential role of advanced data analysis techniques in identifying ADRs.

Evolutionary history of the angiosperm flora of China.

High species diversity may result from recent rapid speciation in a 'cradle' and/or the gradual accumulation and preservation of species over time in a 'museum'. China harbours nearly 10% of angiosperm species worldwide and has long been considered as both a museum, owing to the presence of many species with hypothesized ancient origins, and a cradle, as many lineages have originated as recent topographic changes and climatic shifts-such as the formation of the Qinghai-Tibetan Plateau and the development of the monsoon-provided new habitats that promoted remarkable radiation. However, no detailed phylogenetic study has addressed when and how the major components of the Chinese angiosperm flora assembled to form the present-day vegetation. Here we investigate the spatio-temporal divergence patterns of the Chinese flora using a dated phylogeny of 92% of the angiosperm genera for the region, a nearly complete species-level tree comprising 26,978 species and detailed spatial distribution data. We found that 66% of the angiosperm genera in China did not originate until early in the Miocene epoch (23 million years ago (Mya)). The flora of eastern China bears a signature of older divergence (mean divergence times of 22.04-25.39 Mya), phylogenetic overdispersion (spatial co-occurrence of distant relatives) and higher phylogenetic diversity. In western China, the flora shows more recent divergence (mean divergence times of 15.29-18.86 Mya), pronounced phylogenetic clustering (co-occurrence of close relatives) and lower phylogenetic diversity. Analyses of species-level phylogenetic diversity using simulated branch lengths yielded results similar to genus-level patterns. Our analyses indicate that eastern China represents a floristic museum, and western China an evolutionary cradle, for herbaceous genera; eastern China has served as both a museum and a cradle for woody genera. These results identify areas of high species richness and phylogenetic diversity, and provide a foundation on which to build conservation efforts in China.

Gigantol protects retinal pigment epithelial cells against high glucose-induced apoptosis, oxidative stress and inflammation by inhibiting MTDH-mediated NF-kB signaling pathway.

OBJECTIVE: As a frequent complication of diabetes mellitus (DM), diabetic retinopathy (DR) is now one of the major causes of blindness. Recent reports have shown that retinal pigment epithelial cell (RPEC) damage plays an essential part in DR development and progression. This work intended to explore the potential effects of Gigantol on high glucose (HG)-stimulated RPEC damage and identify potential mechanisms. METHODS: Cell viability, cell damage, and cell apoptosis were evaluated by CCK-8, lactate dehydrogenase (LDH) and flow cytometry assays. The levels of oxidative stress biomarkers and pro-inflammatory cytokines were assessed using corresponding commercial kits and ELISA. Additionally, the levels of MTDH and NF-kB signaling pathway-related proteins were detected by western blotting. RESULTS: Gigantol dose-dependently enhanced cell viability and decreased apoptosis in HG-challenged ARPE-19 cells. Also, Gigantol notably relieved oxidative stress and inflammatory responses in ARPE-19 cells under HG conditions. Gigantol dose-dependently suppressed MTDH expression. In addition, MTDH restoration partially counteracted the protective effects of Gigantol on ARPE-19 cells subject to HG treatment. Mechanically, Gigantol inactivated the NF-kB signaling pathway, which was partly restored after MTDH overexpression. CONCLUSION: Our findings suggested that Gigantol protected against HG-induced RPEC damage by inactivating the NF-kB signaling via MTDH inhibition, offering a potent therapeutic drug for DR treatment.

Comparison of three-dimensional heads-up system versus traditional microscopic system in medical education for vitreoretinal surgeries: a prospective study.

BACKGROUND: To compare the value and efficiency of the three-dimensional (3D) heads-up surgical system and traditional microscopic (TM) system in teaching and learning vitreoretinal surgeries. METHODS: Twenty ophthalmologists and scrub nurses were recruited as teachers, and 45 junior ophthalmology residents and trainee doctors, trainee nurses, and medical students were recruited as observers. Each teacher and observer were assigned to both a 3D-assisted and TM-assisted vitreoretinal surgery and then asked to complete satisfaction questionnaires for both surgical systems at the end of each surgery. RESULTS: The 3D heads-up surgical system was rated significantly higher in most of the subscales and overall satisfaction score by both teachers and observers (P < 0.05). However, ratings for instrument adjustment were significantly higher in the TM group compared to the 3D group for junior ophthalmology residents and trainee doctors (6.1 ± 1.7 vs. 8.8 ± 1.1, P < 0.001). CONCLUSIONS: The 3D heads-up surgical system has great didactical value in the medical education of vitreoretinal surgeries, but it is important to consider the specific needs of different learners when choosing between the two systems. TRIAL REGISTRATION: Not applicable.

Phase-Transition Microcavity Laser.

Liquid-crystal microcavity lasers have attracted considerable attention because of their extraordinary tunability and sensitive response to external stimuli, and because they operate generally within a specific phase. Here, we demonstrate a liquid-crystal microcavity laser operated in the phase transition in which the reorientation of liquid-crystal molecules occurs from aligned to disordered states. A significant wavelength shift of the microlaser is observed, resulting from the dramatic changes in the refractive index of liquid-crystal microdroplets during the phase transition. This phase-transition microcavity laser is then exploited for sensitive thermal sensing, enabling a two-order-of-magnitude enhancement in sensitivity compared with the nematic-phase microlaser operated far from the transition point. Experimentally, we demonstrate an exceptional sensitivity of -40 nm/K and an ultrahigh resolution of 320 µK. The phase-transition microcavity laser features compactness, softness, and tunability, showing great potential for high-performance sensors, optical modulators, and soft matter photonics.

Mitochondrial Genome Characteristics and Phylogenetic Analysis of Fulmekiola serrata (Kobus) (Thysanoptera: Thripidae).

Sugarcane thrips, Fulmekiola serrata (Kobus) (Thysanoptera: Thripidae), is a common foliar pest that infests sugarcane and is found throughout tropical and subtropical countries. In this study, we obtained and analyzed the complete mitochondrial genome of F. serrata for the first time and explored the phylogenetic relationships of the higher-order elements of Thysanoptera members at the mitochondrial level. The complete mitochondrial genome of F. serrata is 16,596 bp in length and includes 13 protein-coding genes (PCGs), 22 transfer RNA (tRNA) genes, 2 ribosomal RNA (rRNA) genes, and 1 noncoding control region. A+T accounted for 75% of the total bases in the mitochondrial genome of F. serrata, revealing an obvious AT bias. Among the 13 PCGs, except for nad5, which had a start codon of TTG, the remaining genes had ATNs typical of insects (ATA, ATT, ATC, and ATG); nad1, nad2, nad3, and atp8 had incomplete termination codons of TA or T. The remaining nine PCGs were complete with the termination codon TAA. Of the 22 tRNA secondary structures, all were typical cloverleaf secondary structures except for trnS1, which was missing the DHU arm. Compared with the hypothetical ancestral gene arrangement of arthropods, F. serrata presented extensive gene rearrangement, with 23 translocated genes, 8 inverted genes, and 5 shuffled genes. Both maximum likelihood (ML) and Bayesian inference (BI) phylogenetic trees resulted in similar topologies: ((Thripidae + (Stenurothripidae + Aeolothripidae)) + Phlaeothripidae), with Thripidae, Aeolothripidae and Phlaeothripidae being monophyletic groups, whereas F. serrata is closely related to Thrips palmi, and the two are sister groups.