Ecology

1.Batjes, N. H. Total carbon and nitrogen in the soils of the world. Eur. J. Soil Sci. 65, 10–21 (1996).Article 
CAS 

Google Scholar 
2.Lal, R. Soil carbon sequestration impacts on global climate change and food security. Science 304, 1623–1627 (2004).CAS 
Article 

Google Scholar 
3.Buringh, P. in The role of terrestrial vegetation in the global carbon cycle: Measurement by remote sensing, 91–109 (Wiley, 1984).4.Hiederer, R. & Köchy, M. Global soil organic carbon estimates and the harmonized world soil database. EUR 79, 25225 (2011).
Google Scholar 
5.Smith, P. et al. Global change pressures on soils from land use and management. Glob. Chang. Biol. 22, 1008–1028 (2016).Article 

Google Scholar 
6.Schlesinger, W. H. The Role of Terrestrial Vegetation in the Global Carbon Cycle: Measurement by Remote Sensing (Wiley, 1984).7.Conant, R. T., Cerri, C. E., Osborne, B. B. & Paustian, K. Grassland management impacts on soil carbon stocks: a new synthesis. Ecol. Appl. 27, 662–668 (2017).Article 

Google Scholar 
8.Köchy, M., Hiederer, R. & Freibauer, A. Global distribution of soil organic carbon–Part 1: masses and frequency distributions of SOC stocks for the tropics, permafrost regions, wetlands, and the world. Soil 1, 351–365 (2015).Article 
CAS 

Google Scholar 
9.Nahlik, A. M. & Fennessy, M. S. Carbon storage in US wetlands. Nat. Commun. 7, 1–9 (2016).Article 
CAS 

Google Scholar 
10.Dixon, R. K. et al. Carbon pools and flux of global forest ecosystems. Science 263, 185–190 (1994).CAS 
Article 

Google Scholar 
11.Atwood, T. B. et al. Global patterns in mangrove soil carbon stocks and losses. Nat. Clim. Chang. 7, 523–528 (2017).CAS 
Article 

Google Scholar 
12.Maclean, J. L., Dawe, D. C., Hardy, B. & Hettel, G. P. Rice Almanac: Source book for the most important economic activity on earth, 3rd edn. (CABI Publishing, 2002).13.Kögel-Knabner, I. et al. Biogeochemistry of paddy soils. Geoderma 157, 1–14 (2010).Article 
CAS 

Google Scholar 
14.Wu, J. Carbon accumulation in paddy ecosystems in subtropical China: evidence from landscape studies. Eur. J. Soil Sci. 62, 29–34 (2011).CAS 
Article 

Google Scholar 
15.Carlson, K. M. et al. Greenhouse gas emissions intensity of global croplands. Nat. Clim. Chang. 7, 63–68 (2017).CAS 
Article 

Google Scholar 
16.FAO (Food and Agriculture Organization of the United Nations). FAOSTAT: FAO Statistical Databases. http://faostat.fao.org/default.aspx (2018).17.Gattinger, A. et al. Enhanced top soil carbon stocks under organic farming. Proc. Natl Acad. Sci. USA 109, 18226–18231 (2012).CAS 
Article 

Google Scholar 
18.Xie, Z. et al. Soil organic carbon stocks in China and changes from 1980s to 2000s. Glob. Chang. Biol. 13, 1989–2007 (2007).Article 

Google Scholar 
19.Qin, Z., Huang, Y. & Zhuang, Q. Soil organic carbon sequestration potential of cropland in China. Glob. Biogeochem. Cycles 27, 711–722 (2013).CAS 
Article 

Google Scholar 
20.Jobbágy, E. G. & Jackson, R. B. The vertical distribution of soil organic carbon and its relation to climate and vegetation. Ecol. Appl. 10, 423–436 (2000).Article 

Google Scholar 
21.Haefele, S. M., Nelson, A. & Hijmans, R. J. Soil quality and constraints in global rice production. Geoderma 235, 250–259 (2014).Article 
CAS 

Google Scholar 
22.Pan, G., Li, L., Wu, L. & Zhang, X. Storage and sequestration potential of topsoil organic carbon in China’s paddy soils. Glob. Chang. Biol. 10, 79–92 (2004).Article 

Google Scholar 
23.Wei, L. et al. Comparing carbon and nitrogen stocks in paddy and upland soils: Accumulation, stabilization mechanisms, and environmental drivers. Geoderma 398, 115121 (2021).Article 

Google Scholar 
24.Wang, P. et al. Long-term rice cultivation stabilizes soil organic carbon and promotes soil microbial activity in a salt marsh derived soil chronosequence. Sci. Rep. 5, 15704 (2015).CAS 
Article 

Google Scholar 
25.Li, Y. et al. Oxygen availability determines key regulators in soil organic carbon mineralisation in paddy soils. Soil Biol. Biochem. 153, 108106 (2021).CAS 
Article 

Google Scholar 
26.Evans, C. D. et al. Acidity controls on dissolved organic carbon mobility in organic soils. Glob. Chang. Biol. 18, 3317–3331 (2012).Article 

Google Scholar 
27.Liu, Y. et al. Impact of prolonged rice cultivation on coupling relationship among C, Fe, and Fe-reducing bacteria over a 1000-year paddy soil chronosequence. Biol. Fertil. Soils 55, 589–602 (2019).CAS 
Article 

Google Scholar 
28.Sinsabaugh, R. L. et al. Stoichiometry of soil enzyme activity at global scale. Ecol. Lett. 11, 1252–1264 (2008).Article 

Google Scholar 
29.Liu, Y. et al. Microbial activity promoted with organic carbon accumulation in macroaggregates of paddy soils under long-term rice cultivation. Biogeosciences 13, 6565–6586 (2016).CAS 
Article 

Google Scholar 
30.Liu, Y. et al. Methanogenic abundance and changes in community structure along a rice soil chronosequence from east China. Eur. J. Soil Sci. 67, 443–455 (2016).CAS 
Article 

Google Scholar 
31.Malik, A. A. et al. Land use driven change in soil pH affects microbial carbon cycling processes. Nat. Commun. 9, 1–10 (2018).CAS 
Article 

Google Scholar 
32.Don, A., Schumacher, J. & Freibauer, A. Impact of tropical land‐use change on soil organic carbon stocks-a meta‐analysis. Glob. Chang. Biol. 17, 1658–1670 (2011).Article 

Google Scholar 
33.Piao, S. et al. The carbon balance of terrestrial ecosystems in China. Nature 458, 1009–1013 (2009).CAS 
Article 

Google Scholar 
34.Davidson, E. A. & Janssens, I. A. Temperature sensitivity of soil carbon decomposition and feedbacks to climate change. Nature 440, 165–173 (2006).CAS 
Article 

Google Scholar 
35.Kirk, G. The Biogeochemistry of Submerged Soils (Wiley, 2004).36.Kramer, M. G., Sanderman, J., Chadwick, O. A., Chorover, J. & Vitousek, P. M. Long‐term carbon storage through retention of dissolved aromatic acids by reactive particles in soil. Glob. Chang. Biol. 18, 2594–2605 (2012).Article 

Google Scholar 
37.Scharpenseel, H. W., Pfeiffer, E. M. & Becker-Heidmann, P. in Advances in Soil Science (eds. Carter, MR, Stewart, BA) (Lewis Publishers, 1996).38.Liao, Q. et al. Increase in soil organic carbon stock over the last two decades in China’s Jiangsu Province. Glob. Chang. Biol. 15, 861–875 (2009).Article 

Google Scholar 
39.Keiluweit, M., Wanzek, T., Kleber, M., Nico, P. & Fendorf, S. Anaerobic microsites have an unaccounted role in soil carbon stabilization. Nat. Commun. 8, 1–10 (2017).CAS 
Article 

Google Scholar 
40.Ghimire, R., Lamichhane, S., Acharya, B. S., Bista, P. & Sainju, U. M. Tillage, crop residue, and nutrient management effects on soil organic carbon in rice-based cropping systems: a review. J. Integr. Agric. 16, 1–15 (2017).Article 

Google Scholar 
41.Maillard, É. & Angers, D. A. Animal manure application and soil organic carbon stocks: a meta‐analysis. Glob. Chang. Biol. 20, 666–679 (2014).Article 

Google Scholar 
42.Tian, K. et al. Effects of long-term fertilization and residue management on soil organic carbon changes in paddy soils of China: a meta-analysis. Agric. Ecosyst. Environ. 204, 40–50 (2015).CAS 
Article 

Google Scholar 
43.Liu, Y. et al. Initial utilization of rhizodeposits with rice growth in paddy soils: rhizosphere and N fertilization effects. Geoderma 338, 30–39 (2019).CAS 
Article 

Google Scholar 
44.Chen, J. et al. A keystone microbial enzyme for nitrogen control of soil carbon storage. Sci. Adv. 4, eaaq1689 (2018).CAS 
Article 

Google Scholar 
45.Zhu, Z. et al. Rice rhizodeposits affect organic matter decomposition in paddy soil: the role of N fertilization and rice growth for enzyme activities, CO2 and CH4 emissions. Soil Biol. Biochem. 116, 369–377 (2018).CAS 
Article 

Google Scholar 
46.Moorhead, D. L. & Sinsabaugh, R. L. A theoretical model of litter decay and microbial interaction. Ecol. Monogr. 76, 151–174 (2006).Article 

Google Scholar 
47.Li, X. et al. Nitrogen fertilization decreases the decomposition of soil organic matter and plant residues in planted soils. Soil Biol. Biochem. 112, 47–55 (2017).CAS 
Article 

Google Scholar 
48.Cui, J. et al. Carbon and nitrogen recycling from microbial necromass to cope with C:N stoichiometric imbalance by priming. Soil Biol. Biochem. 142, 107720 (2020).CAS 
Article 

Google Scholar 
49.Geisseler, D., Linquist, B. A. & Lazicki, P. A. Effect of fertilization on soil microorganisms in paddy rice systems—a meta-analysis. Soil Biol. Biochem. 115, 452–460 (2017).CAS 
Article 

Google Scholar 
50.Sun, W. et al. Climate drives global soil carbon sequestration and crop yield changes under conservation agriculture. Glob. Chang. Biol. 26, 3325–3335 (2020).Article 

Google Scholar 
51.Wissing, L. et al. Management-induced organic carbon accumulation in paddy soils: the role of organo-mineral associations. Soil Tillage Res. 126, 60–71 (2013).Article 

Google Scholar 
52.Baker, J. M., Ochsner, T. E., Venterea, R. T. & Griffis, T. J. Tillage and soil carbon sequestration—-what do we really know? Agric. Ecosyst. Environ. 118, 1–5 (2007).CAS 
Article 

Google Scholar 
53.Lal, R. Challenges and opportunities in soil organic matter research. Eur. J. Soil Sci. 60, 158–169 (2009).CAS 
Article 

Google Scholar 
54.Lal, R. Soil carbon sequestration in India. Clim. Change 65, 277–296 (2004).CAS 
Article 

Google Scholar 
55.Liu, Y. et al. Carbon input and allocation by rice into paddy soils: a review. Soil Biol. Biochem. 133, 97–107 (2019).CAS 
Article 

Google Scholar 
56.Zhao, Y. et al. Economics-and policy-driven organic carbon input enhancement dominates soil organic carbon accumulation in Chinese croplands. Proc. Natl Acad. Sci. USA 115, 4045–4050 (2018).CAS 
Article 

Google Scholar 
57.Wei, X., Zhu, Z., Wei, L., Wu, J. & Ge, T. Biogeochemical cycles of key elements in the paddy-rice rhizosphere: microbial mechanisms and coupling processes. Rhizosphere 10, 100145 (2019).Article 

Google Scholar 
58.Alexandratos, N. & Bruinsma, J. World agriculture towards 2030/2050: the 2012 revision. https://doi.org/10.22004/ag.econ.288998. (2012).59.Rui, W. & Zhang, W. Effect size and duration of recommended management practices on carbon sequestration in paddy field in Yangtze Delta Plain of China: a meta-analysis. Agric. Ecosyst. Environ. 135, 199–205 (2010).CAS 
Article 

Google Scholar 
60.Song, K. et al. Wetland degradation: its driving forces and environmental impacts in the Sanjiang Plain, China. Environ. Manage. 54, 255–271 (2014).Article 

Google Scholar 
61.Dong, J. et al. Northward expansion of paddy rice in northeastern Asia during 2000–2014. Geophys. Res. Lett. 43, 3754–3761 (2016).CAS 
Article 

Google Scholar 
62.Chaturvedi, V. et al. Climate mitigation policy implications for global irrigation water demand. Mitig. Adapt. Strat. Glob. Chang. 20, 389–407 (2015).Article 

Google Scholar 
63.Gathorne-Hardy, A. A life cycle assessment (LCA) of greenhouse gas emissions from SRI and flooded rice production in SE India. Taiwan Water Conserv. J. 61, 111–125 (2013).
Google Scholar 
64.Linquist, B., Van Groenigen, K. J., Adviento‐Borbe, M. A., Pittelkow, C. & Van Kessel, C. An agronomic assessment of greenhouse gas emissions from major cereal crops. Glob. Chang. Biol. 18, 194–209 (2012).Article 

Google Scholar 
65.IPCC. in Contribution of working group II to the fifth assessment report of the Intergovernmental Panel on Climate Change. (eds. Field, C. B. et al) (Cambridge University Press, 2014).66.Xie, Z. et al. CO2 mitigation potential in farmland of China by altering current organic matter amendment pattern. Sci. China Earth Sci. 53, 1351–1357 (2010).CAS 
Article 

Google Scholar 
67.Yan, X. et al. Carbon sequestration efficiency in paddy soil and upland soil under long-term fertilization in southern China. Soil Tillage Res. 130, 42–51 (2013).Article 

Google Scholar 
68.Shang, Q. et al. Net annual global warming potential and greenhouse gas intensity in Chinese double rice‐cropping systems: a 3‐year field measurement in long‐term fertilizer experiments. Glob. Chang. Biol. 17, 2196–2210 (2011).Article 

Google Scholar 
69.Ma, Y. et al. Net global warming potential and greenhouse gas intensity of annual rice–wheat rotations with integrated soil–crop system management. Agric. Ecosyst. Environ. 164, 209–219 (2013).Article 

Google Scholar 
70.Xiong, Z. et al. Differences in net global warming potential and greenhouse gas intensity between major rice-based cropping systems in China. Sci. Rep. 5, 1–9 (2015).CAS 

Google Scholar 
71.Jiang, Y. et al. Acclimation of methane emissions from rice paddy fields to straw addition. Sci. Adv. 5, eaau9038 (2019).Article 
CAS 

Google Scholar 
72.Liu, C., Lu, M., Cui, J., Li, B. & Fang, C. Effects of straw carbon input on carbon dynamics in agricultural soils: a meta‐analysis. Glob. Chang. Biol. 20, 1366–1381 (2014).Article 

Google Scholar 
73.Shakoor, A. et al. A global meta-analysis of greenhouse gases emission and crop yield under no-tillage as compared to conventional tillage. Sci. Total Environ. 750, 142299 (2021).CAS 
Article 

Google Scholar 
74.Zhao, X. et al. Methane and nitrous oxide emissions under no‐till farming in China: a meta‐analysis. Glob. Chang. Biol. 22, 1372–1384 (2016).Article 

Google Scholar 
75.Kim, S. Y., Gutierrez, J. & Kim, P. J. Unexpected stimulation of CH4 emissions under continuous no-tillage system in mono-rice paddy soils during cultivation. Geoderma 267, 34–40 (2016).CAS 
Article 

Google Scholar 
76.Ball, B. C., Scott, A. & Parker, J. P. Field N2O, CO2 and CH4 fluxes in relation to tillage, compaction and soil quality in Scotland. Soil Tillage Res. 53, 29–39 (1999).Article 

Google Scholar 
77.Linquist, B. A., Adviento-Borbe, M. A., Pittelkow, C. M., van Kessel, C. & van Groenigen, K. J. Fertilizer management practices and greenhouse gas emissions from rice systems: a quantitative review and analysis. Field Crop. Res. 135, 10–21 (2012).Article 

Google Scholar 
78.Schlesinger, W. H. Carbon sequestration in soils: some cautions amidst optimism. Agric. Ecosyst. Environ. 82, 121–127 (2000).CAS 
Article 

Google Scholar 
79.Choudhury, A. T. M. A. & Kennedy, I. R. Nitrogen fertilizer losses from rice soils and control of environmental pollution problems. Commun. Soil Sci. Plan. 36, 1625–1639 (2005).CAS 
Article 

Google Scholar 
80.Jiang, Y. et al. Water management to mitigate the global warming potential of rice systems: a global meta-analysis. Field Crop. Res. 234, 47–54 (2019).Article 

Google Scholar 
81.Suryavanshi, P., Singh, Y. V., Prasanna, R., Bhatia, A. & Shivay, Y. S. Pattern of methane emission and water productivity under different methods of rice crop establishment. Paddy Water Environ. 11, 321–329 (2013).Article 

Google Scholar 
82.Yan, X., Akiyama, H., Yagi, K. & Akimoto, H. Global estimations of the inventory and mitigation potential of methane emissions from rice cultivation conducted using the 2006 Intergovernmental Panel on Climate Change Guidelines. Glob. Biogeochem. Cycles https://doi.org/10.1029/2008GB003299 (2009).83.Jiang, Y. et al. Higher yields and lower methane emissions with new rice cultivars. Glob. Chang. Biol. 23, 4728–4738 (2017).Article 

Google Scholar 
84.Li, C. et al. Modeling greenhouse gas emissions from rice-based production systems: sensitivity and upscaling. Glob. Biogeochem. Cycles https://doi.org/10.1029/2003GB002045 (2004).85.Yin, S. et al. Carbon sequestration and emissions mitigation in paddy fields based on the DNDC model: a review. Artif. Intell. Agric. 4, 140–149 (2020).
Google Scholar 
86.FAO, IIASA, ISRIC, ISSCAS, and JRC: Harmonized World Soil Database (version 1.2), Tech. Rep., FAO, Rome, Italy and IIASA, Laxenburg, Austria (2012).87.Allison, L. in Organic carbon. Methods of Soil Analysis: Part 2 Chemical and Microbiological Properties, (ed. A.g. Norman). (American Society of Agronomy, 1965).88.Fang, C. & Moncrieff, J. B. The variation of soil microbial respiration with depth in relation to soil carbon composition. Plant Soil 268, 243–253 (2005).CAS 
Article 

Google Scholar 
89.Yan, X., Cai, Z., Wang, S. & Smith, P. Direct measurement of soil organic carbon content change in the croplands of China. Glob. Chang. Biol. 17, 1487–1496 (2011).Article 

Google Scholar 
90.Rosenberg, M. S., Adams, D. C. & Gurevitch, J. MetaWin 2.0: statistical software for meta-analysis (Sinauer, 2000).91.Yue, Q. et al. Deriving emission factors and estimating direct nitrous oxide emissions for crop cultivation in China. Environ. Sci. Technol. 53, 10246–10257 (2019).CAS 
Article 

Google Scholar 
92.Hedges, L. V., Gurevitch, J. & Curtis, P. S. The meta‐analysis of response ratios in experimental ecology. Ecology 80, 1150–1156 (1999).Article 

Google Scholar 
93.Adams, D. C., Gurevitch, J. & Rosenberg, M. S. Resampling tests for meta‐analysis of ecological data. Ecology 78, 1277–1283 (1997).Article 

Google Scholar 
94.Van Groenigen, K. J., Osenberg, C. W. & Hungate, B. A. Increased soil emissions of potent greenhouse gases under increased atmospheric CO2. Nature 475, 214–216 (2011).Article 
CAS 

Google Scholar  More

1.Ghosh, S. K., Podder, D., Panja, S., & Mukherjee, S. In target areas where human mosquito-borne diseases are diagnosed, the inclusion of the pre-adult mosquito aquatic niches parameters will improve the integrated mosquito control program. PLos Neg. Trop. Dis. 14(8), e0008605 (2020).Article 

Google Scholar 
2.Becker, B. N. et al. Mosquitoes and Their Control 499 (Springer, 2010).Book 

Google Scholar 
3.Hyde, K. D. et al. The amazing potential of fungi: 50 ways we can exploit fungi industrially. Fungal Divers. 97, 1–136 (2019).Article 

Google Scholar 
4.Clark, T. B., Kellen, W. R., Fukuda, T. & Lindegren, J. E. Field and laboratory studies on the pathogenicity of the fungus Beauveria bassiana to three genera of mosquitoes. J. Invertebr. Pathol. 11(1), 1–7 (1968).PubMed 
Article 
PubMed Central 

Google Scholar 
5.Scholte, E. J., Knols, B. G. & Takken, W. Infection of the malaria mosquito Anopheles gambiae with the entomopathogenic fungus Metarhizium anisopliae reduces blood feeding and fecundity. J. Invertebr. Pathol. 91(1), 43–49 (2006).PubMed 
Article 
PubMed Central 

Google Scholar 
6.Bukhari, T., Takken, W. & Koenraadt, C. J. Development of Metarhizium anisopliae and Beauveria bassiana formulations for control of malaria mosquito larvae. Parasit. Vectors 4(1), 23 (2011).PubMed 
PubMed Central 
Article 

Google Scholar 
7.Mukherjee, A., Debnath, P., Ghosh, S. K. & Medda, P. K. Biological control of papaya aphid (Aphis gossypii Glover) using entomopathogenic fungi. Vegetos 33, 1–10 (2020).Article 

Google Scholar 
8.Fernández-Grandon, G. M., Harte, S. J., Ewany, J., Bray, D. & Stevenson, P. C. Additive effect of botanical insecticide and entomopathogenic fungi on pest mortality and the behavioral response of its natural enemy. Plants 9, 173 (2020).PubMed Central 
Article 
CAS 

Google Scholar 
9.Sobczak, J. F. et al. Manipulation of wasp (Hymenoptera: Vespidae) behavior by the entomopathogenic fungus Ophiocordyceps humbertii in the Atlantic forest in Ceará, Brazil. Entomol. News 129, 98–104 (2020).Article 

Google Scholar 
10.Ghosh, S. K. & Pal, S. Entomopathogenic potential of Trichoderma longibrachiatum and its comparative evaluation with malathion against the insect pest Leucinodes orbonalis. Environ. Monit. Assess. 188(1), 37 (2016).PubMed 
Article 
CAS 
PubMed Central 

Google Scholar 
11.Podder, D. & Ghosh, S. K. A new application of Trichoderma asperellum as an anopheline larvicide for eco friendly management in medical science. Sci. Reps. 9(1), 1108 (2019).ADS 
Article 
CAS 

Google Scholar 
12.Jones, E. B. G. Fungal adhesion. Mycol. Res. 98(9), 961–981 (1994).Article 

Google Scholar 
13.Shah, P. A. & Pell, J. K. Entomopathogenic fungi as biological control agents. Appl. Microbiol. Biotechnol. 61, 413–423 (2003).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
14.Rudall, K. M. The chitin/protein complexes of insect cuticles. Adv. Insect Physiol. 1, 257–313 (1963).ADS 
CAS 
Article 

Google Scholar 
15.Shah, F. A., Wang, C. S. & Butt, T. M. Nutrition influences growth and virulence of the insect-pathogenic fungus Metarhizium anisopliae. FEMS Microbiol. Lett. 251(2), 259–266 (2005).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
16.Jackson, M. A., Dunlap, C. A. & Jaronski, S. T. Ecological considerations in producing and formulating fungal entomopathogens for use in insect biocontrol. Biocontrol 55(1), 129–145 (2010).Article 

Google Scholar 
17.Vega, F.E.; Meyling, N., Luangsa-ard, J.& Blackwell, M. Fungal entomopathogens. In: edit Vega, F. and Kaya, H. A. Insect pathology, 2nd edn , San Diego, CA, Academic Press, pp 171–220 (2012).18.Gaugler, R. Entomopathogenic nematodes in biological control. CRC press (2018).19.McKinnon, A. C. et al. Detection of the entomopathogenic fungus Beauveria bassiana in the rhizosphere of wound-stressed zea mays plants. Front. Microbiol. 9, 1161 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
20.Zimmermann, G. Review on safety of the entomopathogenic fungus Metarhizium anisopliae. Biocontrol Sci. Technol. 17(9), 879–920 (2007).Article 

Google Scholar 
21.Hamer, J. E., Howard, R. J., Chumley, F. G. & Valent, B. A mechanism for surface attachment in spores of a plant pathogenic fungus. Science 239(4837), 288–290 (1988).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
22.Dhawan, M. & Joshi, N. (Enzymatic comparison and mortality of Beauveria bassiana against cabbage caterpillar Pieris brassicae LINN. Braz. J. Microbiol. 48(3), 522–529 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
23.Mora, M. A. E., Castilho, A. M. C. & Fraga, M. E. Classification and infection mechanism of entomopathogenic fungi. Arq. Inst. Biol. 84, 0552015 (2017).
Google Scholar 
24.Li, J., Tracy, J. W. & Christensen, B. M. Phenol oxidase activity in hemolymph compartments of Aedes aegypti during melanotic encapsulation reactions against microfilariae. Dev. Comp. Immunol. 16(1), 41–48 (1992).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
25.Hillyer, J. F. & Strand, M. R. Mosquito hemocyte-mediated immune responses. Curr. Opin. Insect Sci. 3, 14–21 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
26.Nanda, K. P. Chronic lead (Pb) exposure results in diminished hemocyte count and increased susceptibility to bacterial infection in Drosophila melanogaster. Chemosphere 236, 124349 (2019).ADS 
CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
27.Ghosh, S. K., Chatterjee, T., Chakravarty, A. & Basak, A. K. Sodium and potassium nitrite-induced developmental genotoxicity in Drosophila melanogaster—effects in larval immune and brain stem cells. Interdiscip. Toxicol. 13(4), 101–105 (2020).
Google Scholar 
28.Chatterjee, T., Ghosh, S. K., Paik, S., Chakravarty, A. & Basak, A. K. Benzoic acid treated Drosophila melanogaster the genetic disruption of larval brain stem cells and non-neural cells during metamorphosis. Toxicol. Environ. Health Sci. https://doi.org/10.1007/s13530-021-00082-w (2021).Article 

Google Scholar 
29.Campos, R. A. Boophilus microplus infection by Beauveria amorpha and Beauveria bassiana: SEM analysis and regulation of subtilisin-like proteases and chitinases. Curr. Microbiol. 50(5), 257–261 (2005).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
30.McFarlane, H. E., Gendre, D. & Western, T. L. Seed coat ruthenium red staining assay. Bio-Protoc. 4, 1096 (2014).Article 

Google Scholar 
31.Bhosale, R. R., Osmani, R. A. M. & Moin, A. Natural gums and mucilages: A review on multifaceted excipients in pharmaceutical science and research. Int. J. Res. Phytochem. Pharmacol 6(4), 901–912 (2014).
Google Scholar 
32.Shah, F. A., Allen, N., Wright, C. J. & Butt, T. M. Repeated in vitro subculturing alters spore surface properties and virulence of Metarhizium anisopliae. FEMS Microbiol. Lett. 276(1), 60–66 (2007).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
33.Hsu, S. C. & Lockwood, J. L. Powdered chitin agar as a selective medium for enumeration of actinomycetes in water and soil. Appl. Environ. Microbiol. 29(3), 422–426 (1975).CAS 
Article 

Google Scholar 
34.Parida, D., Jena, S. K. & Rath, C. C. Enzyme activities of bacterial isolates from iron mine areas of Barbil, Keonjhar district, Odisha, India. Int. J. Pure Appl. Biosci. 2(3), 265–271 (2014).
Google Scholar 
35.Kasana, R. C., Salwan, R., Dhar, H., Dutt, S. & Gulati, A. A rapid and easy method for the detection of microbial cellulases on agar plates using Gram’s iodine. Curr. Microbiol. 57(5), 503–507 (2008).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
36.Medina, P. & Baresi, L. Rapid identification of gelatin and casein hydrolysis using TCA. J. Microbiol. Methods 69(2), 391–393 (2007).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
37.Al-Nahdi, H. S. Isolation and screening of extracellular proteases produced by new isolated Bacillus sp. J. Appl. Pharm. Sci. 2(9), 71–74 (2012).CAS 

Google Scholar 
38.Murthy, N. K. & Bleakley, B. H. Simplified method of preparing colloidal chitin used for screening of chitinase-producing microorganisms. Int. J. Microbiol. 10(2), 1937–8289 (2012).
Google Scholar 
39.Park, S. H., Lee, J. H. & Lee, H. K. Purification and characterization of chitinase from a marine bacterium, Vibrio sp. 98CJ11027. J. Microbiol 38, 224–229 (2000).CAS 

Google Scholar 
40.Roberts, W. K. & Selitrennikoff, C. P. Plant and bacterial chitinases differ in antifungal activity. Microbiology 134(1), 169–176 (1986).Article 

Google Scholar 
41.Tsuchida, O. et al. An alkaline proteinase of an alkalophilic Bacillus sp. Curr. Microbiol. 14(1), 7–12 (1986).CAS 
Article 

Google Scholar 
42.Crowell, A. M., Wall, M. J. & Doucette, A. A Maximizing recovery of water-soluble proteins through acetone precipitation. Anal. Chim. Acta. 796, 48–54 (2013).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
43.He, F. BCA (Bicinchoninic Acid) protein assay. Bio Protocol 1(5), 44 (2011).Article 

Google Scholar 
44.Sierra, L.M., Carmona, E.R., Aguado, L. & Marcos, R. The comet assay in Drosophila: neuroblast and hemocyte cells. In Genotoxicity and DNA Repair. Methods in Pharmacology and Toxicology. Humana Press, New York, NY. 269–82 (2014).45.Xu, T. et al. (2012) HMGB in mollusk Crassostrea ariakensis Gould: structure, pro-inflammatory cytokine function characterization and anti-infection role of its antibody. PLoS ONE 7(11), e50789 (2012).ADS 
CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
46.Basak, A. K., Chatterjee, T., Chakravarty, A. & Ghosh, S. K. Silver nanoparticle-induced developmental inhibition of Drosophila melanogaster accompanies disruption of genetic material of larval neural stem cells and non-neuronal cells. Environ. Monit. Assess. 191(8), 497 (2019).PubMed 
Article 
CAS 
PubMed Central 

Google Scholar  More

1.McFall-Ngai, M. et al. Animals in a bacterial world: a new imperative for the life sciences. Proc. Natl Acad. Sci. 110, 3229–3236 (2013).CAS 
Article 

Google Scholar 
2.Pita, L., Rix, L., Slaby, B. M., Franke, A. & Hentschel, U. The sponge holobiont in a changing ocean: from microbes to ecosystems. Microbiome 6, 46 (2018).CAS 
Article 

Google Scholar 
3.Caruso, R., Lo, B. C. & Núñez, G. Host–microbiota interactions in inflammatory bowel disease. Nat. Rev. Immunol. 20, 411–426 (2020).CAS 
Article 

Google Scholar 
4.Hughes, T. P. et al. Spatial and temporal patterns of mass bleaching of corals in the Anthropocene. Science 359, 80–83 (2018).CAS 
Article 

Google Scholar 
5.Scheele, B. C. et al. Amphibian fungal panzootic causes catastrophic and ongoing loss of biodiversity. Science 363, 1459–1463 (2019).CAS 
Article 

Google Scholar 
6.Bell, J. J., Bennett, H. M., Rovellini, A. & Webster, N. S. Sponges to be winners under near-future climate scenarios. BioScience 68, 955–968 (2018).Article 

Google Scholar 
7.Bosch, T. C. G., Guillemin, K. & McFall-Ngai, M. Evolutionary “experiments” in symbiosis: the study of model animals provides insights into the mechanisms underlying the diversity of host–microbe interactions. Bioessays 41, e1800256 (2019).8.Nyholm, S. V. & McFall-Ngai, M. J. A lasting symbiosis: how the Hawaiian bobtail squid finds and keeps its bioluminescent bacterial partner. Nat. Rev. Microbiol. https://doi.org/10.1038/s41579-021-00567-y (2021).Article 
PubMed 

Google Scholar 
9.Visick, K. L., Stabb, E. V. & Ruby, E. G. A lasting symbiosis: how Vibrio fischeri finds a squid partner and persists within its natural host. Nat. Rev. Microbiol. https://doi.org/10.1038/s41579-021-00557-0 (2021).Article 
PubMed 

Google Scholar 
10.Peixoto, R. S. et al. Coral probiotics: premise, promise, prospects. Annu. Rev. Anim. Biosci. 9, 265–288 (2021).Article 

Google Scholar  More

1.Serebrovsky, A. S. On the possibility of a new method for the control of insect pests. Zool. Zh. 19, 618–630 (1940).
Google Scholar 
2.Curtis, C. F. Possible use of translocations to fix desirable genes in insect pest populations. Nature 218, 368–369 (1968). This paper is one of the first to describe how reciprocal chromosomal translocations could be used to drive a favoured linked trait in a threshold-dependent fashion.CAS 
PubMed 
Article 

Google Scholar 
3.Dawkins, R. The Selfish Gene Vol. 345 (Oxford University Press, 1976).4.Bastide, H. et al. Rapid rise and fall of selfish sex-ratio X chromosomes in Drosophila simulans: spatiotemporal analysis of phenotypic and molecular data. Mol. Biol. Evol. 28, 2461–2470 (2011).CAS 
PubMed 
Article 

Google Scholar 
5.Corbett-Detig, R., Medina, P., Frerot, H., Blassiau, C. & Castric, V. Bulk pollen sequencing reveals rapid evolution of segregation distortion in the male germline of Arabidopsis hybrids. Evol. Lett. 3, 93–103 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
6.Kingan, S. B., Garrigan, D. & Hartl, D. L. Recurrent selection on the Winters sex-ratio genes in Drosophila simulans. Genetics 184, 253–265 (2010).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
7.McLaughlin, R. N. Jr. & Malik, H. S. Genetic conflicts: the usual suspects and beyond. J. Exp. Biol. 220, 6–17 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
8.Presgraves, D. C., Gerard, P. R., Cherukuri, A. & Lyttle, T. W. Large-scale selective sweep among segregation distorter chromosomes in African populations of Drosophila melanogaster. PLoS Genet. 5, e1000463 (2009).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
9.Seymour, D. K., Chae, E., Arioz, B. I., Koenig, D. & Weigel, D. Transmission ratio distortion is frequent in Arabidopsis thaliana controlled crosses. Heredity 122, 294–304 (2019).CAS 
PubMed 
Article 

Google Scholar 
10.Courret, C., Chang, C. H., Wei, K. H., Montchamp-Moreau, C. & Larracuente, A. M. Meiotic drive mechanisms: lessons from Drosophila. Proc. Biol. Sci. 286, 20191430 (2019).CAS 
PubMed 
PubMed Central 

Google Scholar 
11.Kusano, A., Staber, C., Chan, H. Y. & Ganetzky, B. Closing the (Ran)GAP on segregation distortion in Drosophila. Bioessays 25, 108–115 (2003).CAS 
PubMed 
Article 

Google Scholar 
12.Merel, V., Boulesteix, M., Fablet, M. & Vieira, C. Transposable elements in Drosophila. Mob. DNA 11, 23 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
13.Boulesteix, M. & Biemont, C. Transposable elements in mosquitoes. Cytogenet. Genome Res. 110, 500–509 (2005).CAS 
PubMed 
Article 

Google Scholar 
14.Lee, Y. C. & Langley, C. H. Transposable elements in natural populations of Drosophila melanogaster. Philos. Trans. R. Soc. Lond. B Biol. Sci. 365, 1219–1228 (2010).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
15.Kelleher, E. S. Reexamining the P-element invasion of Drosophila melanogaster through the lens of piRNA silencing. Genetics 203, 1513–1531 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
16.Majumdar, S. & Rio, D. C. P transposable elements in drosophila and other eukaryotic organisms. Microbiol. Spectr. 3, MDNA3–0004-2014 (2015).PubMed 
PubMed Central 

Google Scholar 
17.Burns, K. H. & Boeke, J. D. Human transposon tectonics. Cell 149, 740–752 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
18.Doring, H. P., Tillmann, E. & Starlinger, P. DNA sequence of the maize transposable element Dissociation. Nature 307, 127–130 (1984).CAS 
PubMed 
Article 

Google Scholar 
19.Wallau, G. L., Capy, P., Loreto, E. & Hua-Van, A. Genomic landscape and evolutionary dynamics of mariner transposable elements within the Drosophila genus. BMC Genomics 15, 727 (2014).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
20.Hawkins, J. S., Hu, G., Rapp, R. A., Grafenberg, J. L. & Wendel, J. F. Phylogenetic determination of the pace of transposable element proliferation in plants: copia and LINE-like elements in Gossypium. Genome 51, 11–18 (2008).CAS 
PubMed 
Article 

Google Scholar 
21.Biemont, C., Vieira, C., Borie, N. & Lepetit, D. Transposable elements and genome evolution: the case of Drosophila simulans. Genetica 107, 113–120 (1999).CAS 
PubMed 
Article 

Google Scholar 
22.Buchman, A. B., Ivy, T., Marshall, J. M., Akbari, O. S. & Hay, B. A. Engineered reciprocal chromosome translocations drive high threshold, reversible population replacement in drosophila. ACS Synth. Biol. 7, 1359–1370 (2018).CAS 
PubMed 
Article 

Google Scholar 
23.Akbari, O. S. et al. Novel synthetic Medea selfish genetic elements drive population replacement in Drosophila; a theoretical exploration of Medea-dependent population suppression. ACS Synth. Biol. 3, 915–928 (2014).CAS 
PubMed 
Article 

Google Scholar 
24.Buchman, A., Marshall, J. M., Ostrovski, D., Yang, T. & Akbari, O. S. Synthetically engineered Medea gene drive system in the worldwide crop pest Drosophila suzukii. Proc. Natl Acad. Sci. USA 115, 4725–4730 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
25.Champer, J., Zhao, J., Champer, S. E., Liu, J. & Messer, P. W. Population dynamics of underdominance gene drive systems in continuous space. ACS Synth. Biol. 9, 779–792 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
26.Chen, C. C. et al. EXO1 suppresses double-strand break induced homologous recombination between diverged sequences in mammalian cells. DNA Repair. 57, 98–106 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
27.Leftwich, P. T. et al. Recent advances in threshold-dependent gene drives for mosquitoes. Biochem. Soc. Trans. 46, 1203–1212 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
28.Raban, R. R., Marshall, J. M. & Akbari, O. S. Progress towards engineering gene drives for population control. J. Exp. Biol. 223 (Suppl. 1), jeb208181 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
29.Ward, C. M. et al. Medea selfish genetic elements as tools for altering traits of wild populations: a theoretical analysis. Evolution 65, 1149–1162 (2011).PubMed 
Article 

Google Scholar 
30.Oberhofer, G., Ivy, T. & Hay, B. A. Gene drive and resilience through renewal with next generation Cleave and Rescue selfish genetic elements. Proc. Natl Acad. Sci. USA 117, 9013–9021 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
31.Oberhofer, G., Ivy, T. & Hay, B. A. Cleave and Rescue, a novel selfish genetic element and general strategy for gene drive. Proc. Natl Acad. Sci. USA 116, 6250–6259 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
32.Champer, J. et al. A toxin-antidote CRISPR gene drive system for regional population modification. Nat. Commun. 11, 1082 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
33.Yen, P. S. & Failloux, A. B. A review: Wolbachia-based population replacement for mosquito control shares common points with genetically modified control approaches. Pathogens 9, 404 (2020).PubMed Central 
Article 
PubMed 

Google Scholar 
34.O’Neill, S. L. The use of wolbachia by the world mosquito program to interrupt transmission of aedes aegypti transmitted viruses. Adv. Exp. Med. Biol. 1062, 355–360 (2018).PubMed 
Article 
CAS 

Google Scholar 
35.Niang, E. H. A., Bassene, H., Fenollar, F. & Mediannikov, O. Biological control of mosquito-borne diseases: the potential of wolbachia-based interventions in an IVM framework. J. Trop. Med. 2018, 1470459 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
36.Chevalier, B. S. & Stoddard, B. L. Homing endonucleases: structural and functional insight into the catalysts of intron/intein mobility. Nucleic Acids Res. 29, 3757–3774 (2001).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
37.Macreadie, I. G., Scott, R. M., Zinn, A. R. & Butow, R. A. Transposition of an intron in yeast mitochondria requires a protein encoded by that intron. Cell 41, 395–402 (1985).CAS 
PubMed 
Article 

Google Scholar 
38.Rong, Y. S. & Golic, K. G. The homologous chromosome is an effective template for the repair of mitotic DNA double-strand breaks in Drosophila. Genetics 165, 1831–1842 (2003).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
39.Chan, Y. S., Huen, D. S., Glauert, R., Whiteway, E. & Russell, S. Optimising homing endonuclease gene drive performance in a semi-refractory species: the Drosophila melanogaster experience. PLoS ONE 8, e54130 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
40.Windbichler, N. et al. A synthetic homing endonuclease-based gene drive system in the human malaria mosquito. Nature 473, 212–215 (2011). This study is the first demonstration of nuclease-mediated gene drive in mosquitoes based on the homing endonuclease gene I-SceI.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
41.Carroll, D. Genome engineering with targetable nucleases. Annu. Rev. Biochem. 83, 409–439 (2014).CAS 
PubMed 
Article 

Google Scholar 
42.Barrangou, R. et al. CRISPR provides acquired resistance against viruses in prokaryotes. Science 315, 1709–1712 (2007).CAS 
Article 
PubMed 

Google Scholar 
43.Jinek, M. et al. A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity. Science 337, 816–821 (2012). This foundational study developed the most widely used dual synthetic CRISPR system consisting of Cas9 endonuclease and gRNA components.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
44.Doudna, J. A., Sternberg, S. H. A Crack in Creation: Gene Editing and the Unthinkable Power to Control Evolution 281 (Houghton Mifflin Harcourt, 2017).45.Gantz, V. M. & Bier, E. The mutagenic chain reaction: a method for converting heterozygous to homozygous mutations. Science 348, 442–444 (2015). This study reported the first CRISPR-based gene drive in a metazoan organism (D. melanogaster) with a specialized germline.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
46.Gantz, V. M. et al. Highly efficient Cas9-mediated gene drive for population modification of the malaria vector mosquito Anopheles stephensi. Proc. Natl Acad. Sci. USA 112, E6736–E6743 (2015). This study describes the first efficient CRISPR-based gene drive system in mosquitoes, which carried a dual anti-malarial effector cassette.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
47.Hammond, A. et al. A CRISPR-Cas9 gene drive system targeting female reproduction in the malaria mosquito vector Anopheles gambiae. Nat. Biotechnol. 34, 78–83 (2016). This study describes the first efficient CRISPR-based suppression gene drive system in mosquitoes.CAS 
PubMed 
Article 

Google Scholar 
48.Kyrou, K. et al. A CRISPR-Cas9 gene drive targeting doublesex causes complete population suppression in caged Anopheles gambiae mosquitoes. Nat. Biotechnol. 36, 1062–1066 (2018). This study describes a highly efficient suppression gene drive system in mosquitoes targeting an invariant genome target site in the doublesex locus.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
49.Li, M. et al. Development of a confinable gene drive system in the human disease vector Aedes aegypti. eLife 9, e51701 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
50.Grunwald, H. A. et al. Super-Mendelian inheritance mediated by CRISPR-Cas9 in the female mouse germline. Nature 566, 105–109 (2019). This study provided the first proof-of-principle gene drive system in mammals, which selectively sustained drive via the female germline.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
51.DiCarlo, J. E., Chavez, A., Dietz, S. L., Esvelt, K. M. & Church, G. M. Safeguarding CRISPR-Cas9 gene drives in yeast. Nat. Biotechnol. 33, 1250–1255 (2015). This study demonstrated CRISPR-based gene conversion in diploid yeast, which could then be transmitted meiotically.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
52.Valderrama, J. A., Kulkarni, S. S., Nizet, V. & Bier, E. A bacterial gene-drive system efficiently edits and inactivates a high copy number antibiotic resistance locus. Nat. Commun. 10, 5726 (2019). This study generalizes the concept of gene drive to bacteria, where it is applied to efficiently reduce the frequency of antibiotic reistance.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
53.Esvelt, K. M., Smidler, A. L., Catteruccia, F. & Church, G. M. Concerning RNA-guided gene drives for the alteration of wild populations. eLife 3, e03401 (2014).PubMed 
PubMed Central 
Article 

Google Scholar 
54.Adolfi, A. et al. Efficient population modification gene-drive rescue system in the malaria mosquito Anopheles stephensi. Nat. Commun. 11, 5553 (2020). This study reports on the first recoded gene drive in mosquitoes that drove efficiently through both males and females based on the process of lethal/sterile mosaicism.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
55.Champer, J. et al. A CRISPR homing gene drive targeting a haplolethal gene removes resistance alleles and successfully spreads through a cage population. Proc. Natl Acad. Sci. USA 117, 24377–24383 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
56.Kandul, N. P., Liu, J., Bennett, J. B., Marshall, J. M. & Akbari, O. S. A confinable home-and-rescue gene drive for population modification. eLife 10, e65939 (2021).PubMed 
PubMed Central 
Article 

Google Scholar 
57.Terradas, G. et al. Inherently confinable split-drive systems in Drosophila. Nat. Commun. 12, 1480 (2021). This study further develops the strategy of inserting a recoded gene drive in genes essential for viability or reproduction in the context of split drive systems.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
58.Xu, X. S., Gantz, V. M., Siomava, N. & Bier, E. CRISPR/Cas9 and active genetics-based trans-species replacement of the endogenous Drosophila kni-L2 CRM reveals unexpected complexity. eLife 6, e30281 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
59.Lopez Del Amo, V. et al. A transcomplementing gene drive provides a flexible platform for laboratory investigation and potential field deployment. Nat. Commun. 11, 352 (2020). This study reports on the reconstitution of a full gene drive from split constituent parts.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
60.Guichard, A. et al. Efficient allelic-drive in Drosophila. Nat. Commun. 10, 1640 (2019). The study develops two allelic drive systems, copy-cutting and copy-grafting, to propagate favoured alleles of an essential gene.PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
61.Kandul, N. P. et al. Assessment of a split homing based gene drive for efficient knockout of multiple genes. G3 10, 827–837 (2020).CAS 
PubMed 
Article 

Google Scholar 
62.Xu, X.-R. S. et al. Active-genetic neutralizing elements for halting or deleting gene-drives. Mol. Cell 80, 246–262 (2020). This study reports on two drive-neutralizing systems that either inactivate (e-CHACR) or delete and replace (ERACR) a gene drive.CAS 
PubMed 
Article 

Google Scholar 
63.Burt, A. Site-specific selfish genes as tools for the control and genetic engineering of natural populations. Proc. Biol. Sci. 270, 921–928 (2003). This seminal modelling study provides the theoretical underpinnings for the modern gene-drive field.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
64.North, A. R., Burt, A. & Godfray, H. C. J. Modelling the potential of genetic control of malaria mosquitoes at national scale. BMC Biol. 17, 26 (2019). This study provides a comprehensive analysis of the perfomance of suppressive gene drives following iterative releases across various topographies.PubMed 
PubMed Central 
Article 

Google Scholar 
65.North, A. R., Burt, A. & Godfray, H. C. J. Modelling the suppression of a malaria vector using a CRISPR-Cas9 gene drive to reduce female fertility. BMC Biol. 18, 98 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
66.Collins, C. M., Bonds, J. A. S., Quinlan, M. M. & Mumford, J. D. Effects of the removal or reduction in density of the malaria mosquito, Anopheles gambiae s.l., on interacting predators and competitors in local ecosystems. Med. Vet. Entomol. 33, 1–15 (2019).CAS 
PubMed 
Article 

Google Scholar 
67.James, A. A. Gene drive systems in mosquitoes: rules of the road. Trends Parasitol. 21, 64–67 (2005).CAS 
PubMed 
Article 

Google Scholar 
68.Gantz, V. M. & Bier, E. The dawn of active genetics. Bioessays 38, 50–63 (2016).PubMed 
Article 

Google Scholar 
69.Macias, V. M. & James, A. A. in Genetic Control of Malaria and Dengue (ed. Adelman, Z. N.) 423–444 (Elsevier Academic Press, 2015).70.Eckhoff, P. A., Wenger, E. A., Godfray, H. C. & Burt, A. Impact of mosquito gene drive on malaria elimination in a computational model with explicit spatial and temporal dynamics. Proc. Natl Acad. Sci. USA 114, E255–E264 (2017). This study provides a detailed analysis of drive parameters relevant to both suppression-based and modification-based drives and is the first to model a drive in the context of a two-dimensional environment.CAS 
PubMed 
Article 

Google Scholar 
71.Hammond, A. M. et al. The creation and selection of mutations resistant to a gene drive over multiple generations in the malaria mosquito. PLoS Genet. 13, e1007039 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
72.Joyce, E. F., Paul, A., Chen, K. E., Tanneti, N. & McKim, K. S. Multiple barriers to nonhomologous DNA end joining during meiosis in Drosophila. Genetics 191, 739–746 (2012).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
73.Bozas, A., Beumer, K. J., Trautman, J. K. & Carroll, D. Genetic analysis of zinc-finger nuclease-induced gene targeting in Drosophila. Genetics 182, 641–651 (2009).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
74.Do, A. T., Brooks, J. T., Le Neveu, M. K. & LaRocque, J. R. Double-strand break repair assays determine pathway choice and structure of gene conversion events in Drosophila melanogaster. G3 4, 425–432 (2014).PubMed 
Article 
CAS 

Google Scholar 
75.Wei, D. S. & Rong, Y. S. A genetic screen for DNA double-strand break repair mutations in Drosophila. Genetics 177, 63–77 (2007).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
76.Lin, C. C. & Potter, C. J. Non-Mendelian dominant maternal effects caused by CRISPR/Cas9 transgenic components in Drosophila melanogaster. G3 6, 3685–3691 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
77.Champer, J. et al. Novel CRISPR/Cas9 gene drive constructs reveal insights into mechanisms of resistance allele formation and drive efficiency in genetically diverse populations. PLoS Genet. 13, e1006796 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
78.Anopheles gambiae 1000 Genomes Consortiumet al. Genetic diversity of the African malaria vector Anopheles gambiae. Nature 552, 96–100 (2017).Article 
CAS 

Google Scholar 
79.Deredec, A., Burt, A. & Godfray, H. C. The population genetics of using homing endonuclease genes in vector and pest management. Genetics 179, 2013–2026 (2008).PubMed 
PubMed Central 
Article 

Google Scholar 
80.Fasulo, B. et al. A fly model establishes distinct mechanisms for synthetic CRISPR/Cas9 sex distorters. PLoS Genet. 16, e1008647 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
81.Galizi, R. et al. A synthetic sex ratio distortion system for the control of the human malaria mosquito. Nat. Commun. 5, 3977 (2014).CAS 
PubMed 
Article 

Google Scholar 
82.Galizi, R. et al. A CRISPR-Cas9 sex-ratio distortion system for genetic control. Sci. Rep. 6, 31139 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
83.Turner, J. M. Meiotic sex chromosome inactivation. Development 134, 1823–1831 (2007).CAS 
PubMed 
Article 

Google Scholar 
84.Simoni, A. et al. A male-biased sex-distorter gene drive for the human malaria vector Anopheles gambiae. Nat. Biotechnol. 38, 1054–1060 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
85.Carballar-Lejarazu, R. & et al. Next-generation gene drive for population modification of the malaria vector mosquito, Anopheles gambiae.Proc. Natl Acad. Sci. USA 117, 22805–22814 (2020). This study describes a modification gene drive that propagates with high efficiency through both males and females.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
86.Pham, T. B. et al. Experimental population modification of the malaria vector mosquito, Anopheles stephensi. PLoS Genet. 15, e1008440 (2019).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
87.Dong, Y., Simoes, M. L. & Dimopoulos, G. Versatile transgenic multistage effector-gene combinations for Plasmodium falciparum suppression in Anopheles. Sci. Adv. 6, eaay5898 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
88.Dong, Y. et al. Engineered anopheles immunity to Plasmodium infection. PLoS Pathog. 7, e1002458 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
89.Isaacs, A. T. et al. Transgenic Anopheles stephensi coexpressing single-chain antibodies resist Plasmodium falciparum development. Proc. Natl Acad. Sci. USA 109, E1922–E1930 (2012). This study demonstrates 100% protection against parasite transmission in transgenic mosquitoes carrying a dual anti-parasite effector cassette.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
90.Haber, J. E. TOPping off meiosis. Mol. Cell 57, 577–581 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
91.Hammond, A. et al. Regulating the expression of gene drives is key to increasing their invasive potential and the mitigation of resistance. PLoS Genet. 17, e1009321 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
92.Lee, Y. et al. Genome-wide divergence among invasive populations of Aedes aegypti in California. BMC Genomics 20, 204 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
93.Callaway, E. Gene drives thwarted by emergence of resistant organisms. Nature 542, 15 (2017).CAS 
PubMed 
Article 

Google Scholar 
94.Unckless, R. L., Clark, A. G. & Messer, P. W. Evolution of resistance against CRISPR/Cas9 gene drive. Genetics 205, 827–841 (2017).PubMed 
Article 

Google Scholar 
95.Drury, D. W., Dapper, A. L., Siniard, D. J., Zentner, G. E. & Wade, M. J. CRISPR/Cas9 gene drives in genetically variable and nonrandomly mating wild populations. Sci. Adv. 3, e1601910 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
96.Schmidt, H. et al. Abundance of conserved CRISPR-Cas9 target sites within the highly polymorphic genomes of Anopheles and Aedes mosquitoes. Nat. Commun. 11, 1425 (2020). This study provides computational evidence that conserved CRISPR cleavage sites are abundant in the genome.CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
97.Akbari, O. S. et al. Safeguarding gene drive experiments in the laboratory. Science 349, 927–929 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
98.Li, J. et al. Genome-block expression-assisted association studies discover malaria resistance genes in Anopheles gambiae. Proc. Natl Acad. Sci. USA 110, 20675–20680 (2013).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
99.Niu, G. et al. The fibrinogen-like domain of FREP1 protein is a broad-spectrum malaria transmission-blocking vaccine antigen. J. Biol. Chem. 292, 11960–11969 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
100.Zhang, G. et al. Anopheles midgut FREP1 mediates plasmodium invasion. J. Biol. Chem. 290, 16490–16501 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
101.Dong, Y., Simoes, M. L., Marois, E. & Dimopoulos, G. CRISPR/Cas9 -mediated gene knockout of Anopheles gambiae FREP1 suppresses malaria parasite infection. PLoS Pathog. 14, e1006898 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
102.Simoes, M. L., Caragata, E. P. & Dimopoulos, G. Diverse host and restriction factors regulate mosquito-pathogen interactions. Trends Parasitol. 34, 603–616 (2018).PubMed 
Article 

Google Scholar 
103.Nash, A. et al. Integral gene drives for population replacement. Biol. Open 8, bio037762 (2019). This study describes a bipartite drive system that can enable testing of anti-parasite effector cassettes under standard mosquito confinement protocols.CAS 
PubMed 

Google Scholar 
104.Enayati, A., Hanafi-Bojd, A. A., Sedaghat, M. M., Zaim, M. & Hemingway, J. Evolution of insecticide resistance and its mechanisms in Anopheles stephensi in the WHO Eastern Mediterranean Region. Malar. J. 19, 258 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
105.Ffrench-Constant, R. H., Williamson, M. S., Davies, T. G. & Bass, C. Ion channels as insecticide targets. J. Neurogenet. 30, 163–177 (2016).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
106.Silva, J. J. & Scott, J. G. Conservation of the voltage-sensitive sodium channel protein within the Insecta. Insect Mol. Biol. 29, 9–18 (2020).CAS 
PubMed 
Article 

Google Scholar 
107.Casida, J. E. & Durkin, K. A. Novel GABA receptor pesticide targets. Pestic. Biochem. Physiol. 121, 22–30 (2015).CAS 
PubMed 
Article 

Google Scholar 
108.Ihara, M., Buckingham, S. D., Matsuda, K. & Sattelle, D. B. Modes of action, resistance and toxicity of insecticides targeting nicotinic acetylcholine receptors. Curr. Med. Chem. 24, 2925–2934 (2017).CAS 
PubMed 
Article 

Google Scholar 
109.Thapa, S., Lv, M. & Xu, H. Acetylcholinesterase: a primary target for drugs and insecticides. Mini Rev. Med. Chem. 17, 1665–1676 (2017).CAS 
PubMed 
Article 

Google Scholar 
110.Kleinstiver, B. P. et al. Engineered CRISPR-Cas12a variants with increased activities and improved targeting ranges for gene, epigenetic and base editing. Nat. Biotechnol. 37, 276–282 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
111.Walton, R. T., Christie, K. A., Whittaker, M. N. & Kleinstiver, B. P. Unconstrained genome targeting with near-PAMless engineered CRISPR-Cas9 variants. Science 368, 290–296 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
112.Committee on Gene Drive Research in Non-Human Organisms: Recommendations for Responsible Conduct; Board on Life Sciences; Division on Earth and Life Studies; National Academies of Sciences, Engineering, and Medicine. Gene Drives on the Horizon: Advancing Science, Navigating Uncertainty, and Aligning Research with Public Values (The National Academies Press, 2016). This comprehensive advisory and historical review document summarizes consensus views for how to safely rear and study gene-drive systems in the laboratory.113.Adelman, Z. et al. Rules of the road for insect gene drive research and testing. Nat. Biotechnol. 35, 716–718 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
114.James, S. et al. Pathway to deployment of gene drive mosquitoes as a potential biocontrol tool for elimination of malaria in Sub-Saharan Africa: recommendations of a scientific working group(dagger). Am. J. Trop. Med. Hyg. 98, 1–49 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
115.James, S. L., Marshall, J. M., Christophides, G. K., Okumu, F. O. & Nolan, T. Toward the definition of efficacy and safety criteria for advancing gene drive-modified mosquitoes to field testing. Vector Borne Zoonotic Dis. 20, 237–251 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
116.Warmbrod, K. L. et al. Gene Drives: Pursuing Opportunities, Minimizing Risk – A Johns Hopkins University Report on Responsible Governance (Johns Hopkins Bloomberg School of Public Health, Center for Health Security, Johns Hopkins University, 2020).117.Vella, M. R., Gunning, C. E., Lloyd, A. L. & Gould, F. Evaluating strategies for reversing CRISPR-Cas9 gene drives. Sci. Rep. 7, 11038 (2017).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
118.Rode, N. O., Courtier-Orgogozo, V. & Debarre, F. Can a population targeted by a CRISPR-based homing gene drive be rescued? G3 10, 3403–3415 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
119.Fedoroff, N., Wessler, S. & Shure, M. Isolation of the transposable maize controlling elements Ac and Ds. Cell 35, 235–242 (1983).CAS 
PubMed 
Article 

Google Scholar 
120.Paix, A. et al. Precision genome editing using synthesis-dependent repair of Cas9-induced DNA breaks. Proc. Natl Acad. Sci. USA 114, E10745–E10754 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
121.Wu, B., Luo, L. & Gao, X. J. Cas9-triggered chain ablation of cas9 as a gene drive brake. Nat. Biotechnol. 34, 137–138 (2016).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
122.Taxiarchi, C. et al. A genetically encoded anti-CRISPR protein constrains gene drive spread and prevents population suppression. Nat. Commun. 12, 3977 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
123.Conklin, B. R. On the road to a gene drive in mammals. Nature 566, 43–45 (2019).CAS 
PubMed 
Article 

Google Scholar 
124.Salkeld, D. J. Vaccines for conservation: plague, prairie dogs & black-footed ferrets as a case study. Ecohealth 14, 432–437 (2017).PubMed 
Article 

Google Scholar 
125.Teem, J. L. et al. Genetic biocontrol for invasive species. Front. Bioeng. Biotechnol. 8, 452 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
126.Godwin, J. et al. Rodent gene drives for conservation: opportunities and data needs. Proc. Biol. Sci. 286, 20191606 (2019).PubMed 
PubMed Central 

Google Scholar 
127.McFarlane, G. R., Whitelaw, C. B. A. & Lillico, S. G. CRISPR-based gene drives for pest control. Trends Biotechnol. 36, 130–133 (2018).CAS 
PubMed 
Article 

Google Scholar 
128.Klompe, S. E., Vo, P. L. H., Halpin-Healy, T. S. & Sternberg, S. H. Transposon-encoded CRISPR-Cas systems direct RNA-guided DNA integration. Nature 571, 219–225 (2019).CAS 
PubMed 
Article 

Google Scholar 
129.Koonin, E. V., Makarova, K. S., Wolf, Y. I. & Krupovic, M. Evolutionary entanglement of mobile genetic elements and host defence systems: guns for hire. Nat. Rev. Genet. 21, 119–131 (2020).CAS 
PubMed 
Article 

Google Scholar 
130.Strecker, J. et al. RNA-guided DNA insertion with CRISPR-associated transposases. Science 365, 48–53 (2019).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
131.Peters, J. E., Makarova, K. S., Shmakov, S. & Koonin, E. V. Recruitment of CRISPR-Cas systems by Tn7-like transposons. Proc. Natl Acad. Sci. USA 114, E7358–E7366 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
132.Wiegand, T. & Wiedenheft, B. CRISPR Surveillance Turns Transposon Taxi. CRISPR J. 3, 10–12 (2020).PubMed 
PubMed Central 
Article 

Google Scholar 
133.Hamilton, T. A. et al. Efficient inter-species conjugative transfer of a CRISPR nuclease for targeted bacterial killing. Nat. Commun. 10, 4544 (2019).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
134.Price, V. J. et al. Enterococcus faecalis CRISPR-cas is a robust barrier to conjugative antibiotic resistance dissemination in the murine intestine. mSphere 4, e00464-19 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
135.Rodrigues, M., McBride, S. W., Hullahalli, K., Palmer, K. L. & Duerkop, B. A. Conjugative delivery of CRISPR-Cas9 for the selective depletion of antibiotic-resistant enterococci. Antimicrob. Agents Chemother. 63, e01454-19 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
136.Carraro, N. et al. Plasmid-like replication of a minimal streptococcal integrative and conjugative element. Microbiology 162, 622–632 (2016).CAS 
PubMed 
Article 

Google Scholar 
137.Brophy, J. A. N. et al. Engineered integrative and conjugative elements for efficient and inducible DNA transfer to undomesticated bacteria. Nat. Microbiol. 3, 1043–1053 (2018).CAS 
PubMed 
Article 

Google Scholar 
138.Bikard, D. et al. Exploiting CRISPR-Cas nucleases to produce sequence-specific antimicrobials. Nat. Biotechnol. 32, 1146–1150 (2014).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
139.Citorik, R. J., Mimee, M. & Lu, T. K. Sequence-specific antimicrobials using efficiently delivered RNA-guided nucleases. Nat. Biotechnol. 32, 1141–1145 (2014).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
140.Yosef, I., Manor, M., Kiro, R. & Qimron, U. Temperate and lytic bacteriophages programmed to sensitize and kill antibiotic-resistant bacteria. Proc. Natl Acad. Sci. USA 112, 7267–7272 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
141.Park, J. Y. et al. Genetic engineering of a temperate phage-based delivery system for CRISPR/Cas9 antimicrobials against Staphylococcus aureus. Sci. Rep. 7, 44929 (2017).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
142.Pazda, M., Kumirska, J., Stepnowski, P. & Mulkiewicz, E. Antibiotic resistance genes identified in wastewater treatment plant systems – a review. Sci. Total. Env. 697, 134023 (2019).CAS 
Article 

Google Scholar 
143.Kraemer, S. A., Ramachandran, A. & Perron, G. G. Antibiotic pollution in the environment: from microbial ecology to public policy. Microorganisms 7, 180 (2019).CAS 
PubMed Central 
Article 
PubMed 

Google Scholar 
144.Ram, G., Ross, H. F., Novick, R. P., Rodriguez-Pagan, I. & Jiang, D. Conversion of staphylococcal pathogenicity islands to CRISPR-carrying antibacterial agents that cure infections in mice. Nat. Biotechnol. 36, 971–976 (2018).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
145.Bier, E. & Nizet, V. Driving to safety: CRISPR-based genetic approaches to reducing antibiotic resistance. Trends Genet. https://doi.org/10.1016/j.tig.2021.02.007 (2021).Article 
PubMed 

Google Scholar 
146.Rossati, A. et al. Climate, environment and transmission of malaria. Infez. Med. 24, 93–104 (2016).PubMed 

Google Scholar 
147.Fontenille, D. & Powell, J. R. From anonymous to public enemy: how does a mosquito become a feared arbovirus vector? Pathogens 9, 265 (2020).PubMed Central 
Article 
PubMed 

Google Scholar 
148.Lidani, K. C. F. et al. Chagas disease: from discovery to a worldwide health problem. Front. Public Health 7, 166 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
149.Buscher, P., Cecchi, G., Jamonneau, V. & Priotto, G. Human African trypanosomiasis. Lancet 390, 2397–2409 (2017).PubMed 
Article 

Google Scholar 
150.Desjeux, P. Leishmaniasis: current situation and new perspectives. Comp. Immunol. Microbiol. Infect. Dis. 27, 305–318 (2004).CAS 
PubMed 
Article 

Google Scholar 
151.Saxena, V., Bolling, B. G. & Wang, T. West nile virus. Clin. Lab. Med. 37, 243–252 (2017).PubMed 
Article 

Google Scholar 
152.Simon, L. V., Kong, E. L. & Graham, C. in St. Louis Encephalitis (StatPearls, 2020).153.Feng, X. et al. Optimized CRISPR tools and site-directed transgenesis towards gene drive development in Culex quinquefasciatus mosquitoes. Nat. Commun. 12, 2960 (2021).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
154.Nepomichene, T. N., Andrianaivolambo, L., Boyer, S. & Bourgouin, C. Efficient method for establishing F1 progeny from wild populations of Anopheles mosquitoes. Malar. J. 16, 21 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
155.Marchand, R. P. A new cage for observing mating behavior of wild Anopheles gambiae in the laboratory. J. Am. Mosq. Control. Assoc. 1, 234–236 (1985).CAS 
PubMed 

Google Scholar 
156.Nunes-da-Fonseca, R., Berni, M., Tobias-Santos, V., Pane, A. & Araujo, H. M. Rhodnius prolixus: from classical physiology to modern developmental biology. Genesis https://doi.org/10.1002/dvg.22995 (2017).Article 
PubMed 

Google Scholar 
157.Chaverra-Rodriguez, D. et al. Targeted delivery of CRISPR-Cas9 ribonucleoprotein into arthropod ovaries for heritable germline gene editing. Nat. Commun. 9, 3008 (2018).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
158.Macias, V. M. et al. Cas9-mediated gene-editing in the malaria mosquito anopheles stephensi by ReMOT Control. G3 10, 1353–1360 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
159.Chaverra-Rodriguez, D. et al. Germline mutagenesis of Nasonia vitripennis through ovarian delivery of CRISPR-Cas9 ribonucleoprotein. Insect Mol. Biol. 29, 569–577 (2020).CAS 
PubMed 
Article 

Google Scholar 
160.Heu, C. C., McCullough, F. M., Luan, J. & Rasgon, J. L. CRISPR-Cas9-based genome editing in the silverleaf whitefly (Bemisia tabaci). CRISPR J. 3, 89–96 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
161.Prowse, T. A., Adikusuma, F., Cassey, P., Thomas, P. & Ross, J. V. A Y-chromosome shredding gene drive for controlling pest vertebrate populations. eLife 8, e41873 (2019).PubMed 
PubMed Central 
Article 

Google Scholar 
162.Carballar-Lejarazu, R. & James, A. A. Population modification of Anopheline species to control malaria transmission. Pathog. Glob. Health 111, 424–435 (2017).PubMed 
Article 
PubMed Central 

Google Scholar 
163.Annas, G. J. et al. A code of ethics for gene drive research. CRISPR J. 4, 19–24 (2021).PubMed 
PubMed Central 
Article 

Google Scholar 
164.Bier, E. & Sober, E. Gene editing and the war against malaria. Am. Sci. 108, 162–169 (2020).Article 

Google Scholar 
165.Long, K. C. et al. Core commitments for field trials of gene drive organisms. Science 370, 1417–1419 (2020).CAS 
PubMed 
Article 

Google Scholar 
166.Kormos, A. et al. Application of the relationship-based model to engagement for field trials of genetically engineered malaria vectors.Am. J. Trop. Med. Hyg. 104, 805–811 (2020).PubMed Central 
PubMed 

Google Scholar 
167.World Health Organization. Guidance framework for testing of genetically modified mosquitoes. WHO http://apps.who.int/iris/bitstream/10665/127889/1/9789241507486_eng.pdf (2014).168.Smith, D. L., McKenzie, F. E., Snow, R. W. & Hay, S. I. Revisiting the basic reproductive number for malaria and its implications for malaria control. PLoS Biol. 5, e42 (2007).PubMed 
PubMed Central 
Article 
CAS 

Google Scholar 
169.Brauer, F., Castillo-Chavez, C., Mubayi, A. & Towers, S. Some models for epidemics of vector-transmitted diseases. Infect. Dis. Model. 1, 79–87 (2016).PubMed 
PubMed Central 

Google Scholar 
170.Deredec, A., Godfray, H. C. & Burt, A. Requirements for effective malaria control with homing endonuclease genes. Proc. Natl Acad. Sci. USA 108, E874–E880 (2011).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
171.Escalante, A. A. & Pacheco, M. A. Malaria molecular epidemiology: an evolutionary genetics perspective. Microbiol. Spectr. https://doi.org/10.1128/microbiolspec.AME-0010-2019 (2019).Article 
PubMed 
PubMed Central 

Google Scholar 
172.Selvaraj, P. et al. Vector genetics, insecticide resistance and gene drives: An agent-based modeling approach to evaluate malaria transmission and elimination. PLoS Comput. Biol. 16, e1008121 (2020).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar  More

JSATS acoustic receiver deploymentsTwo types of acoustic receivers were deployed for this study: cabled and autonomous acoustic receivers (Table 1). These systems were deployed in the Snake River in eastern Washington State and in the lower Columbia River along the Oregon/Washington border. Two hydropower dams were outfitted with cabled acoustic receiver arrays: LGS and LMN. This allowed passage routes to be defined at each of these dams for each of the tagged fish. Transects of autonomous acoustic receivers (i.e., receiver arrays) were deployed immediately upstream and downstream of each of these dams for computing dam passage survival, and at several other locations to estimate cumulative mortality. From the release site to the final autonomous acoustic receiver array, fish detected at the final array will have passed through LGS, LMN, Ice Harbor, McNary, John Day, The Dalles, and Bonneville Dams. Each of these dams feature powerhouses with large, vertical Kaplan runners19 and spillways utilizing radial or vertical lift spillway gates. Most also include juvenile bypass systems (JBS; Table 1) and surface spill weirs20 to provide safer passage routes for fish.Before deployment, all hydrophones and receivers were evaluated in an acoustic tank lined with anechoic materials at the PNNL Bio-Acoustics & Flow Laboratory (BFL21). The BFL is accredited by the American Association for Laboratory Accreditation to ISO/IEC 17025:2005, which is the international standard for calibration and testing laboratories. The accreditation scope (Certificate Number 3267.01) includes hydrophone sensitivity measurements and power-level measurements of sound sources for frequencies from 50 to 500 kHz for both military equipment and commercial components. The evaluation involved simulating transmissions from tags located at increasing distances. This allowed the performance of each receiver to be validated prior to deployment to ensure the expected detection range will likely be achieved.Cabled receivers – deployment, hardware, and data processingA JSATS cabled acoustic receiver system2 consists of up to four narrowband hydrophones; various types of hydrophone cables (e.g., four-channel “deck” cables, “y-blocks” that split the “deck” cables to individual connectors, and “wet” cables that run from the surface down into the water); a signal-conditioning, variable-gain amplifier; a data acquisition card that features a high-speed, analog-to-digital converter, a digital signal processor and field-programmable gate array; a GPS receiver for synchronizing time among multiple systems; a data acquisition computer; and software for detecting22 and decoding23,24 the acoustic waveforms. Deploying these systems within the forebay of a hydropower dam typically involves rigidly mounting slotted pipes to the upstream edge of the pier nose between the powerhouse and spillway bays. The cabled hydrophones are mounted on “trolleys” that have an L-shaped arm that protrudes and rides in the slot in the pipe and allows the hydrophone to be mounted pointing upstream. Conical baffles containing anechoic material are installed around the hydrophones to block noise coming from behind the hydrophone—either noise from the dam or reflections off the concrete. To obtain the locations of hydrophones that have been lowered below the water surface, survey equipment is used to sight the tip of the hydrophone as it is lowered down the pier nose. This provides the true direction and slope of the pier nose, which is used, along with the length of braided steel cable attached to the trolley to lower it into place, to calculate the 3-D location in space for each hydrophone. The individual “wet” cables for each hydrophone are routed to the forebay deck, where they are combined into “deck” cables carrying four signals using the “y-block” cables. The deck cables are routed to mobile trailers that house the acoustic receiver equipment. Acoustic beacons that send out JSATS tag-code signals every 15 to 60 s are deployed alongside several of the hydrophones across the array. These beacons are used primarily for quality control, to monitor (typically through the internet from an off-site location) the performance of each hydrophone to determine whether there is a reduction in performance so that any malfunctioning hydrophones can be repaired as soon as possible.Two main programs run on the JSATS cabled receiver data acquisition computer. The first program is an energy-based detector software22 that collects the raw acoustic waveforms whenever the hydrophone signal meets a prescribed set of criteria. The second is decoder software24 that processes the waveforms saved by the detector to determine whether there is a valid detection (i.e., a decoded signal that has a valid CRC). If there is a valid detection, the decoded tag-code is saved to a text file, along with the detection time and other metadata. The detector software writes the binary waveform files to the hard drive using the *.bwm file type. The decoder is configured to wait for *.bwm files to be generated. Once a new *.bwm file is detected, the decoder will open the file, decode the data contained in the file, and then change the file extension to *.com to indicate that this file was decoded. If the detector saves data faster than the decoder can process it, for example at a hydraulic structure that is generating large amounts of acoustic noise (e.g., spillways with vertical lift gates), the decoder is configured to skip waveform files to avoid falling behind. Although these files may be skipped, the data contained within these files will still be used because the two different file extensions allow for readily identifying which files were not decoded in real time so that they can be decoded offline in a separate processing step after retrieving the data. The data is physically collected every 1–2 weeks, by swapping out the data collection hard drives. To make sure that all detection waveforms are processed by the decoder, the collected hard drives are put into data processing machines to decode any files that still have the *.bwm file extension. After confirming that all files have been decoded, either in real time or through post-processing after data retrieval, the decoded data from every hydrophone is checked for gaps in data. If a hydrophone is functional, there should be no large gaps in decoded data, since there are multiple stationary acoustic beacons deployed with the cabled hydrophone receiver array.Data is filtered to remove potential false positive decodes. Data filtering for the JSATS cabled acoustic receivers (Fig. 2a) begins with a multipath filter, which removes decodes from multipath signal propagations (e.g., acoustic reflections off the surface/bottom). The multipath filter is used on the data from each individual hydrophone. Any decodes of the same tag-code that occur a very short time (e.g., typically 45 m), too far away temporally from adjacent points ( >10 min), and that result in unrealistic velocities (~2 m/s for the size of fish we studied). Further quality assurance filters remove points that occur before a transmitter was released, after the PIT tag in the fish is detected downstream, and after the acoustic transmitter is detected by downstream autonomous receiver arrays.Once the 3-D trajectories have been computed, they are used to assign passage routes through the dam for each tagged fish (Fig. 3a). The route assignment for each tagged fish is divided into three parts: main route, subroute, and hole (Table 2). The main route describes the part of the general dam structure through which the tagged fish passed. This includes the powerhouse, the spillway, and a generic category, “dam,” which is used for rare scenarios where there is confidence that the fish was physically present but a lack of confidence in specifically where the fish passed the dam. The subroute further divides the main passage route into different subcategories. For a main route of “spillway,” the two subroutes are the traditional (deep) spillbays and special surface weirs (surface spillbays20). The surface weirs reside within one of the spillbays and assignment to either of these two subroutes is made directly using the acoustic telemetry results. For a main route of “powerhouse,” the two subroutes are turbine and JBS. Assignment of the JBS subroute requires that the PIT tag of a fish assigned to the powerhouse was detected by the PIT tag readers within the JBS system; otherwise the fish is assigned the turbine subroute. PIT tag detections at dams where cabled hydrophone arrays have been deployed can be used to assign the JBS subroute, and PIT tag detections at the other dams along the migration route can serve as additional detection events. The hole assignment defines the specific powerhouse intake or spillbay where the passage occurred.Table 2 Dam passage routes (main route, subroute, and hole) at LGS and LMN Dams.Full size tablePassage routes are assigned using the last 3D tracked location and the last detection. Two methods are used because the ability to consistently track the transmitter can diminish as the tagged fish approaches or passes through the plane containing the hydrophones, and the last decoded transmission could be later than the last 3-D tracked location.When the last 3-D tracked location is used, a route is assigned according to whether the last tracked point is within a specific area. This area spans the entire dam plus 25 m on each side and extends from the dam face to 30 m upstream into the forebay. If the last tracked point is within this boundary, the 3-D track passage route is assigned to the bay corresponding to the Y coordinate in the local dam coordinate system. If the last tracked point is outside the piers on either side of the dam, the passage route is assigned to the nearest bay.Route assignment based on the last detection uses the last transmission that was detected by multiple hydrophones. The detections associated with this transmission are sorted by time, and the pier numbers for the two hydrophones on different piers that first detected this transmission are averaged; the passage bay corresponding to this average pier is assigned as the last detection passage route. The default final route assignment is the 3-D tracked route assignment. However, when the two methods indicate different main routes, subroutes, or a different hole that is more than two bays away, the 3-D tracks are manually reviewed, and a decision is made regarding which method should be relied on for the final route assignment.After the final route assignment, a final quality assurance step is to compare the final route assignment to the dam operations. In case a tag-code has been assigned to a closed passage route, the 3-D tracks are reviewed to consider the trajectory of the tagged fish and the location of the nearest open passage route. As previously mentioned, the ability to consistently track a transmitter is diminished as it approaches or passes through the plane containing the hydrophones. An example of when a tagged fish could initially be assigned to a closed passage route would be when a passage route with a strong attractive flow (e.g., surface spill weir) is adjacent to a closed passage bay.Autonomous receivers – deployment, hardware, and data processingA JSATS autonomous receiver (SR5000, Advanced Telemetry Systems [ATS], USA), along with the necessary deployment accessories, consists of a hydrophone that is connected to a cylindrical, positively buoyant, self-contained, battery-powered, autonomous acoustic receiver; a submerged buoy line; an acoustic release; a braided stainless-steel cable; and a steel anchor. These receivers are typically deployed according to the methods described by Titzler et al.10, which involves using a 34 kg or 57 kg (depending on flow) steel anchor to deploy the autonomous receiver system to the river bottom. The anchor is attached to the release side of an acoustic release using a braided stainless-steel cable. The fixed end of the acoustic release is attached to the autonomous receiver using a submerged buoy line. When the acoustic release is remotely triggered, it detaches from the anchor line, and the combined buoyancy of the acoustic receiver and the submerged buoy line bring the system up to the surface. To maintain the receiver orientation in the water column during deployment, a thin plastic sheet is folded around the cylindrical body of the receiver, creating an airfoil-like shape that keeps the receiver oriented in the flow direction. Attached to each autonomous receiver is a JSATS beacon that is similar to the JSATS beacons deployed with the cabled hydrophone receiver arrays and used in the same way. The length of time that the JSATS autonomous receivers can be deployed is largely dependent on the battery life, with data retrieval and battery changes typically done every 2–3 weeks.Although it is possible to use autonomous receivers to conduct 3-D tracking, the process is much more challenging than using the cabled hydrophone receivers because the receivers are not fixed at well-defined locations (e.g., changes in river currents could change the receiver’s depth and horizontal location relative to the anchor) and are not time synchronized with each other. The autonomous receivers are typically used to detect presence of tagged fish, although recent research has investigated methods to improve the ability to conduct 3-D tracking25. Deploying an autonomous receiver array entails deploying several receivers in a line spanning the width of a river with the detection ranges of the receivers overlapping slightly. This creates a virtual detection gate that can be used to determine when a tagged fish passes through this location in the river. In addition to simply determining the migration timing of tagged fish, the autonomous receiver arrays are typically used for analyzing dam passage survival and near-dam behavior (e.g., forebay residence time, tailrace egress time).The data collected by the autonomous receivers is filtered similarly to data from the cabled receivers (Fig. 2b). The primary difference is that the data from each individual autonomous receiver is processed entirely by itself. As a result, the single-hydrophone filter is not used, and a single-node PRI filter is used instead of the message PRI filter.After the event histories for both the cabled and the autonomous acoustic receivers have been determined, individual routes were cross-checked by tracing the chronology of detections of every tagged fish as it was detected along the river in the sequence of acoustic receiver arrays. Upstream movement past a dam or out-of-sequence detections were deemed anomalous detection events. These anomalous detection events could be a few receptions resulting from noise or repeated detections of a transmitter that had been dropped near a receiver array after fish or bird predation. If the apparent behavior was impossible for a live fish, the anomalous detection was excluded from the detection history used for subsequent analysis.JSATS transmittersThe injectable transmitters used in this study (Fig. 1b) were manufactured by PNNL. Each transmitter (model microV26, which is licensed to, and currently commercially available from, Advanced Telemetry Systems as Model SS400) was 15 mm long, had an outside diameter of 3.35 mm, a volume of 0.111 mL, and a mass of 0.216 g in air and 0.105 g in water. The transmitters are generally cylindrical; excess epoxy was eliminated to reduce the weight, and epoxy surrounding the transducer element was minimized. The transmitters had a nominal transmission rate of 1 pulse every 4.2 s. Nominal transmitter life was expected to be about 28 d at a 4.2 s pulse rate. The acoustic signal is transmitted using a carrier frequency of 416.7 kHz, a source level of approximately 156 dB (ref. to 1 µPa at 1 m), and a total signal duration of 477 µs. The transmitter emits a uniquely coded 31-bit signal2, resulting in more than 65,000 individual tag-codes, using binary phase-shift keyed (BPSK) signal encoding.Each fish also bore a PIT tag (HPT12, Biomark, USA; 12.5 mm x Ø2.03 mm; 0.106 g in air). PIT tag detections were used to assign fish to passage through the JBS at LGS and LMN to distinguish between fish that were assigned to a main route of powerhouse and the turbine or JBS subroutes.Tagged fishFor this study, 682 subyearling Chinook salmon (Oncorhynchus tshawytscha) were tagged with the injectable acoustic transmitter and released upstream of LGS Dam on the Snake River in Washington State, USA (Fig. 1a). The fish were obtained from the JBS at LMN Dam and selected using existing fish screening criteria utilized in previous juvenile salmon survival studies26. The fish selected for the study were held in holding tanks for 18 to 30 hours prior to tagging, and for 10 to 25 hours after tagging prior to release. The size criteria for tagged fish was also identical to other recent juvenile salmon survival studies26. For this study the fork-lengths ranged from 95 to 143 mm, and the weights ranged from 7.5 to 29.3 g (see Tagged Fish Data for information on each tagged fish).Tagging procedureWhile each anesthetized fish was at the data station for recording physical parameters, a second person inserted both a disinfected PIT tag and an injectable acoustic transmitter, assigned to a specific fish, into a sterilized 8-gauge stainless-steel hypodermic needle17. First, the injectable transmitter was placed into the needle, battery-end first. The PIT tag was then also inserted in the same needle. A sanitized plastic cap was then placed over each end of the needle to retain the tags. Once both tags had been placed in the needle, the tag loaded needle was handed to the surgeon working at the tagging station. Additional details for the tagging procedure are documented by Deng et al.11.Release procedureThe fish implanted with the injectable acoustic transmitters were released using the same methods as fish tagged with commercially available acoustic transmitters for a separate large-scale survival study26. All fish were tagged at LMN and transported in insulated totes by truck to the single release site (Fig. 1a). There were five release locations across the river at the release site, and equal numbers of fish were released at each of the five locations. Releases occurred for 11 consecutive days (between 22 June and 2 July, 2013) and were staggered between day and night.Data managementUse of JSATS can generate a large volume of data. An integrated suite of science-based tools known as the Hydropower Biological Evaluation Toolset (HBET; https://hydropassage.org/hbet)27 was developed to assist the characterization of hydraulic conditions at hydropower structures and to understand the potential impacts on aquatic life. HBET was initially developed to be utilized to facilitate use of the autonomous sensor technology known as Sensor Fish28. HBET allows researchers to use previously collected Sensor Fish data to design studies to evaluate hypotheses, archive field-collected data, process raw sensor data, compare different hydraulic structures or operating conditions, and to estimate the biological response for species with known dose-response relationships. More recently, HBET was adapted to also provide the functionality of archiving new or previously collected acoustic telemetry data and to produce visualizations from that data. Although it is not necessary to visualize the data set associated with this manuscript, PNNL offers free government and academic use of the HBET software package in the U.S. and a free 90-day trial version of the package to interested parties. More

SamplingSamples were collected on 28 August 2018 during a field campaign to Lake Cadagno29, Switzerland. In situ measurements and water collection was performed at the deepest part of the lake (21 m). Water was collected using a pump CTD system as described in Di Nezio et al.36. Online in situ data were obtained during a continuous downcast of the CTD-system from the water surface down to ~17.5 m depth. During the upcast, discrete water samples were collected from a total of 20 depths (between 12 m and 17 m) above, in, and below the chemocline for chemical analyses and from 3 depths for incubation experiments (13.7 m, 14 m, and 15.5 m).In Lake Cadagno, wind-driven internal waves lead to vertical shifts of the water masses and their corresponding physicochemical parameters56. While sampling, it was apparent that the depths of the individual water masses had slightly shifted between the down- and the upcast. Therefore, we corrected the water depths of the samples collected during the upcast so that the physicochemical parameters during sampling best matched those of the continuous downcast, to ensure that samples were assigned to the respective water mass that they originated from. A custom R script was employed for the depth correction. In brief, all parameters measured by the CTD-system during the upcast and the downcast were normalized to percent (with 100% as the maximum observed value, and 0% the minimum observed value). Per individual sampling depth (during the upcast, where the pump cast CTD remained stationary for some time), average values of conductivity, temperature, and pressure were calculated and converted to percent values. Then, the depth from the downcast profile was identified that best matched all calculated percent values. This was achieved by subtracting the percent values per parameter from all respective data points of the downcast profile. Absolute values of the calculated differences per data row were summed. The depth with the lowest resulting sum, i.e., with the most similar physicochemical parameters, was then chosen as the corrected depth.Chemical analyses, flux calculations, and rate determinationsFor chemical analyses, lake water from the individual sampling depths was sterile-filtered (0.2 µm, cellulose acetate filter) and frozen at −20 °C until analysis. Samples were analyzed with a QuAAtro39 autoanalyzer (Seal Analytical) using the methods described in Strickland and Parsons57 to determine concentrations of dissolved inorganic phosphorus (PO43−), nitrite (NO2−), nitrate (NO3−), and reactive silica (Si). Ammonium concentrations were determined from the same filtered samples using the colorimetric analysis described in Kempers et al.58. Molybdenum concentrations were determined from filtered samples after acidification with 1% HNO3 (69%, ROTIPURAN®, Roth) using an ICP-MS 7900 (Agilent, Santa Clara, USA). Molybdenum was analyzed on mass 95 in He-mode using a multi-element calibration SRM (21 elements, Bernd Kraft). The SRM NIST 1643f was analyzed in parallel to guarantee the quality of analyses. Concentrations of sulfide were determined colorimetrically from unfiltered Lake water samples, following Cline59.To calculate the turbulent flux (J) of ammonium into the chemocline, we assumed a steady-state using Fick’s first law: J = −D∂C/∂x. A turbulent diffusion coefficient (D) of 1.6 × 10−6 m2 s−1 was used, corresponding to turbulence at the Lake Cadagno chemocline boundaries60. The change in concentration (∂C) was calculated over 14.25 m to 14.77 m depth, where the steepest ammonium gradient was observed. Ammonium uptake rates were calculated for the chemocline by integrating this flux over the chemocline from 13.45 m to 14.45 m depth.To quantify N2 fixation and primary production (i.e., CO2 fixation) rates, stable isotope incubations with 15N2 and 13CO2 were performed using established protocols61. Briefly, lake water from three different depths of the chemocline was sampled directly from the CTD pump system into five 250 ml serum bottles per depth. Water was filled into the bottles from bottom to top, allowing 1–2 bottle volumes to overflow to minimize oxygen contamination before crimp-sealing the bottles headspace-free with butyl rubber stoppers. Back in the field laboratory, no more than 8 h after sampling, one bottle per depth was filtered onto pre-combusted (460 °C, 6 h) glass microfiber filters (GF/F, Whatman®, UK) for in situ natural abundance of C and N. 13C-labeled sodium bicarbonate (NaH13CO3, 98 atom% 13C, dissolved in autoclaved MilliQ water; Sigma-Aldrich) was injected (320 µL) into three bottles per depth, to achieve a final concentration of 160 µmol L−1. Then, a volume of 5 ml 15N2 gas (Cambridge Isotope Laboratories, >98 atom% 15N, Lot #: I-19197/AR0586172) was injected as a bubble into the same bottles and shaken for 20 min to equilibrate the 15N2 gas. Sulfide solution was injected aiming for a final concentration of approximately 2 µM to remove trace oxygen contamination in the incubation bottles. Finally, the 15N2 gas bubble was replaced by anoxic in situ lake water from the respective depth. The bottles, together with one untreated control bottle per depth (containing unamended lake water), were incubated for a full light-dark cycle (13 h light, 11 h dark) under natural light conditions (0–8267 Lux, average: 247 Lux, median: 10.8 Lux, as determined by a HOBO pendant data logger, Onset Computer Corporation, Bourne, USA) in a water bath kept at ~12 °C.After incubation, samples were filtered onto pre-combusted GF/F filters. The filters were dried at room temperature and frozen at −20 °C for transport and storage. In addition, subsamples for nanoscale secondary ion mass spectrometry (nanoSIMS) analysis and for the determination of 13C and 15N enrichments in the substrate pools were taken from all bottles amended with 13C and 15N. NanoSIMS samples were fixed with 2% (final w/v) formaldehyde solution for 1 h at room temperature, prior to filtration onto gold-sputtered 0.22 µm polycarbonate membrane filters (GTTP IsoporeTM, Merck Millipore, USA). Subsamples for label% determinations were taken in gas-tight glass vials (Exetainer Labco, UK) and biological activity was terminated with HgCl2.Samples on GF/F filters were analyzed for C and N content and the respective isotopic composition by an elemental analyzer (Thermo Flash EA, 1112 Series) coupled to a continuous-flow isotope ratio mass spectrometer (Delta Plus XP IRMS; Thermo Finnigan, Dreieich, Germany). Enrichment of 15N in the N2 pool was determined using a membrane inlet mass spectrometer (MIMS; GAM200, IPI). Enrichment of 13C in the dissolved inorganic carbon pool was determined from 13C/12C-CO2 ratios after sample acidification with phosphoric acid using cavity ring-down spectroscopy (G2201-I coupled to a Liaison A0301, Picarro Inc., connected to an AutoMate Prep Device, Bushnell, USA). In addition, we tested the used 15N2 gas bottle for contamination with 15N-ammonia62. Briefly, a 2 ml subsample of the used 15N2 gas was injected into a 12 ml gas-tight glass vial (Exetainer) filled with MilliQ (pH 95% sequence identity to any of the MAG NifD/NifK sequences were identified with a blastp search74 to the NCBI nr database. Multiple sequence alignments were obtained with MAFFT87. All full-length sequences were used to construct base trees with RAxML88 and 100 bootstraps in ARB90. The ARB Parsimony function was employed to add partial sequences to the base trees.The resulting trees were visualized in iTOL91.FISH, cell counts, and cell sizesFrom each incubation depth, 10–30 ml lake water was filtered onto 0.22 µm polycarbonate membrane filters (GTTP IsoporeTM, Merck Millipore, USA). The filters were fixed in 2% formaldehyde solution in sterile-filtered lake water for 10–12 h at 4 °C and then washed with MilliQ water. The filters were frozen and stored at −20 °C until further processing.The 16S rRNA FISH probe “Thiosyn459” (Table S7), exclusively targeting T. syntrophicum Cad16, was designed in ARB90. In addition, two competitor probes and four helper probes92 were designed (Table S7) to ensure efficient and specific binding of the probe to the target. All FISH probes and respective formamide concentrations are listed in Table S7. Probes, but not helpers and competitors, were double-labeled with either Atto488 or Atto594 fluorophores. Samples were embedded in 0.05% low melting point agarose. Cells were permeabilized with lysozyme (1.5 mg ml−1) for 30 min at 37 °C. Hybridization was performed for 2–4.5 h at 46 °C. Washing included 15 min in washing buffer at 48 °C and 20 min in 1× PBS buffer at room temperature. We used the hybridization and washing buffers described in Barrero-Canosa et al.93 to reduce background fluorescence. Cells were counterstained with DAPI.Samples were analyzed using a Zeiss Axio Imager.M2 microscope equipped with a Zeiss Axiocam 506 mono camera. Z-stack images were taken and the number of fluorescently labeled cells per image was counted for the individual probes. For each PSB population, we analyzed ≥38 randomly selected fields of view and ≥54 target cells, on one filter replicate each (see Supplementary File S1). Total cell counts were obtained in triplicates through flow cytometry as described in Danza et al.94.For cluster-forming organisms (Thiodictyon syntrophicum, Lamprocystis purpurea, Lamprocystis roseopersicina, and Lamprocystis spp.), the cell size (length and width, for biovolume and C-content calculations, see section below) of 100 cells per population was determined from the maximum-intensity projection of the z-stack images using the Zeiss Zen blue software 3.2.Single-cell analysis with nanoSIMSFor nanoSIMS analyses, we chose the replicate sample from 13.7 m depth that exhibited the highest bulk N2 fixation rate. Random spots were marked with a laser microdissection microscope (6000 B, Leica) on the gold-sputtered GTTP filter covered with cells incubated with 15N2 and 13CO2. After laser marking, FISH was performed as described above. For analysis of Thiodictyon cells, no permeabilization was performed, while for analysis of the other population’s permeabilization was reduced to 15 min at 37 °C using 2 mg ml−1 Lysozyme. Within one hybridization reaction, we simultaneously applied Apur453 with S453D and Laro453 with Cmok453, each probe double labeled with different fluorescent dyes (Atto488 and Atto594).Single-cell 15N- and 13C-assimilation from incubation experiments with 15N2 and 13CO2 was measured using a nanoSIMS 50 L instrument (CAMECA), as described in Martínez-Pérez et al.53. Briefly, instrument precision was monitored regularly on graphite planchet. Samples were pre-sputtered with a Cs+ beam (~300 pA) before the measurements with a beam current of around 1.5 pA. The diameter of the primary beam was tuned More

Modeling frameworkWe consider S species distributed in L distinct habitat patches. The patches form a linear latitudinal chain going around the globe, with dispersal between adjacent patches (Fig. 1). The state variables are species’ local densities and local temperature optima (the temperature at which species achieve maximum intrinsic population growth). This temperature optimum is a trait whose evolution is governed by quantitative genetics18,19,20,21,22: each species, in every patch, has a normally distributed temperature optimum with a given mean and variance. The variance is the sum of a genetic and an environmental contribution. The genetic component is given via the infinitesimal model23,24, whereby a very large number of loci each contribute a small additive effect to the trait. This has two consequences. First, a single round of random mating restores the normal shape of the trait distribution, even if it is distorted by selection or migration. Second, the phenotypic variance is unchanged by these processes, with only the mean being affected25 (we apply a reduction in genetic variance at very low population densities to prevent such species from evolving rapidly; see the Supplementary Information [SI], Section 3.4). Consequently, despite selection and the mixing of phenotypes from neighboring patches, each species retains a normally-shaped phenotypic distribution with the same phenotypic variance across all patches—but the mean temperature optimum may evolve locally and can therefore differ across patches (Fig. 1).Fig. 1: Illustration of our modeling framework.There are several patches hosting local communities, arranged linearly along a latitudinal gradient. Patch color represents the local average temperature, with warmer colors corresponding to higher temperatures. The graph depicts the community of a single patch, with four species present. They are represented by the colored areas showing the distributions of their temperature optima, with the area under each curve equal to the population density of the corresponding species. The green species is highlighted for purposes of illustration. Each species has migrants to adjacent patches (independent of local adaptedness), as well as immigrants from them (arrows from and to the green species; the distributions with dashed lines show the trait distributions of the green species’ immigrant individuals). The purple line is the intrinsic growth rate of a phenotype in the patch, as a function of its local temperature optimum (this optimum differs across patches, which is why the immigrants are slightly maladapted to the temperature of the focal patch.) Both local population densities and local adaptedness are changed by the constant interplay of temperature-dependent intrinsic growth, competition with other species in the same patch, immigration to or emigration from neighboring patches, and (in certain realizations of the model) pressure from consumer species.Full size imageSpecies in our setup may either be resources or consumers. Their local dynamics are governed by the following processes. First, within each patch, we allow for migration to and from adjacent patches (changing both local population densities and also local adaptedness, due to the mixing of immigrant individuals with local ones). Second, each species’ intrinsic rate of increase is temperature-dependent, influenced by how well their temperature optima match local temperatures (Fig. 2a). For consumers, metabolic loss and mortality always result in negative intrinsic growth, which must be compensated by sufficient consumption to maintain their populations. Third, there is a local competition between resource species, which can be thought of as exploitative competition for a set of shared substitutable lower-level resources26. Consumers, when present, compete only indirectly via their shared resource species. Fourth, each consumer has feeding links to five of the resource species (pending their presence in patches where the consumer is also present), which are randomly determined but always include the one resource which matches the consumer’s initial mean temperature optimum. Feeding rates follow a Holling type II functional response. Consumers experience growth from consumption, and resource species experience loss due to being consumed.Fig. 2: Temperature optima and climate curves.a Different growth rates at various temperatures. Colors show species with different mean temperature optima, with warmer colors corresponding to more warm-adapted species. The curves show the maximum growth rate achieved when a phenotype matches the local temperature, and how the growth rate decreases with an increased mismatch between a phenotype and local temperature, for each species. The dashed line shows zero growth: below this point, the given phenotype of a species mismatches the local temperature to the extent that it is too maladapted to be able to grow. Note the tradeoff between the width and height of the growth curves, with more warm-tolerant species having larger maximum growth at the cost of being viable for only a narrower range of temperatures62,63. b Temperature changes over time. After an initial establishment phase of 4000 years during which the pre-climate change community dynamics stabilize, temperatures start increasing at t = 0 for 300 years (vertical dotted line, indicating the end of climate change). Colors show temperature change at different locations along the spatial gradient, with warmer colors indicating lower latitudes. The magnitude and latitudinal dependence of the temperature change is based on region-specific predictions by 2100 CE, in combination with estimates giving an approximate increase by 2300 CE, for the IPCC intermediate emission scenario27.Full size imageFollowing the previous methodology, we derive our equations in the weak selection limit22 (see also the Discussion). We have multiple selection forces acting on the different components of our model. Species respond to local climate (frequency-independent directional selection, unless a species is at the local environmental optimum), to consumers and resources (frequency-dependent selection), and competitors (also frequency-dependent selection, possibly complicated by the temperature-dependence of the competition coefficients mediating frequency dependence). These different modes of selection do not depend on the parameterization of evolution and dispersal, which instead are used to adjust the relative importance of these processes.Communities are initiated with 50 species per trophic level, subdividing the latitudinal gradient into 50 distinct patches going from pole to equator (results are qualitatively unchanged by increasing either the number of species or the number of patches; SI, Section 5.9–5.10). We assume that climate is symmetric around the equator; thus, only the pole-to-equator region needs to be modeled explicitly (SI, Section 3.5). The temperature increase is based on predictions from the IPCC intermediate emission scenario27 and corresponds to predictions for the north pole to the equator. The modeled temperature increase is represented by annual averages and the increase is thus smooth. Species are initially equally spaced, and adapted to the centers of their ranges. We then integrate the model for 6500 years, with three main phases: (1) an establishment period from t = −4000 to t = 0 years, during which local temperatures are constant; (2) climate change, between t = 0 and t = 300 years, during which local temperatures increase in a latitude-specific way (Fig. 2b); and (3) the post-climate change period from t = 300 to t = 2500 years, where temperatures remain constant again at their elevated values.To explore the influence and importance of dispersal, evolution, and interspecific interactions, we considered the fully factorial combination of high and low average dispersal rates, high and low average available genetic variance (determining the speed and extent of species’ evolutionary responses), and four different ecological models. These were: (1) the baseline model with a single trophic level and constant, patch- and temperature-independent competition between species; (2) two trophic levels and constant competition; (3) single trophic level with temperature-dependent competition (where resource species compete more if they have similar temperature optima); and (4) two trophic levels as well as temperature-dependent competition. Trophic interactions can strongly influence diversity in a community, either by apparent competition28 or by acting as extra regulating agents boosting prey coexistence29. Temperature-dependent competition means that the strength of interaction between two phenotypes decreases with an increasing difference in their temperature optima. Importantly, while differences in temperature adaptation may influence competition, they do not influence trophic interactions.The combination of high and low genetic variance and dispersal rates, and four model setups, gives a total of 2 × 2 × 4 = 16 scenarios. For each of them, some parameters (competition coefficients, tradeoff parameters, genetic variances, dispersal rates, consumer attack rates, and handling times; SI, Section 6) were randomly drawn from pre-specified distributions. We, therefore, obtained 100 replicates for each of these 16 scenarios. While replicates differed in the precise identity of the species which survived or went extinct, they varied little in the overall patterns they produced.We use the results from these numerical experiments to explore patterns of (1) local species diversity (alpha diversity), (2) regional trends, including species range breadths and turnover (beta diversity), (3) global (gamma) diversity, and global changes in community composition induced by climate change. In addition, we also calculated the interspecific community-wide trait lag (the difference between the community’s density-weighted mean temperature optima and the current temperature) as a function of the community-wide weighted trait dispersion (centralized variance in species’ density-weighted mean temperature optima; see Methods). The response capacity is the ability of the biotic community to close this trait lag over time30 (SI, Section 4). Integrating trait lag through time31 gives an overall measure of different communities’ ability to cope with changing climate over this time period; furthermore, this measure is comparable across communities. The integrated trait lag summarizes, in a single functional metric, the performance and adaptability of a community over space and time. The reason it is related to performance is that species that on average live more often under temperatures closer to their optima (creating lower trait lags) will perform better than species whose temperature optima are far off from local conditions in space and/or time. Thus, a lower trait lag (higher response capacity) may also be related to other ecosystem functions, such as better carbon uptake which in turn has the potential to feedback to global temperatures32.Overview of resultsWe use our framework to explore the effect of species interactions on local, regional, and global biodiversity patterns, under various degrees of dispersal and available genetic variance. For simplicity, we focus on the dynamics of the resource species, which are present in all scenarios. Results for consumers, when present, are in the SI (Section 5.8). First, we display a snapshot of species’ movement across the landscape with time; before, during, and after climate change. Then we proceed with analyzing local patterns, followed by regional trends, and finally, global trends.Snapshots from the time series of species’ range distributions reveal useful information about species’ movement and coexistence (Fig. 3). Regardless of model setup and parameterization, there is a northward shift in species’ ranges: tropical species expand into temperate regions and temperate species into polar regions. This is accompanied by a visible decline in the number of species globally, with the northernmost species affected most. The models do differ in the predicted degree of range overlap: trophic interactions and temperature-dependent competition both lead to broadly overlapping ranges, enhancing local coexistence (the overlap in spatial distribution is particularly pronounced with high available genetic variance). Without these interactions, species ranges overlap to a substantially lower degree, diminishing local diversity. Below we investigate whether these patterns, observed for a single realization of the dynamics for each scenario, play out more generally as well.Fig. 3: Species’ range shift through time, along a latitudinal gradient ranging from polar to tropical climates (ordinate).Species distributions are shown by colored curves, with the height of each curve representing local density in a single replicate (abscissa; note the different scales in the panels), with the color indicating the species’ initial (i.e., at t = 0) temperature adaptation. The model was run with only 10 species, for better visibility. The color of each species indicates its temperature adaptation at the start of the climate change period, with warmer colors belonging to species with a higher temperature optimum associated with higher latitudes. Rows correspond to a specific combination of genetic variance and dispersal ability of species, columns show species densities at different times (t = 0 start of climate change, t = 300 end of climate change, t = 2500 end of simulations). Each panel corresponds to a different model setup; a the baseline model, b an added trophic level of consumers, c temperature-dependent competition coefficients, and d the combined influence of consumers and temperature-dependent competition.Full size imageLocal trendsTrophic interactions and temperature-dependent competition indeed result in elevated local species richness levels (Fig. 4). The fostering of local coexistence by trophic interactions and temperature-dependent competition is in line with general ecological expectations. Predation pressure can enhance diversity by providing additional mechanisms of density regulation and thus prey coexistence through predator partitioning28,29. In turn, temperature-dependent competition means species can reduce interspecific competition by evolving locally suboptimal mean temperature optima22, compared with the baseline model’s fixed competition coefficients. Hence with temperature-dependent competition, the advantages of being sufficiently different from other locally present species can outweigh the disadvantages of being somewhat maladapted to the local temperatures. If competition is not temperature-dependent, interspecific competition is at a fixed level independent of the temperature optima of each species. An important question is how local diversity is affected when the two processes act simultaneously. In fact, any synergy between their effects is very weak, and is even slightly negative when both the available genetic variance and dispersal abilities are high (Fig. 4, top row).Fig. 4: Local species richness of communities over time, from the start of climate change to the end of the simulation, averaged over replicates.Values are given in 100-year steps. At each point in time, the figure shows the mean number of species per patch over the landscape (points) and their standard deviation (shaded region, extending one standard deviation both up- and downwards from the mean). Panel rows show different parameterizations (all four combinations of high and low genetic variance and dispersal ability); columns represent various model setups (the baseline model; an added trophic level of consumers; temperature-dependent competition coefficients; and the combined influence of consumers and temperature-dependent competition). Dotted vertical lines indicate the time at which climate change ends.Full size imageRegional trendsWe see a strong tendency for poleward movement of species when looking at the altered distributions of species over the spatial landscape (Fig. 3). Indeed, looking at the effects of climate change on the fraction of patches occupied by species over the landscape reveals that initially cold-adapted species lose suitable habitat during climate change, and even afterwards (Fig. 5). For the northernmost species, this always eventuate to the point where all habitat is lost, resulting in their extinction. This pattern holds universally in every model setup and parameterization. Only initially warm-adapted species can expand their ranges, and even they only do so under highly restrictive conditions, requiring both good dispersal ability and available genetic variance as well as consumer pressure (Fig. 5, top row, second and third panel).Fig. 5: Range breadth of each species expressed as the percentage of the whole landscape they occupy (ordinate) at three different time stamps (colors).The mean (points) and plus/minus one standard deviation range (colored bands) are shown over replicates. Numbers along the abscissa represent species, with initially more warm-adapted species corresponding to higher values. The range breadth of each species is shown at three time stamps: at the start of climate change (t = 0, blue), the end of climate change (t = 300, green), and at the end of our simulations (t = 2500, yellow). Panel layout as in Fig. 4.Full size imageOne can also look at larger regional changes in species richness, dividing the landscape into three equal parts: the top third (polar region), the middle third (temperate region), and the bottom third (tropical region). Region-wise exploration of changes in species richness (Fig. 6) shows that the species richness of the polar region is highly volatile. It often experiences the greatest losses; however, with high dispersal ability and temperature-dependent competition, the regional richness can remain substantial and even increase compared to its starting level (Fig. 6, first and third rows, last two columns). Of course, change in regional species richness is a result of species dispersing to new patches and regions as well as of local extinctions. Since the initially most cold-adapted species lose their habitat and go extinct, altered regional species richness is connected to having altered community compositions along the spatial gradient. All regions experience turnover in species composition (SI, Section 5.1), but in general, the polar region experiences the largest turnover, where the final communities are at least 50% and sometimes more than 80% dissimilar to the community state right before the onset of climate change—a result in agreement with previous studies as well7,33.Fig. 6: Relative change in global species richness from the community state at the onset of climate change (ordinate) over time (abscissa), averaged over replicates and given in 100-year steps (points).Black points correspond to species richness over the whole landscape; the blue points to richness in the top third of all patches (the polar region), green points to the middle third (temperate region), and yellow points to the last third (tropical region). Panel layout as in Fig. 4; dotted horizontal lines highlight the point of no net change in global species richness.Full size imageGlobal trendsHence, the identity of the species undergoing global extinction is not random, but strongly biased towards initially cold-adapted species. On a global scale, these extinctions cause decreased richness, and the model predicts large global biodiversity losses for all scenarios (Fig. 6). These continue during the post-climate change period with stable temperatures, indicating a substantial extinction debt which has been previously demonstrated34. Temperature-dependent competition reduces the number of global losses compared to the baseline and trophic models.A further elucidating global pattern is revealed by analyzing the relationship between the time-integrated temperature trait lag and community-wide trait dispersion (Fig. 7). There is an overall negative correlation between the two, but more importantly, within each scenario (unique combination of model and parameterization) a negative relationship is evident. Furthermore, the slopes are very similar: the main difference between scenarios is in their mean trait lag and trait dispersion values (note that the panels do not share axis value ranges). The negative trend reveals the positive effect of more varied temperature tolerance strategies among the species on the community’s ability to respond to climate change. This is analogous to Fisher’s fundamental theorem35, stating that the speed of the evolution of fitness r is proportional to its variance: dr/dt ~ var(r). More concretely, this relationship is also predicted by trait-driver theory, a mathematical framework that focuses explicitly on linking spatiotemporal variation in environmental drivers to the resulting trait distributions30. Communities generated by different models reveal differences in the magnitude of this relationship: trait dispersion is much higher in models with temperature-dependent competition (essentially, niche differentiation with respect to temperature), resulting in lower trait lag. The temperature-dependent competition also separates communities based on their spatial dispersal ability, with faster dispersal corresponding to greater trait dispersion and thus lower trait lag. Interestingly, trophic interactions tend to erode the relationship between trait lag and trait dispersion slightly (R2 values are lower in communities with trophic interactions, both with and without temperature-dependent competition). We have additionally explored the relationship between species richness and trait dispersion, finding a positive relationship between the two (SI, Section 4.1).Fig. 7: The ability of communities in four different models (panels) to track local climatic conditions (ordinate), against observed variation in traits within those communities (abscissa).Larger values along the ordinate indicate that species’ temperature optima are lagging behind local temperatures, meaning a low ability of communities to track local climate conditions. Both quantities are averaged over the landscape and time from the beginning to the end of the climate change period, yielding a single number for every community (points). The greater the average local diversity of mean temperature optima in a community, the closer it is able to match the prevailing temperature conditions. Species’ dispersal ability and available genetic variance (colors) are clustered along this relationship.Full size image More

1.Geist, J. Integrative freshwater ecology and biodiversity conservation. Ecol. Indic. 11, 1507–1516 (2011).Article 

Google Scholar 
2.Lopes-Lima, M. et al. Conservation status of freshwater mussels in Europe: State of the art and future challenges. Biol. Rev. 92, 572–607 (2017).PubMed 
Article 
PubMed Central 

Google Scholar 
3.Geist, J. Strategies for the conservation of endangered freshwater pearl mussels (Margaritifera margaritifera L.): A synthesis of conservation genetics and ecology. Hydrobiologia 644, 69–88 (2010).Article 

Google Scholar 
4.Taeubert, J. E. & Geist, J. The relationship between the freshwater pearl mussel (Margaritifera margaritifera) and its hosts. Biol. Bull. 44, 67–73 (2017).Article 

Google Scholar 
5.Salonen, J. K. et al. Atlantic salmon (Salmo salar) and brown trout (Salmo trutta) differ in their suitability as a host for the endangered freshwater pearl mussel (Margaritifera margaritifera) in northern Fennoscandian rivers. Freshw. Biol. 62, 1346–1358 (2017).Article 

Google Scholar 
6.Geist, J. & Auerswald, K. Physicochemical stream bed characteristics and recruitment of the freshwater pearl mussel (Margaritifera margaritifera). Freshw. Biol. 52, 2299–2316 (2007).Article 

Google Scholar 
7.Stoeckl, K., Denic, M. & Geist, J. Conservation status of two endangered freshwater mussel species in Bavaria, Germany: Habitat quality, threats, and implications for conservation management. Aquat. Conserv. 30, 647–661 (2020).Article 

Google Scholar 
8.Auerswald, K. & Geist, J. Extent and cause of siltation in a headwater stream bed: Catchment and soil erosion is less important than internal stream processes. Land Degrad. Dev. 29, 737–748. https://doi.org/10.1002/ldr.2779 (2018).Article 

Google Scholar 
9.Bauer, G. Threats to the freshwater pearl mussel Margaritifera margaritifera L. in central Europe. Biol. Conserv. 45, 239–253 (1988).Article 

Google Scholar 
10.Boon, P. J. et al. Developing a standard approach for monitoring freshwater pearl mussel (Margaritifera margaritifera) populations in European rivers. Aquat. Conserv. 29, 1365–1379 (2019).Article 

Google Scholar 
11.Hruska, J. Nahrungsansprüche der Flußperlmuschel und deren halbnatürliche Aufzucht in der Tschechischen Republik (Dietary requirements and semi-natural rearing of freshwater pearl mussel in the Czech Republic). Heldia 4, 69–79 (1999).
Google Scholar 
12.Preston, S. J., Keys, A. & Roberts, D. Culturing freshwater pearl mussel Margaritifera margaritifera: A breakthrough in the conservation of an endangered species. Aquat. Conserv. 17, 539–549. https://doi.org/10.1002/aqc.799 (2007).Article 

Google Scholar 
13.Thomas, G. R., Taylor, J. & de Leaniz, C. G. Captive breeding of the endangered freshwater pearl mussel, Margaritifera margaritifera. Endanger. Species Res. 12, 1–9 (2010).Article 

Google Scholar 
14.Gum, B., Lange, M. & Geist, J. A critical reflection on the success of rearing and culturing juvenile freshwater mussels with a focus on the endangered freshwater pearl mussel (Margaritifera margaritifera L.). Aquat. Conserv. 21, 743–751 (2011).Article 

Google Scholar 
15.Geist, J., Rottmann, O., Schröder, W. & Kühn, R. Development of microsatellite markers for the endangered freshwater pearl mussel Margaritifera margaritifera L. (Bivalvia: Unionoidea). Mol. Ecol. Resour. 3, 444–446 (2003).CAS 
Article 

Google Scholar 
16.Geist, J. & Kühn, R. Genetic diversity and differentiation of central European freshwater pearl mussel (Margaritifera margaritifera L.) populations: Implications for conservation and management. Mol. Ecol. 14, 425–439 (2005).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
17.Geist, J. & Kuehn, R. Host-parasite interactions in oligotrophic stream ecosystems: The roles of life history strategy and ecological niche. Mol. Ecol. 17, 997–1008 (2008).PubMed 
Article 
PubMed Central 

Google Scholar 
18.Marchordom, A., Araujo, R., Erpenbeck, D. & Ramos, M. A. Phylogeography and conservation genetics of the endangered European Margaritiferidae (Bivalvia: Unionoidea). Biol. J. Linn. Soc. Lond. 78, 235–252 (2003).Article 

Google Scholar 
19.Stoeckle, et al. Strong genetic differentiation and low genetic diversity of the freshwater pearl mussel (Margaritifera margaritifera L.) in the southwestern European distribution range. Conserv. Genet. 18, 147–157 (2017).Article 

Google Scholar 
20.Karlsson, S., Larsen, B. M. & Hindar, K. Host-dependent genetic variation in freshwater pearl mussel (Margaritifera margaritifera L.). Hydrobiologia 735, 179–190 (2014).Article 

Google Scholar 
21.Geist, J., Söderberg, H., Karlberg, A. & Kuehn, R. Drainage-independent genetic structure and high genetic diversity of endangered freshwater pearl mussels (Margaritifera margaritifera) in northern Europe. Conserv. Genet. 11, 1339–1350 (2010).Article 

Google Scholar 
22.Geist, et al. Genetic structure of Irish freshwater pearl mussels (Margaritifera margaritifera and Margaritifera durrovensis): Validity of subspecies, roles of host fish, and conservation implications. Aquat. Conserv. 28, 923–933 (2018).Article 

Google Scholar 
23.Zanatta, et al. High genetic diversity and low differentiation in North American Margaritifera margaritifera (Bivalvia: Unionida: Margaritiferidae). Biol. J. Linn. Soc. Lond. 123, 850–863 (2018).Article 

Google Scholar 
24.Taeubert, J. E., Denic, M., Gum, B., Lange, M. & Geist, J. Suitability of different salmonid strains as hosts for the endangered freshwater pearl mussel (Margaritifera margaritifera). Aquat. Conserv. 20, 728–734 (2010).Article 

Google Scholar 
25.Marwaha, et al. Host (Salmo trutta) age influences resistance to infestation by freshwater pearl mussel (Margaritifera margaritifera) glochidia. Parasitol. Res. 118, 1519–1532 (2019).PubMed 
Article 

Google Scholar 
26.Taeubert, J. E., Gum, B. & Geist, J. Variable development and excystment of freshwater pearl mussel (Margaritifera margaritifera L.) at constant temperature. Limnologica 43, 319–322 (2013).Article 

Google Scholar 
27.Taeubert, J. E. & Geist, J. Critical swimming speed of brown trout (Salmo trutta) infested with freshwater pearl mussel (Margaritifera margaritifera) glochidia and implications for artificial breeding of an endangered mussel species. Parasitol. Res. 112, 1607–1613 (2013).PubMed 
Article 

Google Scholar 
28.Marwaha, J., Jensen, K. H., Jakobsen, P. J. & Geist, J. Duration of the parasitic phase determines subsequent performance in juvenile freshwater pearl mussels (Margaritifera margaritifera). Ecol. Evol. 7, 1375–1383 (2017).PubMed 
PubMed Central 
Article 

Google Scholar 
29.Eybe, T., Thielen, F., Bohn, T. & Sures, B. Influence of the excystment time on the breeding success of juvenile freshwater pearl mussels (Margaritifera margaritifera). Aquat. Conserv. 25, 21–30 (2015).Article 

Google Scholar 
30.Denic, M., Taeubert, J. E. & Geist, J. Trophic relationships between the larvae of two freshwater mussels and their fish hosts. Invertebr. Biol. 134, 129–135 (2015).Article 

Google Scholar 
31.Denic, M. et al. Influence of stock origin and environmental conditions on the survival and growth of juvenile freshwater pearl mussels (Margaritifera margaritifera) in a cross-exposure experiment. Limnologica 50, 67–74 (2015).CAS 
Article 

Google Scholar 
32.Hyvärinen, H. S. H., Chowdhury, M. M. R. & Taskinen, J. Pulsed flow-through cultivation of Margaritifera margaritifera: Effects of water source and food quantity on the survival and growth of juveniles. Hydrobiologia. 3219–3229 (2021).33.Hyvärinen, H., Saarinen-Valta, M., Mäenpää, E. & Taskinen, J. Effect of substrate particle size on burrowing of the juvenile freshwater pearl mussel Margaritifera margaritifera. Hydrobiologia https://doi.org/10.1007/s10750-021-04522-z (2021).Article 

Google Scholar 
34.Taskinen, J. et al. Effect of pH, iron and aluminum on survival of early life history stages of the endangered freshwater pearl mussel, Margaritifera margaritifera. Toxicol. Environ. Chem. 93, 1764–1777 (2011).CAS 
Article 

Google Scholar 
35.Lavictoire, L., Moorkens, E., Ramsay, A. & Sweeting, R. Effects of substrate size and cleaning regime on growth and survival of captive-bred juvenile freshwater pearl mussels, Margaritifera margaritifera (Linnaeus, 1758). Hydrobiologia 766, 89–102 (2016).Article 

Google Scholar 
36.Eybe, T., Thielen, F., Bohn, T. & Sures, B. The first millimetre: Rearing juvenile freshwater pearl mussels (Margaritifera margaritifera L.) in plastic boxes. Aquat. Conserv. 23, 964–975 (2013).Article 

Google Scholar 
37.Strayer, D. L., Geist, J., Haag, W. R., Jackson, J. K. & Newbold, J. D. Essay: Making the most of recent advances in freshwater mussel propagation and restoration. Conserv. Sci. Pract. 1, e53. https://doi.org/10.1111/csp2.53 (2019).Article 

Google Scholar 
38.Patterson, M. A. et al. Freshwater Mussel Propagation for Restoration (Cambridge University Press, 2018).Book 

Google Scholar 
39.Gstöttenmayr, D., Scheder, C. & Gumpinger, C. Conservation de la mulette perlière d’eau douce en Autriche: un système d’élevage contrôlé en progrès. Penn ar Bed 222, 45–49 (2015).
Google Scholar 
40.Gumpinger, C., Pichler-Scheder, C. & Huemer, D. Das oberösterreichische Artenschutzprojekt „Vision Flussperlmuschel“. Österreichs Fischerei 69, 259–273 (2016).
Google Scholar 
41.Rice, W. R. Analyzing tables of statistical tests. Evolution 43, 223–225 (1989).PubMed 
Article 
PubMed Central 

Google Scholar 
42.DeWoody, J. A. et al. Universal method for producing ROXlabeled size standards suitable for automated genotyping. Biotechniques 37, 348–352. https://doi.org/10.2144/04373BM02 (2004).CAS 
Article 
PubMed 
PubMed Central 

Google Scholar 
43.Goudet, J. Fstat (Version 1.2): A computer program to calculate F-statistics. J. Hered. 86, 485–486 (1995).44.Rousset, F. Genepop’007: A complete reimplementation of the Genepop software for Windows and Linux. Mol. Ecol. Resour. 8, 103–106 (2008).PubMed 
Article 
PubMed Central 

Google Scholar 
45.Haldane, J. B. S. An exact test for randomness of mating. J. Genet. 52, 631–635. https://doi.org/10.1007/BF02981502 (1954).Article 

Google Scholar 
46.Weir, B. S. & Cockerham, C. C. Estimating F-statistics for the analysis of population structure. Evolution 38, 1358–1370 (1984).CAS 
PubMed 
PubMed Central 

Google Scholar 
47.Guo, S. W. & Thompson, E. A. Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics 48, 361–372 (1992).CAS 
PubMed 
MATH 
Article 
PubMed Central 

Google Scholar 
48.Raymond, M. & Rousset, F. An exact test for population differentiation. Evolution 49, 1280–1283 (1995).PubMed 
Article 
PubMed Central 

Google Scholar 
49.Kamvar, Z. N., Tabima, J. F. & Grünwald, N. J. Poppr: An R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ 2, e281. https://doi.org/10.7717/peerj.281 (2014).Article 
PubMed 
PubMed Central 

Google Scholar 
50.Jombart, T. Adegenet: A R package for the multivariate analysis of genetic markers. Bioinformatics 24, 1403–1405 (2008).CAS 
PubMed 
Article 

Google Scholar 
51.Ciofi, C., Beaumont, M. A., Swingland, I. R. & Bruford, M. W. Genetic divergence and units for conservation in the Komodo dragon Varanus komodoensis. Proc. Royal Soc. B 266, 2269–2274 (1999).Article 

Google Scholar 
52.Pritchard, J. K., Stephens, M. & Donnelly, P. Inference of population structure using multilocus genotype data. Genetics 155, 945–959 (2000).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
53.Kopelman, N. M., Mayzel, J., Jakobsson, M., Rosenberg, N. A. & Mayrose, I. CLUMPAK: A program for identifying clustering modes and packaging population structure inferences across K. Mol. Ecol. Resour. 15, 1179–1191 (2015).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
54.Evanno, G., Regnaut, S. & Goudet, J. Detecting the number of clusters of individuals using the software structure: A simulation study. Mol. Ecol. 14, 2611–2620 (2005).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
55.Jakobsson, M. & Rosenberg, N. A. CLUMPP: A cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics 23, 1801–1806 (2007).CAS 
PubMed 
PubMed Central 
Article 

Google Scholar 
56.Rosenberg, N. A. DISTRUCT: A program for the graphical display of population structure. Mol. Ecol. Notes 4, 137–138 (2004).Article 

Google Scholar 
57.Kalinowski, S. T. The computer program STRUCTURE does not reliably identify the main genetic clusters within species: Simulations and implications for human population structure. Heredity 106, 625–632 (2011).CAS 
PubMed 
Article 
PubMed Central 

Google Scholar 
58.Puechmaille, S. J. The program structure does not reliably recover the correct population structure when sampling is uneven: Subsampling and new estimators alleviate the problem. Mol. Ecol. Resour. 16, 608–627 (2016).PubMed 
Article 
PubMed Central 

Google Scholar 
59.Jombart, T., Devillard, S. & Balloux, F. Discriminant analysis of principal components: A new method for the analysis of genetically structured populations. BMC Genet. 11, 94 (2010).PubMed 
PubMed Central 
Article 

Google Scholar 
60.R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. URL: https://www.R-project.org/ (2019).61.Trushenski, J. T., Whelan, G. E. & Bowker, J. D. Why keep hatcheries? Weighing the economic cost and value of fish production for public use and public trust purposes. Fisheries 43, 285–293 (2018).
Google Scholar 
62.Wacker, S., Larsen, B. M., Jakobsen, P. & Karlsson, S. High levels of multiple paternity in a spermcast mating freshwater mussel. Ecol. Evol. 8, 8126–8134 (2018).PubMed 
PubMed Central 
Article 

Google Scholar 
63.Wacker, S., Larsen, B. M., Jakobsen, P. & Karlsson, S. Multiple paternity promotes genetic diversity in captive breeding of a freshwater mussel. Glob. Ecol. Conserv. 17, e00564. https://doi.org/10.1016/j.gecco.2019.e00564 (2019).Article 

Google Scholar 
64.Garrison, N. L., Johnson, P. D. & Whelan, N. V. Conservation genomics reveals low genetic diversity and multiple parentage in the threatened freshwater mussel, Margaritifera hembeli. Conserv. Genet. https://doi.org/10.1007/s10592-020-01329-8 (2021).Article 

Google Scholar 
65.Bauer, G. Reproductive strategy of the freshwater pearl mussel Margaritifera margaritifera. J. Anim. Ecol. 56, 691–704 (1987).Article 

Google Scholar 
66.McMurray, S. E. & Roe, K. J. Perspectives on the controlled propagation, augmentation, and reintroduction of freshwater mussels (Mollusca: Bivalvia: Unionoida). Freshw. Mollusk Biol. Conserv. 20, 1–12 (2017).Article 

Google Scholar 
67.Geist, J. Seven steps towards improving freshwater conservation. Aquat. Conserv. 25, 447–453 (2015).Article 

Google Scholar  More