Global Phylogenomic Radiation of Bemisia Whiteflies and Allied Aleyrodidae (Hemiptera: Sternorrhyncha), Native Rust Quartet Coalescent Species-Tree Inference, Symbiotic and Cytochrome P450 Selection, and Cross-Regional vs. Multi-Locus Demographic Skylines

PREPRINT-23 (v1.0.0 Empirical Preprint) • 100% Native Rust bio-suite-supervisor Across 544 Verified NCBI Records (L = 90,311 Aligned Columns) We reconstruct the Triassic-to-Anthropocene global phylogenomic radiation, symbiotic and metabolic enzyme evolution, and multi-compartment demographic history of Bemisia whiteflies and allied Aleyrodidae (Hemiptera: Sternorrhyncha) across 544 active verified NCBI sequence records (plus 6 autonomously quarantined GenBank UNVERIFIED anomalies) comprising a 64-taxon, 18-locus supermatrix (L = 90,311 aligned columns), 48 nuclear functional coding sequences, and 12 single-species regional and multi-locus Bemisia tabaci cohorts (12 x N = 36 = 432 isolates). Executed end-to-end via autonomous Native Rust bio-suite-supervisor optimization across 32 candidate RuNIE chains (total runtime 503.59 s; B = 1,000 UFBoot; 25,000 MCMC generations with 20% burn-in), RY3 purine/pyrimidine recoding under GTR+F+R4 reduces Sternorrhynchan compositional heterogeneity by 88.0% (RCFV = 0.0063 vs. 0.0525, ln L = -631,907.97), dating the Sternorrhyncha root to 234.6 Ma, Crown Aleyrodidae to 129.4 Ma, the Crown Bemisia tabaci cryptic species complex to 20.6 Ma, and the MED-MEAM1-IO crown divergence to 2.86 Ma. Native Rust quartet multispecies coalescent inference (runie-quartet-coalescent across 18 gene trees and 11,436,768 quartets in 89.77 ms) recovers 8,310,800 satisfied quartets (normalized Q = 0.7267, mean LPP = 0.8602) while ABBA-BABA tests quantify asymmetric MED-MEAM1 and SSA1-IO admixture. Sixty-one-state Goldman-Yang codon models identify strong purifying constraint on Orco chemosensation (omega_0 = 0.208) and Portiera bacteriocyte symbiosis alongside diversifying selection on salivary BtFer1/ALP hydrolases (omega_2 = 1.00) and CYP6CM1 substrate-recognition pockets (omega_2 = 1.00). Finally, 12 coupled Bayesian coalescent skylines contrast Old World savanna and agrarian expansions against New World founder bottlenecks and reveal sharp locus-specific effective population compressions at nuclear CYP6CM1 and ChE coding loci relative to neutral mitochondrial COX1 (1Ne) and Katoh E-INS-i-aligned nuclear ITS1-5.8S (4Ne) reservoirs. Autonomous 100% Native Rust bio-suite-supervisor synthesis across 544 verified NCBI records (L = 90,311 aligned columns, 32 RuNIE chains, 11.44M quartets, 300k MCMC generations in 503.59 s) resolves the 234.6 Ma-to-present Aleyrodidae and Bemisia radiation, CYP6CM1 selection, and 12 regional vs. multi-locus skylines. Canonical Preprint URL: https://crvscience.com/labs/preprints/bemisia-global-radiation-preprint Master (Concept) Zenodo DOI: 10.5281/zenodo.24093136 Open-Source Software Repository (GNU GPL-3.0): https://github.com/bpwhite/crvscience-labs