Ecology

Boutaiba, S., Hacene, H., Bidle, K. A. & Maupin-Furlow, J. A. Microbial diversity of the hypersaline Sidi Ameur and Himalatt Salt Lakes of the Algerian Sahara. J. Arid Environ. 75, 909–916. https://doi.org/10.1016/j.jaridenv.2011.04.010 (2011).ADS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Ventosa, A. Unusual micro-organisms from unusual habitats: hypersaline environments. Symposia Society for General Microbiology (2006).Fukuchi, S., Yoshimune, K., Wakayama, M., Moriguchi, M. & Nishikawa, K. Unique amino acid composition of proteins in halophilic bacteria. J. Mol. Biol. 327, 347–357 (2003).CAS 
Article 

Google Scholar 
Pillai, S. D., Nakatsu, C. H., Miller, R. V. & Yates, M. V. Manual of environmental microbiology. Life High-Salinity Environ. https://doi.org/10.1128/9781555818821 (2015).Article 

Google Scholar 
Poli, A. et al. Microbial diversity in extreme marine habitats and their biomolecules. Microorganisms 5, 25. https://doi.org/10.3390/microorganisms5020025 (2017).CAS 
Article 
PubMed Central 

Google Scholar 
Azpiazu-Muniozguren, M., Martinez-Ballesteros, I., Gamboa, J., Seoane, S. & Bikandi, J. Altererythrobacter muriae sp. nov., isolated from hypersaline Aana Salt Valley spring water, a continental thalassohaline-type solar saltern. Int. J. Syst. Evol. Microbiol. 71, 3 (2021).
Google Scholar 
Zhang, J. et al. Bacterial diversity in Bohai Bay Solar Saltworks, China. Curr. Microbiol. 72, 55–63 (2016).CAS 
Article 

Google Scholar 
Highfield, A., Ward, A., Pipe, R. & Schroeder, D. C. Molecular and phylogenetic analysis reveals new diversity of Dunaliella salina from hypersaline environments. J. Mar. Biol. Assoc. UK 101, 27–37. https://doi.org/10.1017/s0025315420001319 (2021).CAS 
Article 

Google Scholar 
Cycil, L. M. et al. Metagenomic insights into the diversity of halophilic microorganisms indigenous to the Karak Salt Mine, Pakistan. Front. Microbiol. https://doi.org/10.3389/fmicb.2020.01567 (2020).Article 
PubMed 
PubMed Central 

Google Scholar 
Jacob, J. H., Hussein, E. I., Shakhatreh, M. A. K. & Cornelison, C. T. Microbial community analysis of the hypersaline water of the Dead Sea using high-throughput amplicon sequencing. Microbiol. Open 6, e00500. https://doi.org/10.1002/mbo3.500 (2017).CAS 
Article 

Google Scholar 
Ben Abdallah, M. et al. Abundance and diversity of prokaryotes in ephemeral hypersaline lake Chott El Jerid using Illumina Miseq sequencing, DGGE and qPCR assay. Extremophiles 22, 811–823. https://doi.org/10.1007/s00792-018-1040-9 (2018).CAS 
Article 
PubMed 

Google Scholar 
Tazi, L., Breakwell, D. P., Harker, A. R. & Crandall, K. A. Life in extreme environments: Microbial diversity in Great Salt Lake, Utah. Extremophiles 18, 525–535. https://doi.org/10.1007/s00792-014-0637-x (2014).Article 
PubMed 

Google Scholar 
Kashi, F. J., Owlia, P., Amoozegar, M. A. & Kazemi, B. Halophilic prokaryotes in Urmia Salt Lake, a hypersaline environment in Iran. Curr. Microbiol. 78(8), 3230–3238 (2021).Article 

Google Scholar 
Sorokin, D. Y., Roman, P. & Kolganova, T. V. Halo(natrono)archaea from hypersaline lakes can utilize sulfoxides other than DMSO as electron acceptors for anaerobic respiration. Extremophiles 25, 173–180 (2021).CAS 
Article 

Google Scholar 
Hwang, K., Choe, H. & Kim, K. M. Complete genome of Nocardioides aquaticus KCTC 9944T isolated from meromictic and hypersaline Ekho Lake, Antarctica. Mar. Genom. 1, 100889 (2021).Article 

Google Scholar 
Didari, M. et al. Diversity of halophilic and halotolerant bacteria in the largest seasonal hypersaline lake (Aran-Bidgol-Iran). J. Environ. Health Sci. Eng. 18, 961–971. https://doi.org/10.1007/s40201-020-00519-3 (2020).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Oren, A. Diversity of halophilic microorganisms: Environments, phylogeny, physiology, and applications. J. Ind. Microbiol. Biotechnol. 28, 56–63 (2002).CAS 
Article 

Google Scholar 
Mutlu, M. B. et al. Prokaryotic diversity in Tuz Lake, a hypersaline environment in Inland Turkey. FEMS Microbiol. Ecol. 65, 474–483. https://doi.org/10.1111/j.1574-6941.2008.00510.x (2008).CAS 
Article 
PubMed 

Google Scholar 
Antón, J. et al. Distribution, abundance and diversity of the extremely halophilic bacterium Salinibacter ruber. Saline Syst. 4, 15. https://doi.org/10.1186/1746-1448-4-15 (2008).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Oren, A. Microbial life at high salt concentrations: phylogenetic and metabolic diversity. Saline Syst. 4, 2. https://doi.org/10.1186/1746-1448-4-2 (2008).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Abdeljabbar, H., Badiaa, E., Jean-Luc, C., Marie-Laure, F. & Najla, S. Prokaryotic biodiversity of halophilic microorganisms isolated from Sehline Sebkha Salt Lake (Tunisia). Afr. J. Microbiol. Res. 8, 355–367. https://doi.org/10.5897/ajmr12.1087 (2014).Article 

Google Scholar 
Najjari, A., Elshahed, M. S., Cherif, A., Youssef, N. H. & Löffler, F. E. Patterns and determinants of halophilic archaea (Class Halobacteria) diversity in Tunisian endorheic salt lakes and Sebkhet systems. Appl. Environ. Microbiol. 81, 4432–4441. https://doi.org/10.1128/aem.01097-15 (2015).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Aharon, O. The ecology of the extremely halophilic archaea. FEMS Microbiol. Rev. 1, 415–440 (1994).
Google Scholar 
Oren, A. Halophilic Archaea. FEMS Microbiol. Rev. https://doi.org/10.1016/b978-0-12-809633-8.20800-5 (2019).Article 

Google Scholar 
Feng, Y. et al. The evolutionary origins of extreme halophilic archaeal lineages. Genome Biol. Evol. 13, 8. https://doi.org/10.1093/gbe/evab166 (2021).CAS 
Article 

Google Scholar 
Ventosa, A., Nieto, J. J. & Oren, A. Biology of moderately halophilic aerobic bacteria. Microbiol. Mol. Biol. Rev. 62, 504–544 (1998).CAS 
Article 

Google Scholar 
Kushner, D. J. Halophilic bacteria. Adv. Appl. Microbiol. 10, 73–99 (1968).CAS 
Article 

Google Scholar 
Ghozlan, H., Deif, H., Kandil, R. A. & Sabry, S. Biodiversity of moderately halophilic bacteria in hypersaline habitats in Egypt. J. Gen. Appl. Microbiol. 52, 63–72 (2006).CAS 
Article 

Google Scholar 
Ali, I., Prasongsuk, S., Akbar, A., Aslam, M. & Rakshit, S. K. Hypersaline habitats and halophilic microorganisms. Maejo Int. J. Sci. Technol. 10, 330–345 (2016).CAS 

Google Scholar 
Margesin, R. & Schinner, F. Biodegradation and bioremediation of hydrocarbons in extreme environments. Appl. Microbiol. Biotechnol. 56, 650–663. https://doi.org/10.1007/s002530100701 (2001).CAS 
Article 
PubMed 

Google Scholar 
Poosarla, V. G. & Ts, C. Xylanase production by halophilic bacterium Gracilibacillus sp. TSCPVG under solid state fermentation. Res. J. Biotechnol. 16, 92–100 (2021).Article 

Google Scholar 
Foti, M. et al. Diversity, activity, and abundance of sulfate-reducing bacteria in saline and hypersaline soda lakes. Appl. Environ. Microbiol. 73, 2093–2100. https://doi.org/10.1128/aem.02622-06 (2007).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Boujelben, I. et al. Spatial and seasonal prokaryotic community dynamics in ponds of increasing salinity of Sfax solar saltern in Tunisia. Antonie Van Leeuwenhoek 101, 845–857. https://doi.org/10.1007/s10482-012-9701-7 (2012).CAS 
Article 
PubMed 

Google Scholar 
García-Maldonado, J. Q., Bebout, B. M., Everroad, R. C. & López-Cortés, A. Evidence of novel phylogenetic lineages of methanogenic archaea from hypersaline microbial mats. Microb. Ecol. 69, 106–117. https://doi.org/10.1007/s00248-014-0473-7 (2014).CAS 
Article 
PubMed 

Google Scholar 
Abed, R. M. M., de Beer, D. & Stief, P. Functional-structural analysis of nitrogen-cycle bacteria in a hypersaline mat from the omani desert. Geomicrobiol. J. 32, 119–129. https://doi.org/10.1080/01490451.2014.932033 (2014).CAS 
Article 

Google Scholar 
Coban, O., Rasigraf, O., Jong, A., Spott, O. & Bebout, B. M. Quantifying potential N turnover rates in hypersaline microbial mats by 15 N tracer techniques. Appl. Environ. Microbiol. 87, 8 (2021).Article 

Google Scholar 
Rodriguez-Valera, F. Introduction to Saline Environments (Springer, 1993).
Google Scholar 
Wei, H. C., Qi-Shun, F., Fu-Yuan, A., Fa-Shou, S. & Qin, Z. J. Chemical elements in core sediments of the qarhan salt lake and palaeoclimate evolution during 94–9 ka. Acta Geosci. Sin. (2016).Yu, S., Liu, X., Tan, H. & Cao, G. Sustainable Utilization of Qarhan Salt Lake Resources 27–265 (Science Press, 2009).
Google Scholar 
Zhu, D. et al. An evaluation of the core bacterial communities associated with hypersaline environments in the Qaidam Basin, China. Arch. Microbiol. 202, 2093–2103. https://doi.org/10.1007/s00203-020-01927-7 (2020).CAS 
Article 
PubMed 

Google Scholar 
Liu, W., Jiang, H., Yang, J. & Wu, G. Gammaproteobacterial diversity and carbon utilization in response to salinity in the lakes on the qinghai-tibetan plateau. Geomicrobiol. J. 35, 392–403. https://doi.org/10.1080/01490451.2017.1378951 (2018).CAS 
Article 

Google Scholar 
Zhong, Z.-P. et al. Prokaryotic community structure driven by salinity and ionic concentrations in plateau lakes of the tibetan plateau. Appl. Environ. Microbiol. 82, 1846–1858. https://doi.org/10.1128/aem.03332-15 (2016).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
He, C. et al. Synergistic effect of magnetite and zero-valent iron on anaerobic degradation and methanogenesis of phenol. Biores. Technol. 291, 121874. https://doi.org/10.1016/j.biortech.2019.121874 (2019).CAS 
Article 

Google Scholar 
Caporaso, J. G. et al. QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. https://doi.org/10.1038/nmeth.f.303 (2010).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Martin, M. Cutadapt removes adapter sequences from high-throughput sequencing reads. Embnet J. 17, 10–12 (2011).Article 

Google Scholar 
Zhang, J., Kassian, K., Tomáš, F. & Alexandros, S. PEAR: a fast and accurate Illumina Paired-End reAd mergeR. Bioinformatics 30, 614 (2014).CAS 
Article 

Google Scholar 
Schmieder, R. & Edwards, R. Quality control and preprocessing of metagenomic datasets. Bioinformatics 27, 863–864 (2011).CAS 
Article 

Google Scholar 
Edgar, R. C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 10, 996–1000 (2013).CAS 
Article 

Google Scholar 
Schloss, P. D. et al. Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 75, 7537 (2009).ADS 
CAS 
Article 

Google Scholar 
Chen, H. & Boutros, P. C. VennDiagram: A package for the generation of highly-customizable Venn and Euler diagrams in R. BMC Bioinform. 12, 35. https://doi.org/10.1186/1471-2105-12-35 (2011).Article 

Google Scholar 
McArdle, B. H. et al. Fitting multivariate models to community data: A comment on distance-based redundancy analysis. Ecology 82, 290–290 (2001).Article 

Google Scholar 
Langille, M. G. I. et al. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat. Biotechnol. 31, 814–821. https://doi.org/10.1038/nbt.2676 (2013).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Louca, S. & Doebeli, M. Efficient comparative phylogenetics on large trees. Bioinformatics 34, 1–3 (2017).
Google Scholar 
Louca, S., Parfrey, L. W. & Doebeli, M. Decoupling function and taxonomy in the global ocean microbiome. Science 353, 1272 (2016).ADS 
CAS 
Article 

Google Scholar 
Junker, B. H. & Schreiber, F. Analysis of Biological Networks 283–304 (Analysis of biological networks, 2008).Book 

Google Scholar 
Faust, K. & Raes, J. Microbial interactions: From networks to models. Nat. Rev. Microbiol. 10, 538–550. https://doi.org/10.1038/nrmicro2832 (2012).CAS 
Article 
PubMed 

Google Scholar 
Behzad, H., Ibarra, M. A., Mineta, K. & Gojobori, T. Metagenomic studies of the Red Sea. Gene 576, 717–723. https://doi.org/10.1016/j.gene.2015.10.034 (2016).CAS 
Article 
PubMed 

Google Scholar 
Naghoni, A. et al. Microbial diversity in the hypersaline Lake Meyghan, Iran. Sci. Rep. https://doi.org/10.1038/s41598-017-11585-3 (2017).Article 
PubMed 
PubMed Central 

Google Scholar 
Kambura, A. K. et al. Bacteria and Archaea diversity within the hot springs of Lake Magadi and Little Magadi in Kenya. BMC Microbiol. https://doi.org/10.1186/s12866-016-0748-x (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Paul, D. et al. Exploration of microbial diversity and community structure of Lonar lake: The only hypersaline meteorite crater lake within basalt rock. Front. Microbiol. https://doi.org/10.3389/fmicb.2015.01553 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Ventosa, A., de la Haba, R. R., Sánchez-Porro, C. & Papke, R. T. Microbial diversity of hypersaline environments: A metagenomic approach. Curr. Opin. Microbiol. 25, 80–87. https://doi.org/10.1016/j.mib.2015.05.002 (2015).CAS 
Article 
PubMed 

Google Scholar 
Liu, F. H. et al. Bacterial and archaeal assemblages in sediments of a large shallow freshwater lake, Lake Taihu, as revealed by denaturing gradient gel electrophoresis. J. Appl. Microbiol. 106, 1022–1032. https://doi.org/10.1111/j.1365-2672.2008.04069.x (2009).CAS 
Article 
PubMed 

Google Scholar 
Song, H., Li, Z., Du, B., Wang, G. & Ding, Y. Bacterial communities in sediments of the shallow Lake Dongping in China. J. Appl. Microbiol. 112, 79–89. https://doi.org/10.1111/j.1365-2672.2011.05187.x (2012).CAS 
Article 
PubMed 

Google Scholar 
Wu, Q. L., Zwart, G., Schauer, M., Agterveld, K. V. & Hahn, M. W. Bacterioplankton community composition along a salinity gradient of sixteen high-mountain lakes located on the Tibetan Plateau, China. Appl. Environ. Microbiol. 72, 5478–5485 (2006).ADS 
CAS 
Article 

Google Scholar 
Xing, P., Hahn, M. W. & Wu, Q. L. Low taxon richness of bacterioplankton in high-altitude lakes of the Eastern Tibetan Plateau, with a predominance of bacteroidetes and Synechococcus spp. Appl. Environ. Microbiol. 75, 7017–7025. https://doi.org/10.1128/aem.01544-09 (2009).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Liu, Y. et al. Bacterial diversity of freshwater Alpine Lake Puma Yumco on the Tibetan Plateau. Geomicrobiol. J. 26, 131–145. https://doi.org/10.1080/01490450802660201 (2009).CAS 
Article 

Google Scholar 
MounÃc, S., Caumette, P., Matheron, R. & Willison, J. C. Molecular sequence analysis of prokaryotic diversity in the anoxic sediments underlying cyanobacterial mats of two hypersaline ponds in Mediterranean salterns. FEMS Microbiol. Ecol. 44, 117–130. https://doi.org/10.1016/s0168-6496(03)00017-5 (2003).Article 

Google Scholar 
Valenzuela-Encinas, C. et al. Changes in the bacterial populations of the highly alkaline saline soil of the former lake Texcoco (Mexico) following flooding. Extremophiles 13, 609–621. https://doi.org/10.1007/s00792-009-0244-4 (2009).Article 
PubMed 

Google Scholar 
Kim, T. J., Lee, E. Y., Kim, Y. J., Cho, K.-S. & Ryu, H. W. Degradation of polyaromatic hydrocarbons by Burkholderia cepacia 2A–12. World J. Microbiol. Biotechnol. 19, 411–417. https://doi.org/10.1023/A:1023998719787 (2003).CAS 
Article 

Google Scholar 
Gales, G. et al. Preservation of ancestral Cretaceous microflora recovered from a hypersaline oil reservoir. Sci. Rep. https://doi.org/10.1038/srep22960 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Kleinsteuber, S., Riis, V., Fetzer, I., Harms, H. & Müller, S. Population dynamics within a microbial consortium during growth on diesel fuel in saline environments. Appl. Environ. Microbiol. 72, 3531–3542. https://doi.org/10.1128/aem.72.5.3531-3542.2006 (2006).ADS 
CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
Valenzuela-Encinas, C. et al. The archaeal diversity and population in a drained alkaline saline soil of the former lake Texcoco (Mexico). Geomicrobiol. J. 29, 18–22. https://doi.org/10.1080/01490451.2010.520075 (2012).Article 

Google Scholar 
He, S., Tan, J., Hu, W. & Mo, C. Diversity of archaea and its correlation with environmental factors in the Ebinur Lake Wetland. Curr. Microbiol. 76, 1417–1424. https://doi.org/10.1007/s00284-019-01768-8 (2019).CAS 
Article 
PubMed 

Google Scholar 
Sandaa, R. A., Enger, O. & Torsvik, V. Abundance and diversity of Archaea in heavy-metal-contaminated soils. Appl. Environ. Microbiol. 65, 3293–3297 (1999).ADS 
CAS 
Article 

Google Scholar 
Dave, B. P. & Soni, A. Diversity of halophilic archaea at salt pans around Bhavnagar Coast, Gujarat. Proc. Natl. Acad. Sci. India B 83, 225–232. https://doi.org/10.1007/s40011-012-0124-z (2012).Article 

Google Scholar 
Zafrilla, B., Martínez-Espinosa, R., Alonso, M. A. & Bonete, M. J. Biodiversity of Archaea and floral of two inland saltern ecosystems in the Alto Vinalopó Valley, Spain. Saline Syst. 6, 10 (2010).Article 

Google Scholar 
Costa, M., Santos, H. & Galinski, E. A. An overview of the role and diversity of compatible solutes in Bacteria and Archaea. Adv. Biochem. Eng. Biotechnol. 61, 117 (1998).PubMed 

Google Scholar 
Williams, R. J., Howe, A. & Hofmockel, K. S. Demonstrating microbial co-occurrence pattern analyses within and between ecosystems. Front. Microbiol. https://doi.org/10.3389/fmicb.2014.00358 (2014).Article 
PubMed 
PubMed Central 

Google Scholar 
Schmidt, T. S. B., MatiasRodrigues, J. F. & von Mering, C. A family of interaction-adjusted indices of community similarity. ISME J. 11, 791–807. https://doi.org/10.1038/ismej.2016.139 (2016).Article 
PubMed 
PubMed Central 

Google Scholar 
Oyewusi, H. A. et al. Functional profiling of bacterial communities in Lake Tuz using 16S rRNA gene sequences. Biotechnol. Biotechnol. Equip. 35, 1–10. https://doi.org/10.1080/13102818.2020.1840437 (2020).CAS 
Article 

Google Scholar  More

Taylor, M. W., Radax, R., Steger, D. & Wagner, M. Sponge-associated microorganisms: Evolution, ecology, and biotechnological potential. Microbiol. Mol. Biol. Rev. 71, 295–347 (2007).CAS 
PubMed 
PubMed Central 

Google Scholar 
Thomas, T. et al. Diversity, structure and convergent evolution of the global sponge microbiome. Nat. Commun. 7, 11870 (2016).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Webster, N. S. et al. Deep sequencing reveals exceptional diversity and modes of transmission for bacterial sponge symbionts. Environ. Microbiol. 12, 2070–2082 (2010).CAS 
PubMed 
PubMed Central 

Google Scholar 
Sipkema, D. et al. Similar sponge-associated bacteria can be acquired via both vertical and horizontal transmission: Microbial transmission in Petrosia ficiformis. Environ. Microbiol. 17, 3807–3821 (2015).CAS 
PubMed 

Google Scholar 
Cleary, D. F. R. et al. The sponge microbiome within the greater coral reef microbial metacommunity. Nat. Commun. 10, 1644 (2019).Björk, J. R., Díez-Vives, C., Astudillo-García, C., Archie, E. A. & Montoya, J. M. Vertical transmission of sponge microbiota is inconsistent and unfaithful. Nat. Ecol. Evol. 3, 1172–1183 (2019).PubMed 
PubMed Central 

Google Scholar 
Webster, N. S. & Taylor, M. W. Marine sponges and their microbial symbionts: Love and other relationships. Environ. Microbiol. 14, 335–346 (2012).CAS 
PubMed 

Google Scholar 
Kennedy, J. et al. Evidence of a putative deep sea specific microbiome in marine sponges. PLoS ONE 9, e91092 (2014).ADS 
PubMed 
PubMed Central 

Google Scholar 
Steinert, G. et al. Compositional and quantitative insights into bacterial and archaeal communities of south pacific deep-sea sponges (Demospongiae and Hexactinellida). Front. Microbiol. 11, 716 (2020).Busch, K. et al. On giant shoulders: How a seamount affects the microbial community composition of seawater and sponges. Biogeosciences 17, 3471–3486 (2020).ADS 
CAS 

Google Scholar 
Olson, J. B. & Gao, X. Characterizing the bacterial associates of three Caribbean sponges along a gradient from shallow to mesophotic depths. FEMS Microbiol. Ecol. 85, 74–84 (2013).PubMed 

Google Scholar 
Steinert, G. et al. In four shallow and mesophotic tropical reef sponges from Guam the microbial community largely depends on host identity. PeerJ 4, e1936 (2016).PubMed 
PubMed Central 

Google Scholar 
Morrow, K. M., Fiore, C. L. & Lesser, M. P. Environmental drivers of microbial community shifts in the giant barrel sponge, Xestospongia muta, over a shallow to mesophotic depth gradient. Environ. Microbiol. 18, 2025–2038 (2016).CAS 
PubMed 

Google Scholar 
Ebada, S. S. & Proksch, P. The chemistry of marine sponges. In Handbook of Marine Natural Products (eds Fattorusso, E. et al.) 191–293 (Springer, 2012). https://doi.org/10.1007/978-90-481-3834-0_4.Chapter 

Google Scholar 
Kornprobst, J.-M. Porifera (Sponges). Encyclopedia of Marine Natural Products (Wiley, 2014).
Google Scholar 
Leal, M. C., Puga, J., Serôdio, J., Gomes, N. C. M. & Calado, R. Trends in the discovery of new marine natural products from invertebrates over the last two decades—Where and what are we bioprospecting?. PLoS ONE 7, e30580 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Blunt, J. W., Copp, B. R., Keyzers, R. A., Munro, M. H. G. & Prinsep, M. R. Marine natural products. Nat. Prod. Rep. 34, 235–294 (2017).CAS 
PubMed 

Google Scholar 
Unson, M. D., Holland, N. D. & Faulkner, D. J. A brominated secondary metabolite synthesized by the cyanobacterial symbiont of a marine sponge and accumulation of the crystalline metabolite in the sponge tissue. Mar. Biol. 119, 1–11 (1994).CAS 

Google Scholar 
Bewley, C. A., Holland, N. D. & Faulkner, D. J. Two classes of metabolites from Theonella swinhoei are localized in distinct populations of bacterial symbionts. Experientia 52, 716–722 (1996).CAS 
PubMed 

Google Scholar 
Wilson, M. C. et al. An environmental bacterial taxon with a large and distinct metabolic repertoire. Nature 506, 58–62 (2014).ADS 
CAS 
PubMed 

Google Scholar 
Tianero, M. D., Balaich, J. N. & Donia, M. S. Localized production of defence chemicals by intracellular symbionts of Haliclona sponges. Nat. Microbiol. 4, 1149–1159 (2019).CAS 
PubMed 

Google Scholar 
Ivanišević, J., Thomas, O. P., Lejeusne, C., Chevaldonné, P. & Pérez, T. Metabolic fingerprinting as an indicator of biodiversity: Towards understanding inter-specific relationships among Homoscleromorpha sponges. Metabolomics 7, 289–304 (2011).
Google Scholar 
Pérez, T. et al. Oscarella balibaloi, a new sponge species (Homoscleromorpha: Plakinidae) from the Western Mediterranean Sea: Cytological description, reproductive cycle and ecology: O. balibaloi: Description, reproductive cycle and ecology. Mar. Ecol. (Berl.) 32, 174–187 (2011).ADS 

Google Scholar 
Reveillaud, J. et al. Relevance of an integrative approach for taxonomic revision in sponge taxa: Case study of the shallow-water Atlanto-Mediterranean Hexadella species (Porifera: Ianthellidae: Verongida). Invertebr. Syst. 26, 230–248 (2012).
Google Scholar 
Olsen, E. K. et al. Marine AChE inhibitors isolated from Geodia barretti: Natural compounds and their synthetic analogs. Org. Biomol. Chem. 14, 1629–1640 (2016).CAS 
PubMed 

Google Scholar 
Reverter, M., Perez, T., Ereskovsky, A. V. & Banaigs, B. Secondary metabolome variability and inducible chemical defenses in the Mediterranean Sponge Aplysina cavernicola. J. Chem. Ecol. 42, 60–70 (2016).CAS 
PubMed 

Google Scholar 
Reverter, M., Tribalat, M.-A., Pérez, T. & Thomas, O. P. Metabolome variability for two Mediterranean sponge species of the genus Haliclona: Specificity, time, and space. Metabolomics 14, 114 (2018).Villegas-Plazas, M. et al. Variations in microbial diversity and metabolite profiles of the tropical marine sponge Xestospongia muta with season and depth. Microb. Ecol. 78, 243–256 (2019).CAS 
PubMed 

Google Scholar 
Mohanty, I. et al. Multi-omic profiling of Melophlus sponges reveals diverse metabolomic and microbiome architectures that are non-overlapping with ecological neighbors. Mar. Drugs 18, 124 (2020).CAS 
PubMed Central 

Google Scholar 
Bowerbank, J. S. On the anatomy and physiology of the Spongiadae. Part I. On the spicula. Philos. Trans. R. Soc. Lond. 148, 279–332 (1858).ADS 

Google Scholar 
Vosmaer, G. C. J. The sponges of the ‘Willem Barents’ expedition 1880 and 1881. Bijdragen tot de Dierkunde 12, 1–47 (1885).
Google Scholar 
Radax, R. et al. Metatranscriptomics of the marine sponge Geodia barretti: Tackling phylogeny and function of its microbial community. Environ. Microbiol. 14, 1308–1324 (2012).CAS 
PubMed 

Google Scholar 
Topsent, E. Spongiaires provenant des campagnes scientifiques de la ‘Princesse Alice’ dans les Mers du Nord (1898–1899—1906–1907). Résultats des campagnes scientifiques accomplies par le Prince Albert I. Monaco 45, 1–67 (1913).
Google Scholar 
Yashayaev, I. & Loder, J. W. Further intensification of deep convection in the Labrador Sea in 2016. Geophys. Res. Lett. 44, 1429–1438 (2017).ADS 

Google Scholar 
Gutleben, J. et al. Diversity of tryptophan halogenases in sponges of the genus Aplysina. FEMS Microbiol. Ecol. 95, fiz108 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Indraningrat, A. et al. Cultivation of sponge-associated bacteria from Agelas sventres and Xestospongia muta collected from different depths. Mar. Drugs 17, 578 (2019).CAS 
PubMed Central 

Google Scholar 
Ramiro-Garcia, J. et al. NG-Tax, a highly accurate and validated pipeline for analysis of 16S rRNA amplicons from complex biomes. F1000 Res. 5, 1791 (2018).
Google Scholar 
Yilmaz, P. et al. The SILVA and “All-species Living Tree Project (LTP)” taxonomic frameworks. Nucl. Acids Res. 42, D643–D648 (2014).CAS 
PubMed 

Google Scholar 
Erngren, I., Smit, E., Pettersson, C., Cárdenas, P. & Hedeland, M. The effects of sampling and storage conditions on the metabolite profile of the marine sponge Geodia barretti. Front. Chem. 9:662659 (2021)Smith, C. A., Want, E. J., O’Maille, G., Abagyan, R. & Siuzdak, G. XCMS: Processing mass spectrometry data for metabolite profiling using nonlinear peak alignment, matching, and identification. Anal. Chem. 78, 779–787 (2006).CAS 
PubMed 

Google Scholar 
Kuhl, C., Tautenhahn, R., Böttcher, C., Larson, T. R. & Neumann, S. CAMERA: An integrated strategy for compound spectra extraction and annotation of liquid chromatography/mass spectrometry data sets. Anal. Chem. 84, 283–289 (2012).CAS 
PubMed 

Google Scholar 
Oksanen, J. et al. vegan: Community Ecology Package (2017).Dat, T. T. H., Steinert, G., Thi Kim Cuc, N., Smidt, H. & Sipkema, D. Archaeal and bacterial diversity and community composition from 18 phylogenetically divergent sponge species in Vietnam. PeerJ 6, e4970 (2018).PubMed 
PubMed Central 

Google Scholar 
Miller, M. A., Pfeiffer, W. & Schwartz, T. Creating the CIPRES science gateway for inference of large phylogenetic trees. In 2010 Gateway Computing Environments Workshop (GCE) 1–8 (IEEE, 2010). https://doi.org/10.1109/GCE.2010.5676129.Letunic, I. & Bork, P. Interactive Tree Of Life (iTOL) v4: Recent updates and new developments. Nucl. Acids Res. 47, W256–W259 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
Thévenot, E. A., Roux, A., Xu, Y., Ezan, E. & Junot, C. Analysis of the human adult urinary metabolome variations with age, body mass index, and gender by implementing a comprehensive workflow for univariate and OPLS statistical analyses. J. Proteome Res. 14, 3322–3335 (2015).PubMed 

Google Scholar 
Weiss, S. et al. Correlation detection strategies in microbial data sets vary widely in sensitivity and precision. ISME J. 10, 1669–1681 (2016).CAS 
PubMed 
PubMed Central 

Google Scholar 
Deng, Y. et al. Molecular ecological network analyses. BMC Bioinform. 13, 113 (2012).
Google Scholar 
Friedman, J. & Alm, E. J. Inferring correlation networks from genomic survey data. PLoS Comput. Biol. 8, e1002687 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Durno, W. E., Hanson, N. W., Konwar, K. M. & Hallam, S. J. Expanding the boundaries of local similarity analysis. BMC Genom. 14, S3 (2013).
Google Scholar 
Reshef, D. N. et al. Detecting novel associations in large data sets. Science 334, 1518–1524 (2011).ADS 
CAS 
PubMed 
PubMed Central 
MATH 

Google Scholar 
Hall, M. M., Torres, D. J. & Yashayaev, I. Absolute velocity along the AR7W section in the Labrador Sea. Deep Sea Res. Part 1 Oceanogr. Res. Pap. 72, 72–87 (2013).
Google Scholar 
Reveillaud, J. et al. Host-specificity among abundant and rare taxa in the sponge microbiome. ISME J. 8, 1198–1209 (2014).CAS 
PubMed 
PubMed Central 

Google Scholar 
Moitinho-Silva, L. et al. Predicting the HMA-LMA status in marine sponges by machine learning. Front. Microbiol. 8, 752 (2017).Lidgren, G., Bohlin, L. & Bergman, J. Studies of Swedish marine organisms VII. A novel biologically active indole alkaloid from the sponge Geodia barretti. Tetrahedron Lett. 27, 3283–3284 (1986).CAS 

Google Scholar 
Sjögren, M. et al. Antifouling activity of brominated cyclopeptides from the marine sponge Geodia barretti. J. Nat. Prod. 67, 368–372 (2004).PubMed 

Google Scholar 
Sölter, S. Identifizierung und Synthese von Naturstoffen aus Borealen Schwämmen (Universität Hamburg, 2004).
Google Scholar 
Di, X. et al. 6-Bromoindole derivatives from the Icelandic marine sponge Geodia barretti: Isolation and anti-inflammatory activity. Mar. Drugs 16, 437 (2018).CAS 
PubMed Central 

Google Scholar 
Carstens, B. B. et al. Isolation, characterization, and synthesis of the barrettides: Disulfide-containing peptides from the marine sponge Geodia barretti. J. Nat. Prod. 78, 1886–1893 (2015).CAS 
PubMed 

Google Scholar 
Hedner, E. et al. Brominated cyclodipeptides from the marine sponge Geodia barretti as selective 5-HT ligands. J. Nat. Prod. 69, 1421–1424 (2006).CAS 
PubMed 

Google Scholar 
Hedner, E. et al. Antifouling activity of a dibrominated cyclopeptide from the marine sponge Geodia barretti. J. Nat. Prod. 71, 330–333 (2008).CAS 
PubMed 

Google Scholar 
Erwin, P. M., Pita, L., López-Legentil, S. & Turon, X. Stability of sponge-associated bacteria over large seasonal shifts in temperature and irradiance. Appl. Environ. Microbiol. 78, 7358–7368 (2012).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Cárdenas, C. A., Bell, J. J., Davy, S. K., Hoggard, M. & Taylor, M. W. Influence of environmental variation on symbiotic bacterial communities of two temperate sponges. FEMS Microbiol. Ecol. 88, 516–527 (2014).PubMed 

Google Scholar 
Glasl, B., Smith, C. E., Bourne, D. G. & Webster, N. S. Exploring the diversity-stability paradigm using sponge microbial communities. Sci. Rep. 8, 8425 (2018).Schöttner, S. et al. Relationships between host phylogeny, host type and bacterial community diversity in cold-water coral reef sponges. PLoS ONE 8, e55505 (2013).ADS 
PubMed 
PubMed Central 

Google Scholar 
Lurgi, M., Thomas, T., Wemheuer, B., Webster, N. S. & Montoya, J. M. Modularity and predicted functions of the global sponge-microbiome network. Nat. Commun. 10, 992 (2019).ADS 
PubMed 
PubMed Central 

Google Scholar 
Luter, H. M. et al. Microbiome analysis of a disease affecting the deep-sea sponge Geodia barretti. FEMS Microbiol. Ecol. 93, fix074 (2017).Thistle, D. Ecosystems of the Deep Oceans (Elsevier, 2003).
Google Scholar 
Pita, L., Erwin, P. M., Turon, X. & López-Legentil, S. Till death do us part: Stable sponge-bacteria associations under thermal and food shortage stresses. PLoS ONE 8, e80307 (2013).ADS 
PubMed 
PubMed Central 

Google Scholar 
Webster, N. S., Cobb, R. E. & Negri, A. P. Temperature thresholds for bacterial symbiosis with a sponge. ISME J. 2, 830–842 (2008).CAS 
PubMed 

Google Scholar 
Gerringer, M. E., Drazen, J. C. & Yancey, P. H. Metabolic enzyme activities of abyssal and hadal fishes: Pressure effects and a re-evaluation of depth-related changes. Deep Sea Res. Part 1 Oceanogr. Res. Pap. 125, 135–146 (2017).CAS 

Google Scholar 
Yashayaev, I. Hydrographic changes in the Labrador Sea, 1960–2005. Prog. Oceanogr. 73, 242–276 (2007).ADS 

Google Scholar 
Rhein, M., Steinfeldt, R., Kieke, D., Stendardo, I. & Yashayaev, I. Ventilation variability of Labrador Sea Water and its impact on oxygen and anthropogenic carbon: A review. Philos. Trans. A Math. Phys. Eng. Sci. 375, 20160321 (2017).ADS 
PubMed 
PubMed Central 

Google Scholar 
Galand, P. E., Potvin, M., Casamayor, E. O. & Lovejoy, C. Hydrography shapes bacterial biogeography of the deep Arctic Ocean. ISME J. 4, 564–576 (2010).PubMed 

Google Scholar 
Frank, A. H., Garcia, J. A. L., Herndl, G. J. & Reinthaler, T. Connectivity between surface and deep waters determines prokaryotic diversity in the North Atlantic Deep Water: North Atlantic dark ocean prokaryotic biogeography. Environ. Microbiol. 18, 2052–2063 (2016).PubMed 
PubMed Central 

Google Scholar 
Agogué, H., Lamy, D., Neal, P. R., Sogin, M. L. & Herndl, G. J. Water mass-specificity of bacterial communities in the North Atlantic revealed by massively parallel sequencing. Mol. Ecol. 20, 258–274 (2011).PubMed 

Google Scholar 
Djurhuus, A., Boersch-Supan, P. H., Mikalsen, S.-O. & Rogers, A. D. Microbe biogeography tracks water masses in a dynamic oceanic frontal system. R. Soc. Open Sci. 4, 170033 (2017).ADS 
PubMed 
PubMed Central 

Google Scholar 
Müller, O. et al. Spatiotemporal dynamics of ammonia-oxidizing Thaumarchaeota in distinct Arctic water masses. Front. Microbiol. 9, 1–13 (2018).ADS 

Google Scholar 
Kraemer, S., Ramachandran, A., Colatriano, D., Lovejoy, C. & Walsh, D. A. Diversity and biogeography of SAR11 bacteria from the Arctic Ocean. ISME J. https://doi.org/10.1038/s41396-019-0499-4 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
Monier, A. et al. Upper Arctic Ocean water masses harbor distinct communities of heterotrophic flagellates. Biogeosciences 10, 4273–4286 (2013).ADS 

Google Scholar 
Monier, A. et al. Oceanographic structure drives the assembly processes of microbial eukaryotic communities. ISME J. 9, 990–1002 (2015).CAS 
PubMed 

Google Scholar 
Corrège, T. The relationship between water masses and benthic ostracod assemblages in the western Coral Sea, Southwest Pacific. Palaeogeogr. Palaeoclimatol. Palaeoecol. 105, 245–266 (1993).
Google Scholar 
Muhling, B. A., Beckley, L. E., Koslow, J. A. & Pearce, A. F. Larval fish assemblages and water mass structure off the oligotrophic south-western Australian coast: SW Australian larval fish assemblages. Fish. Oceanogr. 17, 16–31 (2007).
Google Scholar 
Eerkes-Medrano, D. et al. A community assessment of the demersal fish and benthic invertebrates of the Rosemary Bank Seamount Marine Protected Area (NE Atlantic). Deep Sea Res. Part 1 Oceanogr. Res. Pap. https://doi.org/10.1016/j.dsr.2019.103180 (2019).Article 

Google Scholar 
Puerta, P. et al. Influence of water masses on the biodiversity and biogeography of deep-sea benthic ecosystems in the North Atlantic. Front. Mar. Sci. 7, 239 (2020).Roberts, E. et al. Water masses constrain the distribution of deep-sea sponges in the North Atlantic Ocean and Nordic Seas. Mar. Ecol. Prog. Ser. 659, 75–96 (2021).ADS 

Google Scholar 
Kenchington, E. et al. Connectivity modelling of areas closed to protect vulnerable marine ecosystems in the northwest Atlantic. Deep Sea Res. Part 1 Oceanogr. Res. Pap. 143, 85–103 (2019).
Google Scholar 
Louca, S. et al. Function and functional redundancy in microbial systems. Nat. Ecol. Evol. 2, 936–943 (2018).PubMed 

Google Scholar 
McCauley, M., Chiarello, M., Atkinson, C. L. & Jackson, C. R. Gut microbiomes of freshwater mussels (Unionidae) are taxonomically and phylogenetically variable across years but remain functionally stable. Microorganisms 9, 411 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Page, M., West, L., Northcote, P., Battershill, C. & Kelly, M. Spatial and temporal variability of cytotoxic metabolites in populations of the New Zealand Sponge Mycale hentscheli. J. Chem. Ecol. 31, 1161–1174 (2005).CAS 
PubMed 

Google Scholar 
Ternon, E., Perino, E., Manconi, R., Pronzato, R. & Thomas, O. P. How environmental factors affect the production of guanidine alkaloids by the Mediterranean sponge Crambe crambe. Mar. Drugs 15, 181 (2017).PubMed Central 

Google Scholar 
Sacristán-Soriano, O., Banaigs, B. & Becerro, M. A. Temporal trends in the secondary metabolite production of the sponge Aplysina aerophoba. Mar. Drugs 10, 677–693 (2012).PubMed 
PubMed Central 

Google Scholar 
Ivanisevic, J. et al. Biochemical trade-offs: Evidence for ecologically linked secondary metabolism of the sponge Oscarella balibaloi. PLoS ONE 6, e28059 (2011).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Burg, M. B. & Ferraris, J. D. Intracellular organic osmolytes: Function and regulation. J. Biol. Chem. 283, 7309–7313 (2008).CAS 
PubMed 
PubMed Central 

Google Scholar 
Nau-Wagner, G., Boch, J., Le Good, J. A. & Bremer, E. High-affinity transport of choline-O-sulfate and its use as a compatible solute in Bacillus subtilis. Appl. Environ. Microbiol. 65, 560–568 (1999).ADS 
CAS 
PubMed 
PubMed Central 

Google Scholar 
Popowich, A., Zhang, Q. & Le, X. C. Arsenobetaine: The ongoing mystery. Natl. Sci. Rev. 3, 451–458 (2016).CAS 

Google Scholar 
Connor, K. M. & Gracey, A. Y. High-resolution analysis of metabolic cycles in the intertidal mussel Mytilus californianus. Am. J. Physiol. Regul. Integr. Comp. Physiol. 302, R103–R111 (2012).CAS 
PubMed 

Google Scholar 
Cárdenas, P. Who produces Ianthelline? The Arctic sponge Stryphnus fortis or its sponge Epibiont Hexadella dedritifera: A probable case of sponge–sponge contamination. J. Chem. Ecol. 42, 339–347 (2016).PubMed 

Google Scholar 
Steffen, K. et al. Barrettides: A peptide family specifically produced by the deep-sea sponge Geodia barretti. J. Nat. Prod. 84, 3138–3146 (2021).CAS 
PubMed 
PubMed Central 

Google Scholar 
Abbamondi, G. R., De Rosa, S., Iodice, C. & Tommonaro, G. Cyclic dipeptides produced by marine sponge-associated bacteria as quorum sensing signals. Nat. Prod. Commun. 9, 229–232 (2014).CAS 
PubMed 

Google Scholar 
Kasheverov, I. et al. 6-Bromohypaphorine from Marine Nudibranch Mollusk Hermissenda crassicornis is an agonist of human α7 nicotinic acetylcholine receptor. Mar. Drugs 13, 1255–1266 (2015).CAS 
PubMed 
PubMed Central 

Google Scholar 
Moitinho-Silva, L. et al. The sponge microbiome project. Gigascience 6, 1–7 (2017).CAS 
PubMed 

Google Scholar 
Kielak, A. M., Barreto, C. C., Kowalchuk, G. A., van Veen, J. A. & Kuramae, E. E. The ecology of acidobacteria: Moving beyond genes and genomes. Front. Microbiol. 7, 744 (2016).Crits-Christoph, A., Diamond, S., Butterfield, C. N., Thomas, B. C. & Banfield, J. F. Novel soil bacteria possess diverse genes for secondary metabolite biosynthesis. Nature 558, 440–444 (2018).ADS 
CAS 
PubMed 

Google Scholar  More

To our knowledge, this study provides the first recorded information on the active swimming of Japanese eels and on their transport by currents in the open ocean. Specifically, the strong flow of the KC largely dominated the movements of the eels and transported them northeastward while they swam mainly southward, and active swimming contributed a little to their travel trajectories. In contrast, the swimming of eels made a relatively higher contribution to their travel trajectories in the TS area.Our in situ estimates of the mean swimming speeds of Japanese eels (26–41 cm/s) were similar or slightly lower than those of European eels. In the acoustic tracking experiment of European eels considering environmental current vectors, their swimming speeds were 35–58 cm/s in the coastal midwater26. In a laboratory experiment using stamina tunnels with stable temperatures, the optimal swimming speeds of European eels were estimated to be 61–68 cm/s (0.74–1.02 BL/s)56, which were higher than the in situ estimates. The minimum swimming speed of European eels is considered to be 40 cm/s if they will arrive at their spawning area in the Sargasso Sea (distance of 5500 km) in 6 months, and their optimal swimming speeds were sufficient to migrate over the long distance in time for the near-spawning period after escape from their growth habitats56. However, field studies using PSAT tagging also reported that in situ migration speeds (including transport by currents) were less than the optimal swimming speeds and suggested that some European eels could reach their spawning area within the near-spawning periods and that others only arrive in time for the following spawning season19.Our estimated effective swimming speed of Japanese eels, all day and all night over the tracking periods, ranged from 3 to 30 cm/s with individual variations. These estimates were consistent with the swimming speeds (excluding transport by currents) of 2.2–15.1 km/day (2–18 cm/s) estimated in the PSAT study of Japanese eels14. Silver-phase Japanese eels start migrating from their coastal growth habitats in Japan primarily in October to December57, 58, and spawning near the West Mariana Ridge occurs in April to August33, 35. Numerical models assuming that migrating eels use true navigation (readjusted compass) or a constant compass heading (fixed compass from the departure place to the spawning site) indicate that the minimal swimming speed required to arrive at the spawning area within 8 months is 10–12 cm/s37. Our estimated effective swimming speeds of five out of ten eels during the day and eight out of ten eels during the night were similar or higher than these minimal speeds. The low effective swimming speeds frequently observed during the day might be due to the relatively low values observed in the swimming speed at 10 min intervals and the swimming directions often varying during the day. When eels swim with stable orientation, as observed in three of the eels (WE2999_TS, WE3001_TS, and WE3002_TS) during the night, the effective swimming speeds exceeded 25 cm/s. If such a stable orientation is maintained and compensate the low speeds during the day, the eels that leave during autumn and winter will be able to arrive at the spawning area during the next spring to summer.It should also be noted that the swimming speeds in body length per second were significantly higher in shallow water during the night than in deep water during the day. In the open ocean, anguillid species exhibit DVMs during oceanic migration, swimming at depth during the day and in the shallows during the night9,10,11,12,13,14,15,16,17,18,19,20,21,22. These DVMs are likely related to the possible avoidance from visual predators under light conditions19 or maturation control59. Essentially, through the DVMs, the eels encounter low temperatures ( 20 °C) during the same day. Generally, the swimming speeds of fishes are restricted by the ambient water temperature60, and the water temperature encountered through DVMs might influence the horizontal-swimming speeds of Japanese eels.Other factors besides swimming speed are important for the success of eel migrations, such as adapting to mesopelagic zones that silver eels undergo during their spawning migrations. The most important and obvious morphological adaptation in mesopelagic fish is their well-developed eyes, and migrating eels also seem to use this strategy. These fish often have relatively large pupils61, high photosensitive structures, such as tubular eyes62, a pure rod multibank retina63, and maximum rhodopsin absorption to adapt to the blue-green light in the deep sea64. The eyes of catadromous eels displayed enlargement during their transformation into migrating silver-phase eels65, 66 and potentially increase their retinal surface area, which results in the possibility of increased photon capture. In addition, the rhodopsins in the eyes change from a freshwater type with a maximum absorption of ~ 500 nm to a deep-sea type with a maximum absorption of ~ 480 nm67,68,69. Their extreme sensitivity to light is evident through their DVM in mesopelagic water, where the timing of a large descent and ascent in the DVM demonstrated by migrating catadromous eels is precisely synchronized with sunrise and sunset. Furthermore, eels alter their swimming depth in response to the phase of the Moon9, 15, 20, 21, appearing to be capable of perceiving extremely low-intensity moonlight.This study showed that three eels released in the TS area (mainly 300–400-m depth) and one eel in the KC area (near surface) were found to change their swimming direction around the time of the solar culmination when the Sun’s bearing changed. The clockwise and counterclockwise trajectories of these eels corresponded to whether the Sun moved from the east to west in the southern and northern sky, suggesting that they demonstrated horizontal negative phototaxis swimming to avoid sunlight. They might move to avoid high-intensity sunlight horizontally, not vertically, as they gradually increase the swimming depths possibly due to acclimation to cold deep water after release. The daytime swimming depths of the eels became deeper day-by-day after their release (Fig. 4); a similar phenomenon was observed in European eels12, American eels17, and long fin eels13. Recently, Higuchi et al.20 observed that the daytime swimming depths of Japanese eels released in the TS area gradually became deeper until 13 days after their release. These facts indicate that they gradually acclimate to the cold water at the deep depths after release. Since this tracking study was conducted 2–8 days after their release, the daytime swimming depth of eels would not have reached a steady state yet. The relatively high intensity from sunlight at the shallow depths where eels swam immediately after release in the TS area might cause horizontal avoidance behavior from the light.In other cases, many eels, especially those released in the KC area, did not demonstrate the rotational behavior. The eels in the KC area mostly stayed deeper (500–800 m) during the day than the eels in the TS area (stayed at depths of 300–600 m) even during the periods shortly after their release. This is possibly due to higher water temperatures even at the deeper depths in the KC area (Fig. 4). The eels in the TS area did not demonstrate clear rotational behavior at depths of more than 400 m. The PSAT studies have reported that the steady swimming depths during the day were 500–800 m14, 20. Therefore, it was assumed that the rotational behavior observed in some eels was not a regular behavior during their migration. However, the rotational behavior observed in this study suggests that they surely perceive the horizontal direction of Sun’s bearing at 400 m depths at least. Generally, they exhibit DVM precisely synchronizing with sunrise and sunset and surely perceive the change in sunlight intensity at deeper depths9,10,11,12,13,14,15,16,17,18,19,20,21,22. Even though the rotational behavior were not observed below 400 m, it remains unknown whether the eels could not perceive the Sun’s bearing from the light penetrated at depth; thus, further investigation of response to underwater light is required in future.While possible negative phototaxis behaviors were observed in some eels after release around the time of solar culmination, the trajectories of ten eels during the entire period of tracking experiments implied that each eel tended to swim meridionally toward the bearing of the Sun at culmination. We observed that eels released at middle (20°–34° N) and low (12°–13°N) latitudes tended to swim southward and northward in the meridional direction, respectively (Fig. 6A, B). The tendency to move in a north–south swimming direction corresponded to whether the Sun culminated to the north or south: eels swam southward if the culmination occurred in the southern sky, but they swam northward if it occurred in the northern sky (Fig. 6). In the KC area (33°–35° N), the Sun rose in the southeast, passed celestial meridian in the southern sky, and set in the southwest (Fig. 6C). At 20° N in the summer time when the tracking study was conducted, the Sun also passed a celestial meridian in the southern sky, but rose in the northeast and set in the northwest (Fig. 6C). When Sun culmination occurred in the southern sky, the meridional swimming directions tended to be southward (Fig. 6A). However, at 12° to 13° N in the summer time, the Sun rose in the northeast, passed the celestial meridian in the northern sky, and set in the northwest (Fig. 6D). When the Sun at culmination appeared in the northern sky, the meridional swimming directions tended to be northward (Fig. 6B). Furthermore, the swimming behavior by one eel (WE4264_TS) that was released slightly south (14° 15′ N) from the latitude with the Sun passing through the zenith was also indicative of the meridional swimming traits. This eel moved in a northerly direction on the first day, but then it lost its north–south bias in swimming around 14° 30′ N, where the Sun nearly passed through the zenith (Figs. 1 and 6D). These observations imply that the eels might move toward the latitude with the Sun passing through the zenith.Figure 6Swimming trajectories of eels and solar paths in the celestial sphere viewed from east during each tracking period. Swimming trajectories of eels released at (A) 20°N in the tropical–subtropical area and the Kuroshio Current area, and (B) 12°–14°15′N in the tropical–subtropical area. Solar paths through the north (N)–south (S) axis and the zenith at the time of tracking in (C) 20°N in the tropical–subtropical area and the Kuroshio Current area, and (D) 12°–14°15′N in the tropical–subtropical area.Full size imageTheoretically, it is possible for mesopelagic animals to use solar cues for navigation at depths shallower than the asymptotic depth, below which penetrating light rays are symmetrical around the vertical axis and the polarization plane becomes horizontal. For example, the Sargasso Sea, where the two Atlantic catadromous eels spawn1, 3, has extremely transparent water70, and the major axis of radiance distribution still remains tilted in the mesopelagic zone. The angle of maximum radiance of sunlight at 475 nm was 13° at depths of 400 m when the Sun’s elevation was 60° (Fig. 7)52, 53. In highly transparent water, the asymptotic depth could be as high as 1000–1200 m, and the depths below this cannot be utilized for compass use53. Currently, it is not possible to verify whether the Sun culminating to north or south caused the meridional swimming tendencies of eels in this study. Potentially, these meridional swimming tendencies could be due to other orientation clues, such as the geomagnetic field, as discussed for temperate anguillid eels17, 45. Nevertheless, in future studies, it would be worthwhile considering solar cues as a possible candidate factor in the orientation of eels, even when under faint underwater light conditions.Figure 7Optical features of underwater sunlight. (A) Schematic diagram of sunlight penetrating the deep ocean at 90° to the solar bearing. The line of arrows indicates the major axis of the incident beam in a vertical plane perpendicular to the Sun’s bearing. Blue light (around 475 nm) reaches the lowest depths. With increasing depth, the light field alters its character into a less directed distribution and a lower energetic level through scattering and absorption processes. Penetrating light rays are symmetrical around the region below the asymptotic depth. (B) An example of spectral radiance distribution (e. g. 475 nm) at a certain depth. The radiance distribution is shown by an ellipsoid and the major axis is drawn by a line with arrow. The refracted angular deviation (a) of the major axis of underwater radiance distribution from the vertical axis equals the tilt of the electric vector (ee bar) from the horizontal axis53. When the Sun’ s elevation was 60° in the Sargasso Sea, the radiance distributions were measured at three different depths and the tilt of the electric vector were estimated to be 24° at depths of 100 m and 200 m and 13° at depth of 400 m52, 53.Full size imageGiven that eels might be able to use the Sun’s bearing at culmination to orient their meridional swimming direction, this orientation scheme could support a clockwise eel migration route following a partial subtropical gyre2, 37. Japanese eels that departed from the nursery area first transported northeastward via the strong KC. Maintaining southward swimming in the current, they eventually crossed the current and shifted to the southward migration course. When they enter the KC, movement to the left of the bearing of the Sun at culmination (i.e., south) is the typical pattern for the early migration of eels from Japan. The movements of eels observed in the KC were consistent with the expected route; however, eels released at low latitudes of the TS area often swam northward but also westward, which resulted in their traveling an unreasonable distance from the spawning area. This might be due to their behavior during early migration. In this study, eels were transported from Japan and released into the open ocean at low latitudes. They might have swum toward the expected bearing of the Sun at culmination as if they were in the north and moved to the left of the Sun’s bearing along with the North Equatorial Current, which would mimic the early migration of eels leaving Japan and moving along the KC.Among the eels tracked in this study were individuals with impaired swim bladders, yellow-phase eels in the process of hormone-treatment maturation, and silver-phase eels collected from different rivers in different years. Despite these variations, the swimming characteristics of the eels did not differ in terms of their DVM behavior16 and swimming speed. Nevertheless, confirmation of our results using samples with a uniform status in future research would be highly desirable. In this study, the tracked eel position was assumed to be identical to that of the tracking ship and the errors between these two positions could not be evaluated; thus, the positioning of tracked fish also may need to be improved in future studies. Experimental studies, such as tracking of blind, magnetically disturbed, or olfactory-blocked eels, could help obtain or eliminate alternative candidate clues and enhance our understanding of the navigational system of anguillid eels. Controlled laboratory experiments are required to directly quantify the ability of eels to perceive radiance distribution or polarization, along with any associated behaviors. In addition, the internal clock of eels required to perform celestial navigation should be investigated. Meanwhile, the results obtained from this study can enhance our knowledge of the mechanisms underlying the migratory behaviors of eels in the open ocean. More

Kim, I. R. et al. Genetic diversity and population structure of nutria (Myocastor coypus) in South Korea. Animals 9, 1164. https://doi.org/10.3390/ani9121164 (2019).Article 
PubMed Central 

Google Scholar 
GISD. Of the World’s Worst Invasive Alien Species. Global Invasive Species Database. http://www.iucngisd.org/gisd/100_worst.php. 100, (2021).Hong, S., Do, Y., Kim, J. Y., Kim, D. & Joo, G. Distribution, spread and habitat preferences of nutria (Myocastor coypus) invading the lower Nakdong River, South Korea. Biol. Invas. 17, 1485–1496. https://doi.org/10.1007/s10530-014-0809-8 (2015).Article 

Google Scholar 
Ojeda, R., Bidau, C. & Emmons, L. Myocastor coypus (errata version published in 2017). The IUCN Red List Threat. Species (2016): e.T14085A121734257.Tsiamis, K. et al. Baseline Distribution of Invasive Alien Species of Union Concern (Publications Office of the European Union, 2017).
Google Scholar 
Carter, J. & Leonard, B. P. A review of the literature on the worldwide distribution, spread of, and efforts to eradicate the coypu (Myocastor coypus). Wildl. Soc. Bull. 30, 162–175 (2002).
Google Scholar 
Kim, Y. C. et al. Distribution and management of nutria (Myocastor coypus) populations in South Korea. Sustainability 11, 4169. https://doi.org/10.3390/su11154169 (2019).Article 

Google Scholar 
Park, J. H. et al. The first case of Capillaria hepatica infection in a nutria (Myocastor coypus) in Korea. Korean J. Parasitol. 52, 527–529. https://doi.org/10.3347/kjp.2014.52.5.527 (2014).Article 
PubMed 
PubMed Central 

Google Scholar 
Fratini, F., Turchi, B. E., Ebani, V. V. & Bertelloni, F. The presence of Leptospira in coypus (Myocastor coypus) and rats (Rattus norvegicus) living in a protected wetland in Tuscany (Italy). Vet. Arh. 85, 407–414 (2015).
Google Scholar 
Lee, D. H., Kil, J. H. & Kim, D. E. The study on the distribution and inhabiting status of nutria (Myocastor coypus) in Korea. Korean J. Environ. Ecol. 27, 316–326 (2013).CAS 

Google Scholar 
Guichón, M. L., Doncaster, C. P. & Cassini, M. H. Population structure of coypus (Myocastor coypus) in their region of origin and comparison with introduced populations. J. Zool. 261, 265–272. https://doi.org/10.1017/S0952836903004187 (2003).Article 

Google Scholar 
Bertolino, S., Perrone, A. & Gola, L. Effectiveness of coypu control in small Italian Wetland areas. Wildl. Soc. Bull. 33, 714–720. https://doi.org/10.2193/0091-7648(2005)33[714:EOCCIS]2.0.CO;2 (2005).Article 

Google Scholar 
Schertler, A. et al. The potential current distribution of the coypu (Myocastor coypus) in Europe and climate change induced shifts in the near future. NeoBiota 58, 129–160. https://doi.org/10.3897/neobiota.58.33118 (2020).Article 

Google Scholar 
Hilts, D. J., Belitz, M. W., Gehring, T. M., Pangle, K. L. & Uzarski, D. G. Climate change and nutria range expansion in the Eastern United States. J. Wild. Manaag. 83, 591–598. https://doi.org/10.1002/jwmg.21629’ (2019).Article 

Google Scholar 
Jarnevich, C. et al. Evaluating simplistic methods to understand current distributions and forecast distribution changes under climate change scenarios: An example with coypu (Myocastor coypus). NeoBiota 32, 107–125. https://doi.org/10.3897/neobiota.32.8884 (2017).Article 

Google Scholar 
Korean Metrological Administration, (2020). Korean Climate Change Assessment Report 2020.Guillera-Arroita, G. et al. Is my species distribution model fit for purpose? Matching data and models to applications. Glob. Ecol. Biogeogr. 24, 276–292. https://doi.org/10.1111/geb.12268 (2015).Article 

Google Scholar 
Phillips, S. J., Anderson, R. P. & Schapire, R. E. Maximum entropy modeling of species geographic distributions. Ecol. Modell. 190, 231–259. https://doi.org/10.1016/j.ecolmodel.2005.03.026 (2006).Article 

Google Scholar 
Hong, S., Cowan, P., Do, Y. & Gim, J. S. Seasonal feeding habits of coypu (Myocastor coypus) in South Korea. Hystrix 27, 123–128 (2016).
Google Scholar 
Kim, H. S., Kong, J. Y., Kim, J. H., Yeon, S. C. & Hong, I. H. A Case of Fascioliasis in A Wild Nutria, Myocastor coypus Republic of Korea. Korean J. Parasitol. 56, 375–378. https://doi.org/10.3347/kjp.2018.56.4.375 (2018).Article 
PubMed 
PubMed Central 

Google Scholar 
Do, Y., Kim, J. Y., Im, R. Y. & Kim, S. B. Spatial distribution and social characteristics for wetlands in Gyeongsangnam-do Province. Korean J. Limnol. 45, 252–260 (2012).
Google Scholar 
IPCC. Climate Change 2013: The Physical Science Basis. Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. (2013).Sheffels, T. R. Status of Nutria (Myocastor coypus) Populations in the Pacific Northwest and Development of Associated Control and Management Strategies, with an Emphasis on Metropolitan Habitats, PhD Thesis (Portland State Univ., 2013).Doncaster, C. P. & MlCOL, T. Annual cycle of a coypu (Myocastor coypus) population: Male and female strategies. J. Zool. 217, 227–240. https://doi.org/10.1111/j.1469-7998.1989.tb02484.x (1989).Article 

Google Scholar 
Reggiani, G., Boitani, L. & Stefano, R. Population dynamics and regulation in the coypu Myocastor coypus in Central Italy. Ecography 18, 138–146. https://doi.org/10.1111/j.1600-0587.1995.tb00334.x (1995).Article 

Google Scholar 
Cha, Y., Cho, K. H., Lee, H., Kang, T. & Kim, J. H. The relative importance of water temperature and residence time in predicting cyanobacteria abundance in regulated rivers. Water Res. 124, 11–19. https://doi.org/10.1016/j.watres.2017.07.040 (2017).CAS 
Article 
PubMed 

Google Scholar 
Hellmann, J. J., Byers, J. E., Bierwagen, B. G. & Dukes, J. S. Five potential consequences of climate change for invasive species. Conserv. Biol. 22, 534–543. https://doi.org/10.1111/j.1523-1739.2008.00951.x (2008).Article 
PubMed 

Google Scholar 
Pereira, A. D. et al. Modeling the geographic distribution of Myocastor coypus (Mammalia, Rodentia) in Brazil: Establishing priority areas for monitoring and an alert about the risk of invasion. Stud. Neotrop. Fauna Environ. 55, 139–148. https://doi.org/10.1080/01650521.2019.1707419 (2020).Article 

Google Scholar 
Pearson, R. G. & Dawson, T. P. Predicting the impacts of climate change on the distribution of species: Are bioclimate envelope models useful?. Glob. Ecol. Biogeogr. 12, 361–371. https://doi.org/10.1046/j.1466-822X.2003.00042.x (2003).Article 

Google Scholar 
Rogers, C. E. & McCarty, J. P. Climate change and ecosystems of the mid-atlantic region. Clim. Res. 14, 235–244. https://doi.org/10.3354/cr014235 (2000).Article 

Google Scholar 
Adhikari, P. et al. Potential impact of climate change on plant invasion in the Republic of Korea. J. Ecol. Environ. 43, 36. https://doi.org/10.1186/s41610-019-0134-3 (2019).Article 

Google Scholar 
Welsch, D. J., Smart, D. L., Boyer, J. N. & Minkin, P. Forested Wetlands: Functions, Benefits and the Use of Best Management Practices (US Dept of the Interior Fish and Wildlife Service, 2021).
Google Scholar 
Borgnia, M., Galante, M. L. & Cassini, M. H. Diet of the coypu (nutria, Myocastor coypus) in agro-systems of Argentinean pampas. J. Wildl. Manag. 64, 354–361. https://doi.org/10.2307/3803233 (2000).Article 

Google Scholar 
Colares, I. G., Oliveira, R. N. V., Liveira, R. M. & Colares, E. P. Feeding habits of coypu (Myocastor coypus Molina 1978) in the wetlands of the Southern region of Brazil. An. Acad. Bras. Cienc. 82, 671–678. https://doi.org/10.1590/s0001-37652010000300015 (2010).Article 
PubMed 

Google Scholar 
Corriale, M. J., Arias, S. M., Bó, R. F. & Porini, G. Habitat-use patterns of the coypu (Myocastor coypus) in an urban wetland of its original distribution. Acta Theriol. 51, 295–302. https://doi.org/10.1007/BF03192681 (2006).Article 

Google Scholar 
Linscombe, G., Kinler, N. & Wright, V. Nutria population density and vegetative changes in brackish marsh in coastal Louisiana. In Worldwide Furbearer Conference Proceedings (eds Chapman, J. A. & Pursley, D.) 129–141 (Worlwide Furbearer Conference Inc, 1981).
Google Scholar 
Aliev, F. Contribution to the study of nutria migrations (Myocastor coypus). Saugetierkd. Mitt. 16, 301–303 (1968).
Google Scholar 
Farashi, A. & Najafabadi, M. S. A model to predict dispersion of the alien nutria, Myocastor coypus Molina, 1782 (Rodentia) Northern Iran. Acta Zool. Bulg. 69, 65–70 (2017).
Google Scholar 
Vilà, M. et al. How well do we understand the impacts of alien species on ecosystem services? A pan-European, cross-taxa assessment. Front. Ecol. Environ. 8, 135–144. https://doi.org/10.1890/080083 (2010).Article 

Google Scholar 
Adhikari, P. et al. Seasonal and altitudinal variation in roe deer (Capreolus pygargus tianschanicus) diet on Jeju Island, South Korea. J. Asia Pac. Biodivers. 9, 422–428. https://doi.org/10.1016/j.japb.2016.09.001 (2016).Article 

Google Scholar 
Koo, K. A., Kong, W. S., Nibbelink, N. P., Hopkinson, C. S. & Lee, J. H. Potential effects of climate change on the distribution of cold-tolerant evergreen broadleaved woody plants in the Korean Peninsula. PLoS ONE 10, e0134043. https://doi.org/10.1371/journal.pone.0134043 (2015).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
National Institute of Biological Research. Korean Red List of Threatened Species 2nd edn. (Ministry of Environement of Korea, 2014).
Google Scholar 
Kil, J. et al. Monitoring of Invasive Alien Species Designated by the Wildlife Protection Act (VII) (Natl Inst. of Environmental Research, 2013).
Google Scholar 
Busby, J. R. In Bioclim, a Bioclimatic Analysis and Prediction System in Nature Conservation: Cost Effective Biological Surveys and Data Analysis (eds Margules, C. R. & Austin, M. P.) 64–68 (CSIRO, 1991).
Google Scholar 
Lee I. H., Park S. H., Kang, H. S. & Cho C. H. Regional climate projections using the HadGEM3-RA in Proceedings of the 3rd International Conference on Earth System Modelling; Hamburg, Germany. 17–21 September 2012. (2012).Robert, J. H., Phillips, S., Leathwick, J. & Elith, J. Package ‘dismo’ version 1.3. , https://cran.rproject.org/web/packages/dismo.pdf (2021).Jeon, J. Y., Adhikari, P. & Seo, C. Impact of climate change on potential dispersal of Paeonia obovata (Paeoniaceae), a critically endangered medicinal plant of South Korea. Ecol. Environ. Conserv. 26, S145–S155 (2020).
Google Scholar 
Dormann, C. F. et al. Collinearity: A review of methods to deal with it and a simulation study evaluating their performance. Ecography 36, 27–46. https://doi.org/10.1111/j.1600-0587.2012.07348.x (2013).Article 

Google Scholar 
Shin, M. S., Seo, C., Lee, M. & Kim, J. Y. Prediction of potential species richness of plants adaptable to climate change in the Korean Peninsula. J. Environ. Impact Assess. 27, 562–581 (2018).
Google Scholar 
Adhikari, P. et al. Northward range expansion of southern butterflies according to climate change in South Korea. KSCCR 11, 643–656. https://doi.org/10.15531/KSCCR.2020.11.6.643 (2020).Article 

Google Scholar 
Song, C. et al. Estimation of future land cover considering shared socioeconomic pathways using scenario generators. KSCCR 9, 223–234. https://doi.org/10.15531/KSCCR.2018.9.3.223 (2018).Article 

Google Scholar 
Fick, S. E. & Hijmans, R. J. WorldClim 2: New 1-km spatial resolution climate surfaces for global land areas. Int. J. Climatol. 37, 4302–4315. https://doi.org/10.1002/joc.5086 (2017).Article 

Google Scholar 
Dukes, J. S. & Mooney, H. A. Does global change increase the success of biological invaders?. Trends Ecol. Evol. 14, 135–139. https://doi.org/10.1016/s0169-5347(98)01554-7 (1999).CAS 
Article 
PubMed 

Google Scholar 
Thuiller, W., Georges, D., Gueguen, M., Engler, R. & Breiner, F. Package ‘biomod2’: Ensemble Platform for Species Distribution Modeling, version 3.5.1. https://cran.r-project.org/web/packages/biomod2/biomod2.pdf (2021).Elith, J. et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography 29, 129–151. https://doi.org/10.1111/j.2006.0906-7590.04596.x (2006).Article 

Google Scholar 
Barbet-Massin, M., Jiguet, F., Albert, C. H. & Thuiller, W. Selecting pseudo-absences for species distribution models: how, where and how many?. Methods Ecol. Evol. 3, 327–338. https://doi.org/10.1111/j.2041-210X.2011.00172.x (2012).Article 

Google Scholar 
Brown, J. L. SDM toolbox: A python-based GIS toolkit for landscape genetic, biogeographic and species distribution model analyses. Methods Ecol. Evol. 5, 694–700. https://doi.org/10.1111/2041-210X.12200 (2014).Article 

Google Scholar 
Veloz, S. D. Spatially autocorrelated sampling falsely inflates measures of accuracy for presence–only niche models. J. Biogeogr. 36, 2290–2299. https://doi.org/10.1111/j.1365-2699.2009.02174.x (2009).Article 

Google Scholar 
Adhikari, P., Lee, Y. H., Park, Y.-S. & Hong, S. H. Assessment of the spatial invasion risk of intentionally introduced alien plant species (IIAPS) under environmental change in South Korea. Biology 10, 1169 (2021).Article 

Google Scholar 
Hong, S. H., Lee, Y. H., Lee, G., Lee, D. H. & Adhikari, P. Predicting impacts of climate change on northward range expansion of invasive weeds in South Korea. Plants 10, 1604. https://doi.org/10.3390/plants10081604 (2021).Article 
PubMed 
PubMed Central 

Google Scholar 
Pearsons, R. G. Species distribution modeling for conservation educators and practitioners. Lessons Conserv. 3, 54–58 (2010).
Google Scholar 
Allouche, O., Tsoar, A. & Kadmon, R. Assessing the accuracy of species distribution models: Prevalence, kappa and the true skill statistic (TSS). J. Appl. Ecol. 43, 1223–1232. https://doi.org/10.1111/j.1365-2664.2006.01214.x (2006).Article 

Google Scholar 
Thuiller, W., Lavorel, S. & Araújo, M. B. Niche properties and geographical extent as predictors of species sensitivity to climate change. Glob. Ecol. Biogeogr. 14, 347–357. https://doi.org/10.1111/j.1466-822X.2005.00162.x (2005).Article 

Google Scholar 
Lobo, J. M., Jiménez-Valverde, A. & Real, R. AUC: A misleading measure of the performance of predictive distribution models. Global. Ecol. Biogeography. 17, 145–151. https://doi.org/10.1111/j.1466-8238.2007.00358.x (2008).Article 

Google Scholar 
Swets, J. A. Measuring the accuracy of diagnostic systems. Science 240, 1285–1293 (1988).MathSciNet 
CAS 
Article 
ADS 

Google Scholar 
Baldwin, R. Use of maximum entropy modeling in wildlife research. Entropy 11, 854–866. https://doi.org/10.3390/e11040854 (2009).Article 
ADS 

Google Scholar 
Adhikari, P. et al. Potential impact of climate change on the species richness of subalpine plant species in the mountain national parks of South Korea. J. Ecol. Environ. 42, 36. https://doi.org/10.1186/s41610-018-0095-y (2018).Article 

Google Scholar 
Hijmans, R. J. et al. Package ‘raster’ v 3.5: geographical data analysis and modeling. https://cran.r-project.org/web/packages/raster/raster.pdf, (2021). More

We used publicly available, daily, county-level COVID-19 cases and deaths from the Pennsylvania Department of Health (PA DOH) (https://www.health.pa.gov/topics/disease/coronavirus/pages/cases.aspx)13,14 for Centre County and the six neighboring counties with which it shares borders: Blair, Clearfield, Clinton, Huntingdon, Mifflin, and Union (Table 1, Fig. 1). Official COVID-19 reporting for these counties began on March 1, 2020 and is ongoing.Table 1 Summary statistics. COVID-19 reporting, census data, SafeGraph mobile-device derived data.Full size tableFigure 1(a) The cumulative COVID-19 case trajectory for Centre County minus the student cases (red line) has the same shape as the outbreak for the neighboring counties. When looking at student cases only (blue line), the curve leads other counties. Centre County cumulative cases including the university (purple line) take on the shape of an early increase because of the student cases. (b) When aggregating cases from students and non-students, Centre County (purple dot) reported about the number of cases expected for its population size, relative to the neighboring counties (black dots). When the university-reported student cases are separated from the non-student residents of the county, cases reported in Centre County non-students (red dots show possible range of total cases) fall below the number of cases we would expect for the population size. Student cases only (blue dot) are high for the student population size.Full size imageWithin Centre County, PSU provided COVID-19 testing for UP students from August 7, 2020 onward and reported anonymized weekly (2020) and daily (2021) confirmed cases, negative test results, and total tests completed for each campus in a public dashboard (Figs. 1a, S1) (https://virusinfo.psu.edu/covid-19-dashboard/)8. Two types of testing were conducted: students who were enrolled in on-campus classes were randomly selected for surveillance testing and all students could use on-demand testing. Through March 23, 2021, a total of 45,092 random tests were conducted for surveillance, of which 462, or 1.0%, were infected. Surveillance testing efforts ranged from 2440 to 4020 weekly tests through the Fall 2020 semester and were designed to consistently test approximately 1% of students throughout the school year.During the same time period, 75,436 on-demand tests were conducted, of which 6093, or 8.1%, were infected. Students living in both on-campus dorms and off-campus apartments had equal access to university-provided testing. Both on-campus and off-campus residences are within Centre County so positive and negative tests results were also included in the overall Centre County reports of COVID-19 cases.Pre-arrival testing was required for students returning to campus from transmission hotspots. Students with positive tests from pre-arrival testing were required to isolate for 10–14 days after their positive test before arriving on campus. Results from pre-arrival testing for students returning to campus in the Fall of 2020 are not included in these data.At the county level, PA DOH reports the total positive, probable, and negative tests for each county. Because PSU is within Centre County, we estimated the number of total positive and negative tests for non-student Centre County residents by subtracting the PSU estimates (from the PSU dashboard) from the Centre County estimates provided by PA DOH. However, not all student tests were reported to DOH. A portion of the on-demand tests conducted for PSU UP students were completed by a third-party vendor, which required student registration. At the time of student registration, an estimated 0–25% of students registered with an address for a family home that did not reflect their residence in Centre County. Their test results were reported to the county of their registered address. This impacts a maximum of 1,166 positive student test results and 10,760 negative student tests.We conducted a sensitivity analysis to assess the uncertainty in reporting around the negative and positive students tests that may have been misallocated due to the reported residence of student tests. We have calculated the minimum and maximum number of affected positive and negative student tests. This uncertainty from student tests impacts non-student values, which are calculated by subtracting student values from county level reports. The calculations are based on a range of a possible 0–1166 positive student tests misallocated to other counties and up to 10,760 misallocated negative student tests. We have used the ranges of misallocated student tests to calculate, for non-student Centre County residents, the full possible range of (1) total cases, (2) reported cases per capita, and (3) tests per capita (Table 1, Fig. 1b). As a result, our estimates of cases and per capita testing among non-student residents in Centre County are imprecise (Table 1).We also used publicly available data from PA DOH data and PSU to calculate COVID-19 deaths per 100,000 for Centre County, the six neighboring counties, and PSU UP.We acquired county-level data on median household income, population size, and college enrollment status from the 2019 United States Census Bureau’s American Community Survey (ACS) 5-year data (https://www.census.gov/data/developers/data-sets/acs-5year.html) for all seven previously mentioned counties in central PA15.We divide the census block groups (CBG) of Centre County into two categories. We first designated ‘student-dominated CBGs’ as CBGs where  > 50% of ACS responses report enrollment as undergraduate students. We consider data from the 19 student-dominated CBGs in Centre County to be representative of the student population in Centre County. In addition to off-campus locations, the 19 student-dominated CBGs include all on-campus dorms. These 19 CBGs are either on or adjacent to PSU’s UP campus and occupy exactly 6 census tracts. The remaining 25 county census tracts were designated as non-student dominated areas.SafeGraph16 receives geolocation data from anonymized mobile devices collected from numerous applications. We analyzed SafeGraph’s mobile device-derived daily visit counts to points of interest (POI), which are fixed locations, such as businesses or attractions. SafeGraph data provide daily counts for total numbers of visits by mobile devices while using at least one application that provides geolocation data to SafeGraph. A “visit” indicates that the device entered the building or spatial perimeter designated as a POI. We acquired daily visit counts for POIs in the seven previously mentioned counties in central PA from January 1, 2019 forward (Table 1) and within Centre County grouped counts into student-dominated CBGs and non-student dominated CBGs. From January 1, 2020 forward, we used SafeGraph data on the median daily minutes that devices spent outside of their home in each county and the student- and non-student dominated CBG divisions in Centre County. The “home location” of each device is defined by its location overnight. Finally, we used SafeGraph’s weekly calculated number of devices residing in each county and the CBGs of Centre County for 2019 to measure SafeGraph’s data representation across the seven counties and the CBGs of Centre County.No administrative permissions were required to obtain these data. Academic researchers can register to receive access to SafeGraph data at no charge for non-commercial purposes only. See Data Availability statement below for details. More

Knapp, A. K. & Smith, M. D. Variation among biomes in temporal dynamics of aboveground primary production. Science 291, 481–484 (2001).CAS 
PubMed 
ADS 

Google Scholar 
Collins, S. L. et al. Stability of tallgrass prairie during a 19-year increase in growing season precipitation. Funct. Ecol. 26, 1450–1459 (2012).
Google Scholar 
Maurer, G. E., Hallmark, A. J., Brown, R. F., Sala, O. E. & Collins, S. L. Sensitivity of primary production to precipitation across the United States. Ecol. Lett. 23, 527–536 (2020).PubMed 

Google Scholar 
IPCC. IPCC. (Cambridge University Press, 2013) https://doi.org/10.1017/cbo9781107415324.Knapp, A. K. et al. Differential sensitivity to regional-scale drought in six central US grasslands. Oecologia 177, 949–957 (2015).PubMed 
ADS 

Google Scholar 
Smith, M. D. An ecological perspective on extreme climatic events: A synthetic definition and framework to guide future research. J. Ecol. 99, 656–663 (2011).
Google Scholar 
Zeppel, M. J. B., Wilks, J. V. & Lewis, J. D. Impacts of extreme precipitation and seasonal changes in precipitation on plants. Biogeosciences 11, 3083–3093 (2014).ADS 

Google Scholar 
Frank, D. A. Drought effects on above- and belowground production of a grazed temperate grassland ecosystem. Oecologia 152, 131–139 (2007).PubMed 
ADS 

Google Scholar 
Skinner, R. H., Hanson, J. D., Hutchinson, G. L. & Schuman, G. E. Response of C3 and C4 grasses to supplemental summer precipitation. J. Range Manag. 55, 517–522 (2002).
Google Scholar 
Shi, Z. et al. Dual mechanisms regulate ecosystem stability under decade-long warming and hay harvest. Nat. Commun. 7, 1–6 (2016).ADS 

Google Scholar 
Zavaleta, E. S. et al. Grassland responses to three years of elevated temperature, CO2, precipitation, and N deposition. Ecol. Monogr. 73, 585–604 (2003).
Google Scholar 
Prather, R. M., Castillioni, K., Welti, E. A. R., Kaspari, M. & Souza, L. Abiotic factors and plant biomass, not plant diversity, strongly shape grassland arthropods under drought conditions. Ecology 101, 1–7 (2020).
Google Scholar 
Nippert, J. B., Knapp, A. K. & Briggs, J. M. Intra-annual rainfall variability and grassland productivity: Can the past predict the future?. Plant Ecol. 184, 65–74 (2006).
Google Scholar 
La Pierre, K. J. et al. Explaining temporal variation in above-ground productivity in a mesic grassland: The role of climate and flowering. J. Ecol. 99, 1250–1262 (2011).
Google Scholar 
Cleland, E. E. et al. Sensitivity of grassland plant community composition to spatial vs. temporal variation in precipitation. Ecology 94, 1687–1696 (2013).PubMed 

Google Scholar 
Grant, K., Kreyling, J., Heilmeier, H., Beierkuhnlein, C. & Jentsch, A. Extreme weather events and plant–plant interactions: Shifts between competition and facilitation among grassland species in the face of drought and heavy rainfall. Ecol. Res. 29, 991–1001 (2014).
Google Scholar 
Brooker, R. W. et al. Facilitation in plant communities: The past, the present, and the future. J. Ecol. 96, 18–34 (2008).MathSciNet 

Google Scholar 
Schöb, C., Armas, C. & Pugnaire, F. I. Direct and indirect interactions co-determine species composition in nurse plant systems. Oikos 122, 1371–1379 (2013).
Google Scholar 
Gross, N., Börger, L., Duncan, R. P. & Hulme, P. E. Functional differences between alien and native species: Do biotic interactions determine the functional structure of highly invaded grasslands?. Funct. Ecol. 27, 1262–1272 (2013).
Google Scholar 
van der Merwe, S., Greve, M., Olivier, B. & le Roux, P. C. Testing the role of functional trait expression in plant–plant facilitation. Funct. Ecol. https://doi.org/10.1111/1365-2435.13681 (2020).Article 

Google Scholar 
Tremmel, D. C. & Bazzaz, F. A. How neighbor canopy architecture affects target plant performance. Ecology 74, 2114–2124 (1993).
Google Scholar 
Weiher, E. & Keddy, P. A. In Ecological Assembly Rules: Perspective, Advances, Retreats. (eds. Weiher, E. & Keddy, P. A.) (Cambridge University Press, 2001).Anten, N. P. R. & Hirose, T. Interspecific differences in above-ground growth patterns result in spatial and temporal partitioning of light among species in a tall-grass meadow. J. Ecol. 87, 583–597 (1999).
Google Scholar 
Yann, H., Pascal, A. & Niklaus, A. H. Competition for light causes plant. Science 324, 636–638 (2009).
Google Scholar 
Walker, B., Kinzig, A. & Langridge, J. Plant attribute diversity, resilience, and ecosystem function: The nature and significance of dominant and minor species. Ecosystems 2, 95–113 (1999).
Google Scholar 
Brooker, R. W. Plant–plant interactions and environmental change. New Phytol. 171, 271–284 (2006).PubMed 

Google Scholar 
Michalet, R. & Pugnaire, F. I. Facilitation in communities: Underlying mechanisms, community and ecosystem implications. Funct. Ecol. 30, 3–9 (2016).
Google Scholar 
Maestre, F. T., Callaway, R. M., Valladares, F. & Lortie, C. J. Refining the stress-gradient hypothesis for competition and facilitation in plant communities. J. Ecol. 97, 199–205 (2009).
Google Scholar 
Saccone, P., Delzon, S., Jean-Philippe, P., Brun, J. J. & Michalet, R. The role of biotic interactions in altering tree seedling responses to an extreme climatic event. J. Veg. Sci. 20, 403–414 (2009).
Google Scholar 
Smith, M. D., Knapp, A. K. & Collins, S. L. A framework for assessing ecosystem dynamics in response to chronic resource alterations induced by global change. Ecology 90, 3279–3289 (2009).PubMed 

Google Scholar 
Borer, E. T., Seabloom, E. W., Gruner, D. S., Harpole, W. S. & Hillebrand, H. Herbivores and nutrients control grassland plant diversity via light limitation. Nature 508, 517–520 (2014).CAS 
PubMed 
PubMed Central 
ADS 

Google Scholar 
de Sassi, C. & Tylianakis, J. M. Climate change disproportionately increases herbivore over plant or parasitoid biomass. PLoS One 7, e40557 (2012).PubMed 
PubMed Central 

Google Scholar 
Strauss, S. Y. & Ivalú Cacho, N. Nowhere to run, nowhere to hide: The importance of enemies and apparency in adaptation to harsh soil environments. Am. Nat. 182, E1 (2013).PubMed 

Google Scholar 
Brady, K. U., Kruckeberg, A. R. & Bradshaw, H. D. Evolutionary ecology of plant adaptation to serpentine soils. Annu. Rev. Ecol. Evol. Syst. 36, 243–266 (2005).
Google Scholar 
Moran, M. S. et al. Soil evaporation response to Lehmann lovegrass (Eragrostis lehmanniana) invasion in a semiarid watershed. Agric. For. Meteorol. 149, 2133–2142 (2009).ADS 

Google Scholar 
Pérez-Harguindeguy, N. et al. New handbook for standardised measurement of plant functional traits worldwide. Aust. J. Bot. 61, 167–234 (2013).
Google Scholar 
Díaz, S. et al. The global spectrum of plant form and function. Nature 529, 167–171 (2016).PubMed 
ADS 

Google Scholar 
Gross, N., Suding, K. N. & Lavorel, S. Leaf dry matter content and lateral spread predict response to land use change for six subalpine grassland species. J. Veg. Sci. 18, 289–300 (2007).
Google Scholar 
Quiroga, R., Golluscio, R., Blanco, L. & Fernandez, R. Aridity and grazing as convergent selective forces: An experiment with an Arid Chaco bunchgrass. Ecol. Appl. https://doi.org/10.1890/09-0641 (2010).Article 
PubMed 

Google Scholar 
Blumenthal, D. M. et al. Traits link drought resistance with herbivore defence and plant economics in semi-arid grasslands: The central roles of phenology and leaf dry matter content. J. Ecol. 108, 2336–2351 (2020).
Google Scholar 
Taylor, S. H. et al. Ecophysiological traits in C3 and C4 grasses: A phylogenetically controlled screening experiment. New Phytol. 185, 780–791 (2010).CAS 
PubMed 

Google Scholar 
N’Guessan, M. & Hartnett, D. C. Differential responses to defoliation frequency in little bluestem (Schizachyrium scoparium) in tallgrass prairie: Implications for herbivory tolerance and avoidance. Plant Ecol. 212, 1275–1285 (2011).
Google Scholar 
Castillioni, K. et al. Drought mildly reduces plant dominance in a temperate prairie ecosystem across years. Ecol. Evol. 10, 6702–6713 (2020).PubMed 
PubMed Central 

Google Scholar 
Ivalú Cacho, N. & Strauss, S. Y. Occupation of bare habitats, an evolutionary precursor to soil specialization in plants. Proc. Natl. Acad. Sci. U. S. A. 111, 15132–15137 (2014).PubMed 
PubMed Central 
ADS 

Google Scholar 
Cottingham, K. L., Lennon, J. T. & Brown, B. L. Knowing when to draw the line: Designing more informative ecological experiments. Front. Ecol. Environ. 3, 145–152 (2005).
Google Scholar 
Xu, X., Sherry, R. A., Niu, S., Li, D. & Luo, Y. Net primary productivity and rain-use efficiency as affected by warming, altered precipitation, and clipping in a mixed-grass prairie. Glob. Change Biol. 19, 2753–2764 (2013).ADS 

Google Scholar 
Braun-Blanquet, J. Plant Sociology: The Study of Plant Communities. (1932).Shipley, B. The AIC model selection method applied to path analytic models compared using ad-separation test. Ecology 94, 560–564 (2013).PubMed 

Google Scholar 
Lefcheck, J. S. piecewiseSEM: Piecewise structural equation modelling in r for ecology, evolution, and systematics. Methods Ecol. Evol. 7, 573–579 (2016).
Google Scholar 
Grace, J. B. In Structural Equation Modeling and Natural Systems. (Cambridge University Press, 2006). https://doi.org/10.1017/CBO9780511617799.Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D. & Team, R.C. nlme: Linear and nonlinear mixed effects models. R package version 3.1 111 (2013).Bates, D., Mächler, M., Bolker, B. M. & Walker, S. C. Fitting linear mixed-effects models using lme4. Ecol. Austral 67, 1–48 (2015).
Google Scholar 
Pearson, D. E., Ortega, Y. K. & Maron, J. L. The tortoise and the hare: reducing resource availability shifts competitive balance between plant species. J. Ecol. 105, 999–1009 (2017).CAS 

Google Scholar 
Maron, J. L. & Crone, E. Herbivory: Effects on plant abundance, distribution and population growth. Proc. R. Soc. B Biol. Sci. 273, 2575–2584 (2006).
Google Scholar 
Bertness, M. & Callaway, R. M. Positive interactions in communities. Trends Ecol. Evol. 9, 191–193 (1994).CAS 
PubMed 

Google Scholar 
Ploughe, L. W. et al. Community Response to Extreme Drought (CRED): A framework for drought-induced shifts in plant–plant interactions. New Phytol. 222, 52–69 (2019).PubMed 

Google Scholar 
Klanderud, K., Vandvik, V. & Goldberg, D. The importance of Biotic vs. Abiotic drivers of local plant community composition along regional bioclimatic gradients. PLoS One 10, 1–15 (2015).
Google Scholar 
Maricle, B. R., Caudle, K. L. & Adler, P. B. Influence of water availability on photosynthesis, water potential, leaf δ 13 C, and phenology in dominant C 4 grasses in Kansas, USA. Trans. Kans. Acad. Sci. 118, 173–193 (2015).
Google Scholar 
Collins, S. L., Knapp, A. K., Briggs, J. M., Blair, J. M. & Steinauer, E. M. Modulation of diversity by grazing and mowing in native tallgrass prairie. Science 280, 745–747 (1998).CAS 
PubMed 
ADS 

Google Scholar 
Gornish, E. S. & Tylianakis, J. Community shifts under climate change: Mechanisms at multiple scales. Am. J. Bot. 100, 1422–1434 (2013).PubMed 

Google Scholar  More

Download PDF

I study the Chinese giant salamander (Andrias davidianus), which is native to the Yangtze River Basin of central China. This particular species is critically endangered in the wild owing to habitat loss and overcatching — a particular problem is their use in traditional Chinese medicine. My research focuses on the salamander’s conservation biology and evolutionary ecology.In this photo, I am releasing a Chinese giant salamander at the Golden Whip River in Zhangjiajie National Forest Park on an early morning in September 2021. My team and I caught the salamander the night before, to measure its size and collect tissue samples for genetic analyses.My interest in aquatic animals started as a child. I grew up in a rural village in Hunan province, and I remember spending most of my childhood playing and fishing near my home. Because of this, I knew where each fish species lived in nearby rivers and lakes, and it sparked my interest in river ecology.I’m employed as an associate professor at Jishou University, where I lead a team dedicated to researching this species of salamander. Wild salamanders are quiet, nocturnal animals that live in remote areas. This makes studying them challenging. My team tried many creative ways to track down the animals, including walking along riverbanks with torches and photographing salamanders under water — but these techniques didn’t work as well as we needed them to. We eventually found that the best way to trap wild salamanders is to use small live fish and chicken livers as bait. The research is challenging, but we’ve learnt to be patient and celebrate every small success we have.Studying Chinese giant salamanders has also taught me an important life lesson: adapt to thrive. When food is abundant, the salamanders grow rapidly; when food is scarce, they can go up to 11 months without feeding. In my personal life and work, I have experienced successes and failures, and taking on that lesson has been useful.

Nature 603, 194 (2022)
doi: https://doi.org/10.1038/d41586-022-00564-y

Related Articles

Close-up with a parasite that can blind

Handling snakes for science

Broaden your scientific audience with video animation

Managing up: how to communicate effectively with your PhD adviser

Subjects

Careers

Conservation biology

Ecology

Latest on:

Careers

Smaller science company? Tailor your CV for a manager, not HR
Career Column 25 FEB 22

Female scientists in Africa are changing the face of their continent
Editorial 22 FEB 22

African scientists engage with the public to tackle local challenges
Career Feature 15 FEB 22

Ecology

How colonialism fed the flames of Australia’s catastrophic wildfires
Research Highlight 24 FEB 22

Apply Singapore Index on Cities’ Biodiversity at scale
Correspondence 22 FEB 22

Marching in the streets for climate-crisis action
Career Q&A 22 FEB 22

Jobs

POST-DOC POSITION IN ELECTROPHYSIOLOGY OF FUNGAL NETWORKS

VU University Amsterdam
Amsterdam, Netherlands

Director, Division of Receipt and Referral Center for Scientific Review National Institutes of Health (NIH) Department of Health and Human Services (DHHS)

National Institutes of Health (NIH)
Bethesda, MD, United States

Postdoctoral Fellow

NIH National Heart, Lung, and Blood Institute (NHLBI)
Bethesda, MD, United States

Two postdoctorial researchers in structural biology, with focus on artificial intelligence

University of Gothenburg (GU)
Uppsala, Sweden More

The full-length sequences of PacBio SMRT sequencingBased on PacBio SMRT sequencing, 3,751,089, 3,434,452, 3,900,180, 8,535,019, and 4,435,846 subreads were generated for root, stem, leaf, flower, and seed, with a N50 of 3040, 3367, 2611, 2198, and 4584 bp, respectively (Table S1; Fig. S1). Subreads were processed to generate circular consensus sequences (CCSs). By detecting the primers and poly(A) tail, 238,196, 232,290, 211,535, 257,905, and 231,877 full-length non-chimeric (FLNC) reads were identified for root, stem, leaf, flower, and seed, with a mean length of 2633, 3070, 2561, 1746, and 3762 bp, respectively (Table S2; Fig. S2). After Iterative Clustering for Error Correction (ICE) clustering, polishing, base correction, de-redundancy, and non-plant sequences filtering, 37,789, 34,034, 38,100, 54,937, and 53,906 unigenes were retained for root, stem, leaf, flower, and seed, respectively, with an average unigene length of 1802–3786 bp and N50 of 2238–4707 bp (Table S2). The length of most unigenes from five organs exceeded 2000 bp, accounting for 68.88% of the total number (Table S3; Fig. 1A). Based on Benchmarking Universal Single-Copy Orthologs (BUSCO) assessment, about 88.1% (single-copy: 353; duplicated: 916) of the 1440 core embryophyte genes were found to be complete (90.6% were present when counting fragmented genes), suggesting the high integrity of the M. micrantha transcriptome (Fig. S3).Figure 1Length distribution of unigenes from PacBio SMRT sequencing (A) and Illumina RNA-Seq (B) across five organs.Full size imageDe novo assembly of Illumina RNA-Seq dataBased on Illumina RNA-Seq, 43.23, 40.27, 41.01, 65.85, and 41.09 million clean reads were obtained for root, stem, leaf, flower, and seed, respectively, with Q20 exceeding 96.72%. Using Trinity software, clean reads were de novo assembled into 124,238, 60,232, 63,370, 93,229, and 66,411 unigenes for root, stem, leaf, flower, and seed. After filtering non-plant sequences, 124,233, 60,232, 63,370, 93,228, and 66,410 unigenes were finally retained for the five organs, respectively (Table S4). The length of most unigenes (84.70%) was shorter than 2000 bp (Table S3). In addition, the average length and N50 of unigenes generated by Illumina RNA-Seq were 1067–1312 bp and 1336–1685 bp, respectively, which were shorter than that from PacBio SMRT sequencing (Table S4; Fig. 1B).Functional annotationTo obtain a comprehensive functional annotation of M. micrantha transcriptome, unigenes generated by PacBio SMRT sequencing were annotated in seven public databases, including NCBI non-redundant nucleotide sequences (NT), NCBI non-redundant protein sequences (NR), Gene Ontology (GO), Eukaryotic Orthologous Groups (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG), Swiss-Prot, and Pfam protein families. For root, stem, leaf, flower, and seed, 35,714 (94.51%), 32,614 (95.83%), 36,134 (94.84%), 49,197 (89.55%), and 50,962 (94.54%) unigenes were annotated to at least one database, respectively, suggesting that our transcriptome is well annotated and that most of unigenes may be functional (Table 1).Table 1 Statistics of annotation of full-length transcripts from five M. micrantha organs in seven databases.Full size tableBased on NR database annotation, the top three homologous species for the five organs were Cynara cardunculus, Vitis vinifera, and Daucus carota (Fig. S4). The top homologous species was a plant of the Asteraceae family. For the GO function annotation, “binding”, “catalytic activities”, “metabolic process”, “cellular process”, “cell”, and “cell part” were functional categories with the most abundant unigenes (Fig. S5). In addition, numerous unigenes were assigned to “response to stimulus”, “response to biotic stimulus”, and “response to oxidative stress” category (Table S5). Positive response to stress stimuli is an important strategy for invasive plants to adapt to the environment. In the KEGG annotation, the top two pathways with the most abundant unigenes were “carbohydrate metabolism” and “translation”. Furthermore, “energy metabolism” and “environmental adaptation” were also worthy of attention, which are important pathways responsible for energy supply and stress responses (Fig. S6).TFs identification and AS analysisUsing the iTAK pipeline, 1776 (root), 1293 (stem), 1627 (leaf), 2529 (flower), and 1733 (seed) unigenes were identified as TFs, which were classified into 68 families (Table S6). C3H (884), C2H2 (525), and bHLH (501) were the most abundant TF families (Fig. S7A). In addition, MYB (333) TFs were also found in the five organs. The differential expression levels of the top 15 TF families were further characterized. We found that the top 15 TF families had a certain amount of expression in the five organs of M. micrantha (Fig. S7B).For root, stem, leaf, flower, and seed, 3300, 2324, 3219, 4730, and 3740 unique transcript models (UniTransModels) were constructed, among which the UniTransModels containing two isoforms were the most common (Fig. S8A). There were 329, 270, 358, 336, and 537 AS events identified in root, stem, leaf, flower, and seed, respectively. Retained introns (RIs) were detected as the most abundant AS event in all five organs, followed by alternative 3′ splice sites (A3) and alternative 5′ splice sites (A5). Mutually exclusive exons (MX) were the least frequent event (Fig. S8B).Gene expression analysisThe number of unigenes in different expression level intervals was similar across the five organs (Fig. 2A). Using FPKM  > 0.3 as the threshold for unigene expression, the total number of unigenes expressed in the five organs was 102,464 (Fig. 2B). Among them, 39,227 unigenes were co-expressed in all five organs. The information of differentially expressed genes (DEGs) identified in pairwise comparisons among the five organs is listed in Table S7. In total, 21,161 DEGs were identified among the five organs (Fig. S9). The number of DEGs between the five organs varied from 3469 (root vs stem) to 10,716 (leaf vs seed) (Fig. 2C). Notably, 933, 428, 1410, 1018, and 1292 DEGs showed significant higher expression in root, stem, leaf, flower, and seed, respectively (Figs. S10 and S11).Figure 2Gene expression patterns in five M. micrantha organs. (A) The FPKM interval distribution in the five organs. (B) Venn diagram of the number of unigenes expressed in five organs. (C) Number of differentially expressed genes in each pairwise comparison of five organs.Full size imageKEGG enrichment of unigenes with higher expression in each organAccording to the KEGG enrichment analysis results, there were obvious differences in enriched pathways in the five organs (Table S8; Fig. 3). The unigenes with higher expression in root were mainly enriched to defense response and protein processing pathways, such as “plant-pathogen interaction” and “protein processing in endoplasmic reticulum”. In stem, unigenes with higher expression were predominantly enriched to pathways related to the secondary metabolite, sugar, and terpenoid biosynthesis, such as “phenylpropanoid biosynthesis”, “starch and sucrose metabolism”, and “diterpenoid biosynthesis”. In flower, unigenes with higher expression were mainly related to “starch and sucrose metabolism”, “phenylpropanoid biosynthesis”, and “cutin, suberine, and wax biosynthesis”. The unigenes with higher expression in seed were mainly enriched in three fatty acid and sugar metabolism pathways, namely “biosynthesis of unsaturated fatty acids”, “galactose metabolism”, and “amino sugar and nucleotide sugar metabolism”. The unigenes with higher expression in leaf were significantly enriched in photosynthesis pathways, including “photosynthesis-antenna proteins”, “photosynthesis”, “porphyrin and chlorophyll metabolism”, and “carbon fixation in photosynthetic organisms”, which are important for the photosynthesis of M. micrantha.Figure 3The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis of unigenes with higher expression in each organ. The significantly enriched pathways with corrected p-value (q value)  More