OSCR

Recurrent synapses between CO2-sensitive olfactory sensory neurons enable robust CO2 detection in Aedes aegypti mosquitoes.

Code ↔ Paper

4 matches between paragraphs of the paper and lines of its authors' code, computed by the harvester (lexical-v1). Click a colored paragraph or line to see its counterpart.

The 4 matches
  1. [1] § Methods › Analysis of D. melanogaster neurons in Flywire ↔ retrieve data/retrieve_neuron_data.ipynb, lines 325–414 · score 0.59 · root_id, DM1, hemibrain, VC3, VC5, VM6
  2. [2] § Methods › Network model ↔ model/ReciprocalOSN-model.ipynb, lines 152–169 · score 0.56 · odor responses, background odors, root, model, OSN
  3. [3] § Methods › Quantification and statistical analysis ↔ fig_s7_s8/cld_lettering.py, lines 20–103 · score 0.53 · pairwise comparisons, Piepho, algorithm, median, S8, S7
  4. [4] § Methods › Gene expression of Acetylcholine and GABA in olfactory sensory neurons ↔ fig_s6/fig_s6b_gaba_gene_expression.ipynb, lines 115–119 · score 0.52 · gene expression, GABA, raw, antennae, Cell, neurons

Paper

Loaded from Europe PMC by your browser, not stored by OSCR: doi.org · Europe PMC

The paper is loaded when this pane is shown.

The authors' code

Jupyter notebook · 725 lines · 33 KB · GPL-3.0 · 1 match

  1. # %%
  2. import pandas as pd
  3. import matplotlib.pyplot as plt
  4. from fafbseg import flywire
  5. import pymaid
  6. import navis as nv
  7. import numpy as np
  8. import seaborn as sns
  9. import scipy.stats as stats
  10. import statsmodels
  11. import scikit_posthocs as sp
  12. import sys
  13. import scipy
  14. import random
  15. import string
  16. from sklearn import metrics
  17. from sklearn.cluster import KMeans
  18. from sklearn.preprocessing import StandardScaler
  19. from sklearn.cluster import DBSCAN
  20. import matplotlib.colors as colors
  22. flywire.get_materialization_versions(dataset="production")
  24. # %%
  25. #check available flywire versions
  26. flywire.get_materialization_versions(dataset="production")
  27. #connect your catmaid instance
  28. catmaid_token = ""
  29. instance=pymaid.CatmaidInstance('https://radagast.hms.harvard.edu/catmaidaedes', catmaid_token)
  31. # %%
  32. def get_cable_overlap(OSNs, PN = None, volume = None):
  33. neurons = flywire.skeletonize_neuron(OSNs)
  34. if volume:
  35. neurons = nv.in_volume(neurons, volume)
  37. if PN == None:
  38. cable_overlap = nv.cable_overlap(neurons, neurons, dist= "2 microns")
  39. return cable_overlap
  41. else:
  42. PN = flywire.skeletonize_neuron(PN)
  43. PN_cable_overlap = nv.cable_overlap(neurons, PN, dist= "2 microns")
  44. return PN_cable_overlap
  46. def get_cable_length(OSNs,volume = None):
  47. neurons = flywire.skeletonize_neuron(OSNs)
  48. if volume:
  49. neurons = nv.in_volume(neurons, volume)
  51. return neurons.cable_length
  53. def get_total_synapses(OSNs, volume = None):
  54. neurons = flywire.skeletonize_neuron(OSNs)
  55. if volume:
  56. neurons = nv.in_volume(neurons, volume)
  57. return flywire.get_synapses(neurons, pre=True, post=False, attach=True, dataset='production', min_score=65,progress=False)
  59. def get_num_connectors(neuron1, neuron2, volume = None):
  60. if volume:
  61. neuron1 = nv.in_volume(neuron1, volume)
  62. neuron2 = nv.in_volume(neuron2, volume)
  63. df = flywire.synapses.get_adjacency(neuron1, targets=neuron2, dataset='production', min_score=65, progress=True)
  64. return df.values[0][0]
  66. def synapses_per_micron(cable_overlap_df, PN = None):
  67. cable_overlap_synapse_density = []
  68. for rowIndex, row in cable_overlap_df.iterrows(): #iterate over rows
  69. synapse_count = []
  70. cable_overlap_length = []
  71. for columnIndex, value in row.items():
  72. if value != 0:
  73. synapse_count.append(get_num_connectors(rowIndex, columnIndex))
  74. cable_overlap_length.append(value)
  75. cable_overlap_synapse_density.append(np.nanmean(np.array(synapse_count) / np.array(cable_overlap_length)) * 1000)
  76. return pd.DataFrame(cable_overlap_synapse_density,columns=['overlap_density'], index=cable_overlap_df.index)
  78. def get_neurons_without_overlap(cable_overlap_df):
  79. cable_overlap_df = cable_overlap_df.sum(axis=1)
  80. no_PN_overlap_neurons = [cable_overlap_df.index[index] for index, neuron in enumerate(cable_overlap_df) if neuron == 0.0]
  81. return no_PN_overlap_neurons
  83. # %% [markdown]
  84. # ## Projection Neuron IDs for glomeruli
  85. # |Glomerulus|Right Hemisphere|Left Hemisphere|
  86. # |:---------|:------------:|-----------:|
  87. # |DA4l |720575940608970197 |720575940649045369|
  88. # |DC1| 720575940637056887| 720575940621529435|
  89. # |DL4 |720575940627708688 | 720575940617343316|
  90. # |DL5 |720575940617207185 | 720575940639080765|
  91. # |DM3| 720575940614727903 | 720575940630549370|
  92. # |DM4 |720575940623528925 | 720575940615366055|
  93. # |DP1m |720575940618308825 | 720575940622726271|
  94. # |VA2| 720575940611079236 | 720575940609460491|
  95. # |VA4| 720575940617668747 | 720575940642086389|
  96. # |VA6 |720575940636873791 | 720575940637594334|
  97. # |VA7l | 720575940625784153 | 720575940618277184|
  98. # |VC1 |720575940637526190 | 720575940633190847|
  99. # |VC2 |720575940628739560 | 720575940630097103|
  100. # |VL2a |720575940618295454 | 720575940603464672|
  101. # |VL2p |720575940628467611 | 720575940618865112|
  102. # |VM4 |720575940640971995|
  103. # |DA4m |720575940631930828|
  104. # |VM6 |720575940624057853|
  105. #
  106. # ## Glomeruli with more than one uPN projection neurons
  107. # |Glomerulus|Right Hemisphere|
  108. # |:---------|------------:|
  109. # |D |720575940603985952 720575940633804212 720575940627844734|
  110. # |DA1 |720575940614309535 720575940637208718 720575940613345442 720575940626034819 720575940619385765 720575940605102694 720575940603231916 720575940621185050|
  111. # |DA2 |720575940622364184 720575940611849187 720575940626937617 720575940622734835|
  112. # |DC2 |720575940627160322 720575940616824588|
  113. # |DC3| 720575940630131663 720575940613419695 720575940634939327 720575940634715743|
  114. # |DL1 |720575940616079035 720575940622368792 720575940613809554|
  115. # |DL2d |720575940628316047 720575940630928715 720575940620103590 720575940642046453|
  116. # |DL3 |720575940630778428 720575940621926669 720575940627353426 720575940630778428 720575940620638335|
  117. # |DM2 |720575940638633535 720575940630024566|
  118. # |DM5 |720575940644258327 720575940629926159 720575940640032896|
  119. # |DM6 |720575940631328104 720575940622785062 720575940620437339|
  120. # |DP1l |720575940648650884 720575940606762686|
  121. # |VA1d |720575940622811430 720575940640804853 720575940619620852|
  122. # |VA1v |720575940629097922 720575940620199962 720575940628283560 720575940629733626 720575940632720026 720575940634229615|
  123. # |VA3 |720575940656090017 720575940629009513 720575940604431584|
  124. # |VA5 |720575940619838912 720575940637009636 720575940637231593|
  125. # |VA7m |720575940620068518 720575940640398477 720575940639680957|
  126. # |VC3m |720575940608845525 720575940608199748 720575940624515571|
  127. # |VL1 |720575940624863015 720575940637712243 720575940620595053|
  128. # |VM1 |720575940628077568 720575940623444715|
  129. # |VM2 |720575940640697808 720575940617179451|
  130. # |VM3 |720575940627480586 720575940610048981|
  131. # |VM7d |720575940628524292 720575940613167007|
  132. # |VM7v |720575940619337536 720575940621256212|
  133. # |VM5d |720575940611372714 720575940644961588 720575940635082871 720575940627035388|
  135. # %%
  136. DA4l_pn = 720575940608970197
  137. DC1_pn = 720575940637056887
  138. DL4_pn = 720575940627708688
  139. DL5_pn = [720575940617207185, 720575940630432815]
  140. DM3_pn = 720575940614727903
  141. DM4_pn = [720575940623528925 , 720575940628323919, 720575940606767806]
  142. DP1m_pn = 720575940618308825
  143. VA2_pn = 720575940611079236
  144. VA4_pn = 720575940617668747
  145. VA6_pn = 720575940636873791
  146. VA7l_pn = 720575940625784153
  147. VC1_pn = 720575940637526190
  148. VC2_pn = 720575940628739560
  149. VL2a_pn = 720575940618295454
  150. VL2p_pn = 720575940628467611
  151. VM4_pn = 720575940640971995
  152. DM1_pn = 720575940630770042
  153. V_bilateral_rightsoma_pn=720575940626143806
  154. V_bilateral_leftsoma_pn=720575940630546540
  155. V_unilateral_pn=720575940637910106
  156. V_pn = [720575940626143806,720575940630546540,720575940637910106]
  157. D_pn = [720575940603985952, 720575940633804212, 720575940627844734]
  158. DA1_pn = [720575940614309535, 720575940637208718, 720575940613345442, 720575940626034819, 720575940619385765, 720575940605102694, 720575940603231916, 720575940621185050,720575940646122804]
  159. DA2_pn = [720575940622364184, 720575940611849187, 720575940626937617, 720575940638719104,720575940622734835]
  160. DA4m_pn = 720575940631930828
  161. DC2_pn = [720575940627160322, 720575940616824588]
  162. DC3_pn = [720575940630131663, 720575940613419695, 720575940634939327]
  163. DL1_pn = [720575940616079035, 720575940622368792]
  164. DL2d_pn = [720575940628316047, 720575940630928715, 720575940620103590, 720575940642046453,720575940605421442, 720575940613706866, 720575940623922869]
  165. DL3_pn = [720575940630778428, 720575940621926669, 720575940627353426, 720575940630778428, 720575940620638335]
  166. DM2_pn = [720575940638633535, 720575940630024566]
  167. DM5_pn = [720575940644258327, 720575940629926159, 720575940640032896]
  168. DM6_pn = [720575940631328104, 720575940622785062, 720575940620437339]
  169. DP1l_pn = [720575940648650884, 720575940606762686]
  170. VA1d_pn = [720575940622811430, 720575940640804853, 720575940619620852]
  171. VA1v_pn = [720575940629097922, 720575940620199962, 720575940628283560, 720575940629733626, 720575940632720026, 720575940634229615]
  172. VA3_pn = [720575940656090017, 720575940604431584]
  173. VA5_pn = [720575940619838912, 720575940637009636, 720575940637231593]
  174. VA7m_pn = [720575940620068518, 720575940640398477, 720575940639680957]
  175. VC5_pn = [720575940608845525, 720575940608199748, 720575940624515571] #VC3m PNs due to naming changing VC3m -> VC5
  176. VL1_pn = [720575940624863015, 720575940620595053]
  177. VM1_pn = [720575940628077568, 720575940623444715]
  178. VM2_pn = [720575940640697808, 720575940617179451]
  179. VM3_pn = [720575940627480586, 720575940610048981]
  180. VM6_pn = 720575940624057853
  181. VM7d_pn = [720575940628524292, 720575940613167007]
  182. VM7v_pn = [720575940619337536, 720575940621256212]
  183. VM5d_pn = [720575940611372714, 720575940644961588, 720575940635082871, 720575940627035388]
  185. # added from 2020 Current Biology
  186. DA3_pn = [720575940625924618, 720575940659400577]
  187. DC4_pn = 720575940613579943
  188. DL2v_pn = [720575940620467438,720575940611022515,720575940628181520 ,720575940627772009]
  189. VC3_pn = [720575940615394719,720575940619902598 ,720575940619928429 ,720575940620465904]
  190. VC4_pn = [720575940624001321,720575940635933119,720575940609454603]
  191. VM5v_pn = [720575940610505170 ,720575940620189790 ,720575940625431866 ]
  193. # %%
  194. #load .ply files
  196. dm1_l=nv.read_mesh("../glom_meshes/DM1_R.ply", output='volume') # DM1
  197. v_l=nv.read_mesh("../glom_meshes/V_R.ply", output='volume') # V
  199. gloms = ["DA4l", "DC1", "DL4", "DL5", "DM3","DM4", "DP1m",
  200. "VA2", "VA4", "VA6", "VA7l", "VC1", "VC2", "VL2a",
  201. "VL2p", "VM4","DM1","V",
  202. "D",
  203. "DA1",
  204. "DA2",
  205. "DA4m",
  206. "DC2",
  207. "DC3",
  208. "DL1",
  209. "DL2d",
  210. "DL3",
  211. "DM2",
  212. "DM5",
  213. "DM6",
  214. "DP1l",
  215. "VA1d",
  216. "VA1v",
  217. "VA3",
  218. "VA5",
  219. "VA7m",
  220. "VL1",
  221. "VM1",
  222. "VM2",
  223. "VM3",
  224. "VM7d",
  225. "VM7v",
  226. "VM5d",
  227. "DA3",
  228. "DC4",
  229. "DL2v",
  230. "VC3",
  231. 'VC5',
  232. 'VM6',
  233. "VC4",
  234. "VM5v",
  235. ]
  237. for glom in gloms:
  238. print(glom)
  239. exec(glom +"_mesh = nv.read_mesh(" + "'../glom_meshes/" + glom + "_R.ply' , output='volume')")
  241. # %%
  242. #retrieve neuron_ids
  243. glomeruli_list = ["ORN_DA4l",
  244. "ORN_DC1",
  245. "ORN_DL4",
  246. "ORN_DL5",
  247. "ORN_DM3",
  248. "ORN_DM4",
  249. "ORN_DP1m",
  250. "ORN_VA2",
  251. "ORN_VA4",
  252. "ORN_VA6",
  253. "ORN_VA7l",
  254. "ORN_VC1",
  255. "ORN_VC2",
  256. "ORN_VL2a",
  257. "ORN_VL2p",
  258. "ORN_VM4",
  259. "ORN_DM1",
  260. "ORN_V",
  261. "ORN_D",
  262. "ORN_DA1",
  263. "ORN_DA2",
  264. "ORN_DA4m",
  265. "ORN_DC2",
  266. "ORN_DC3",
  267. "ORN_DL1",
  268. "ORN_DL2d",
  269. "ORN_DL3",
  270. "ORN_DM2",
  271. "ORN_DM5",
  272. "ORN_DM6",
  273. "ORN_DP1l",
  274. "ORN_VA1d",
  275. "ORN_VA1v",
  276. "ORN_VA3",
  277. "ORN_VA5",
  278. "ORN_VA7m",
  279. "ORN_VC3",
  280. 'ORN_VC5',
  281. 'ORN_VM6',
  282. "ORN_VL1",
  283. "ORN_VM1",
  284. "ORN_VM2",
  285. "ORN_VM3",
  286. "ORN_VM7d",
  287. "ORN_VM7v",
  288. "ORN_VM5d",
  289. "ORN_DA3",
  290. "ORN_DC4",
  291. "ORN_DL2v",
  292. "ORN_VC3l",
  293. "ORN_VC4",
  294. "ORN_VM5v",
  295. ]
  297. glomeruli_OSNs = pd.DataFrame()
  298. for glomeruli in glomeruli_list:
  299. if glomeruli in ['ORN_VC3','ORN_VC5','ORN_VM6']:
  300. annotation_search = flywire.search_annotations(flywire.NeuronCriteria(cell_type=glomeruli, side='left', annotation_version="v2.0.0"))
  301. else:
  302. print(glomeruli)
  303. annotation_search = flywire.search_annotations(flywire.NeuronCriteria(hemibrain_type=glomeruli, side='left', annotation_version="v2.0.0"))
  304. print(f' for {glomeruli} {annotation_search.shape[0]} neurons were found')
  305. glomeruli_OSNs = pd.concat([glomeruli_OSNs, annotation_search])
  306. glomeruli_OSNs = glomeruli_OSNs[['root_id','hemibrain_type','cell_type']] #keep just the root_id and the glomeruli
  308. # %%
  309. #the following neurons have no synapses with projection neurons in their glomerulus, so they were excluded from all analysis
  310. neurons_to_exclude = [720575940642936520, 720575940623359466, 720575940633033133, 720575940621862605,
  311. 720575940605724081, 720575940604876977, 720575940630971434, 720575940632824531,
  312. 720575940611791834, 720575940614890642, 720575940611437973, 720575940619969329,
  313. 720575940614019112, 720575940615639199, 720575940639118653, 720575940615629173,
  314. 720575940621082138, 720575940608615363, 720575940630679749, 720575940628398791,
  315. 720575940639089571, 720575940634883179, 720575940638202852, 720575940614006934,
  316. 720575940603870688, 720575940632033875, 720575940632391736, 720575940623750852,
  317. 720575940622248709, 720575940631930415, 720575940613885506, 720575940627787343,
  318. 720575940624671490, 720575940631393823, 720575940619029589, 720575940623655605,
  319. 720575940630462572,
  320. ]
  322. # %%
  323. glomeruli_OSNs = glomeruli_OSNs[~glomeruli_OSNs['root_id'].isin(neurons_to_exclude)]
  325. # %%
  326. glom_list = ["D",
  327. "DA1",
  328. "DA2",
  329. "DA3",
  330. "DA4l",
  331. "DA4m",
  332. "DC1",
  333. "DC2",
  334. "DC3",
  335. "DC4",
  336. "DL1",
  337. "DL2d",
  338. "DL2v",
  339. "DL3",
  340. "DL4",
  341. "DL5",
  342. "DM1",
  343. "DM2",
  344. "DM3",
  345. "DM4",
  346. "DM5",
  347. "DM6",
  348. "DP1l",
  349. "DP1m",
  350. "V",
  351. "VA1d",
  352. "VA1v",
  353. "VA2",
  354. "VA3",
  355. "VA4",
  356. "VA5",
  357. "VA6",
  358. "VA7l",
  359. "VA7m",
  360. "VC1",
  361. "VC2",
  362. "VC3",
  363. "VC4",
  364. "VC5",
  365. "VL1",
  366. "VL2a",
  367. "VL2p",
  368. "VM1",
  369. "VM2",
  370. "VM3",
  371. "VM4",
  372. "VM5d"
  373. "VM5v",
  374. "VM6",
  375. "VM7d",
  376. "VM7v"]
  380. for glom in glom_list:
  381. retries = 5
  382. for retry in range(retries):
  383. try:
  384. print(glom)
  385. if glom in ['VM6']:
  386. exec(f"{glom}_feedforward_synapses = flywire.get_adjacency(glomeruli_OSNs[(glomeruli_OSNs['cell_type']=='ORN_{glom}v') | (glomeruli_OSNs['cell_type']=='ORN_{glom}l') | (glomeruli_OSNs['cell_type']=='ORN_{glom}m')]['root_id'], {glom}_pn, min_score=65, dataset='production')")#.sum(axis=1)")
  387. exec(f"{glom}_recurrent_synapses = flywire.get_adjacency(glomeruli_OSNs[(glomeruli_OSNs['cell_type']=='ORN_{glom}v') | (glomeruli_OSNs['cell_type']=='ORN_{glom}l') | (glomeruli_OSNs['cell_type']=='ORN_{glom}m')]['root_id'], min_score=65, dataset='production')")#.sum(axis=1)")
  388. exec(f"{glom}_total_synapses = get_total_synapses(glomeruli_OSNs[(glomeruli_OSNs['cell_type']=='ORN_{glom}v') | (glomeruli_OSNs['cell_type']=='ORN_{glom}l') | (glomeruli_OSNs['cell_type']=='ORN_{glom}m')]['root_id'], volume = {glom}_mesh)")
  389. exec(f"{glom}_feedforward = get_cable_overlap(glomeruli_OSNs[(glomeruli_OSNs['cell_type']=='ORN_{glom}v') | (glomeruli_OSNs['cell_type']=='ORN_{glom}l') | (glomeruli_OSNs['cell_type']=='ORN_{glom}m')]['root_id'],PN={glom}_pn, volume={glom}_mesh)")
  390. exec(f"{glom}_recurrent = get_cable_overlap(glomeruli_OSNs[(glomeruli_OSNs['cell_type']=='ORN_{glom}v') | (glomeruli_OSNs['cell_type']=='ORN_{glom}l') | (glomeruli_OSNs['cell_type']=='ORN_{glom}m')]['root_id'],volume={glom}_mesh)")
  391. exec(f"{glom}_cable_length = get_cable_length(glomeruli_OSNs[(glomeruli_OSNs['cell_type']=='ORN_{glom}v') | (glomeruli_OSNs['cell_type']=='ORN_{glom}l') | (glomeruli_OSNs['cell_type']=='ORN_{glom}m')]['root_id'], volume={glom}_mesh)")
  392. elif glom in ['VC3','VC5']:
  393. exec(f"{glom}_feedforward_synapses = flywire.get_adjacency(glomeruli_OSNs[glomeruli_OSNs['cell_type']=='ORN_{glom}']['root_id'], {glom}_pn, min_score=65, dataset='production')")#.sum(axis=1)")
  394. exec(f"{glom}_recurrent_synapses = flywire.get_adjacency(glomeruli_OSNs[glomeruli_OSNs['cell_type']=='ORN_{glom}']['root_id'], min_score=65, dataset='production')")#.sum(axis=1)")
  395. exec(f"{glom}_total_synapses = get_total_synapses(glomeruli_OSNs[glomeruli_OSNs['cell_type']=='ORN_{glom}']['root_id'], volume = {glom}_mesh)")
  396. exec(f"{glom}_feedforward = get_cable_overlap(glomeruli_OSNs[glomeruli_OSNs['cell_type']=='ORN_{glom}']['root_id'],PN={glom}_pn, volume={glom}_mesh)")
  397. exec(f"{glom}_recurrent = get_cable_overlap(glomeruli_OSNs[glomeruli_OSNs['cell_type']=='ORN_{glom}']['root_id'],volume={glom}_mesh)")
  398. exec(f"{glom}_cable_length = get_cable_length(glomeruli_OSNs[glomeruli_OSNs['cell_type']=='ORN_{glom}']['root_id'], volume={glom}_mesh)")
  399. else:
  400. exec(f"{glom}_feedforward_synapses = flywire.get_adjacency(glomeruli_OSNs[glomeruli_OSNs['hemibrain_type']=='ORN_{glom}']['root_id'], {glom}_pn, min_score=65, dataset='production')")#.sum(axis=1)")
  401. exec(f"{glom}_recurrent_synapses = flywire.get_adjacency(glomeruli_OSNs[glomeruli_OSNs['hemibrain_type']=='ORN_{glom}']['root_id'], min_score=65, dataset='production')")#.sum(axis=1)")
  402. exec(f"{glom}_total_synapses = get_total_synapses(glomeruli_OSNs[glomeruli_OSNs['hemibrain_type']=='ORN_{glom}']['root_id'], volume = {glom}_mesh)")
  403. exec(f"{glom}_feedforward = get_cable_overlap(glomeruli_OSNs[glomeruli_OSNs['hemibrain_type']=='ORN_{glom}']['root_id'],PN={glom}_pn, volume={glom}_mesh)")
  404. exec(f"{glom}_recurrent = get_cable_overlap(glomeruli_OSNs[glomeruli_OSNs['hemibrain_type']=='ORN_{glom}']['root_id'],volume={glom}_mesh)")
  405. exec(f"{glom}_cable_length = get_cable_length(glomeruli_OSNs[glomeruli_OSNs['hemibrain_type']=='ORN_{glom}']['root_id'], volume={glom}_mesh)")
  406. exec(f"{glom}_feedforward_synapses_per_micron = synapses_per_micron({glom}_feedforward,PN={glom}_pn)")
  407. exec(f"{glom}_recurrent_synapses_per_micron = synapses_per_micron({glom}_recurrent)")
  408. except:
  409. if retry < retries - 1:
  410. print(f'number of retries left = {retries-retry}')
  411. continue
  412. else:
  413. raise
  414. break
  416. # %% [markdown]
  417. # retrieve mosquito data from catmaid
  419. # %%
  420. #MD1
  421. MD1OSNs = pymaid.get_skids_by_annotation(['innervates MD1', 'left palp nerve', 'sensory neuron', 'PSPs done'], allow_partial = False, intersect = True)
  422. MD1OSNsALL = pymaid.get_skids_by_annotation(['innervates MD1', 'left palp nerve', 'sensory neuron'], allow_partial = False, intersect = True)
  423. MD1volume = pymaid.get_volume('MD1 04/06/21')
  424. MD1neurons = pymaid.get_neuron(MD1OSNs)
  425. MD1neurons = pymaid.CatmaidNeuron.prune_by_volume(MD1neurons,v=MD1volume,mode='IN',prevent_fragments=True)
  426. MD1neuronsALL = pymaid.get_neuron(MD1OSNsALL)
  428. # %%
  430. #MD2
  431. MD2OSNs = pymaid.get_skids_by_annotation(['innervates MD2', 'left palp nerve', 'sensory neuron', 'PSPs done'], allow_partial = False, intersect = True)
  432. MD2OSNsALL = pymaid.get_skids_by_annotation(['innervates MD2', 'left palp nerve', 'sensory neuron'], allow_partial = False, intersect = True)
  433. MD2volume = pymaid.get_volume('MD2 2022')
  434. MD2neurons = pymaid.get_neuron(MD2OSNs)
  435. MD2neurons = pymaid.CatmaidNeuron.prune_by_volume(MD2neurons,v=MD2volume,mode='IN',prevent_fragments=True)
  436. MD2neuronsALL = pymaid.get_neuron(MD2OSNsALL)
  440. # %%
  441. #MD3
  442. MD3OSNs = pymaid.get_skids_by_annotation(['innervates MD3', 'left palp nerve', 'sensory neuron', 'PSPs done'], allow_partial = False, intersect = True)
  443. MD3OSNsALL = pymaid.get_skids_by_annotation(['innervates MD3', 'left palp nerve', 'sensory neuron'], allow_partial = False, intersect = True)
  444. MD3volume = pymaid.get_volume('MD3 2022')
  445. MD3neurons = pymaid.get_neuron(MD3OSNs)
  446. MD3neurons = pymaid.CatmaidNeuron.prune_by_volume(MD3neurons,v=MD3volume,mode='IN',prevent_fragments=True)
  447. MD3neuronsALL = pymaid.get_neuron(MD3OSNsALL)
  451. # %%
  452. #initialize empty dataframe
  453. flyneurons=pd.DataFrame(columns=['id', 'total cable length', 'cable length in glom', 'unilateral feedforward', 'bilateral left feedforward', 'bilateral right feedforward', 'recurrent connections', 'glomerulus'])
  455. md1_synapse_count = []
  456. md1_synapse_count_in_volume = []
  457. md1_total_synapses = []
  458. md1_total_synapses_restricted = []
  459. md1_synapse_count_feedforward = []
  460. md1_synapse_count_feedforward_unrestricted = []
  461. md1_cable_overlap = []
  462. md1_feedforward_density = []
  463. md1_feedforward_density_unrestricted = []
  464. md1_recurrent_density = []
  465. md1_recurrent_density_unrestricted = []
  466. md1_cable_length = []
  467. md1_cable_overlap_percentage = []
  468. md1_feedforward_fraction = []
  469. md1_recurrent_fraction = []
  470. md1_recurrent_fraction_restricted = []
  471. md1_recurrent_synapses = []
  472. md1_recurrent_synapses_restricted = []
  474. for n in MD1neurons:
  475. skeletonid=n.skeleton_id
  476. skelid=int(skeletonid)
  477. print(skelid)
  478. labels = pymaid.get_label_list()
  479. neuronlabels=labels[labels.skeleton_id==skelid]
  480. branchpoint=neuronlabels[neuronlabels.tag=='first branch point']
  482. if len(branchpoint) !=1:
  483. print('error, neuron skelid=%i does not have exactly 1 branchpoint tag'% skelid)
  484. else: bpnode=branchpoint.node_id.values[0]
  486. PN = pymaid.get_neuron('295')
  488. cable_overlap_length = nv.cable_overlap(n, PN, dist= "2 microns")
  490. connections = pymaid.get_partners(n).set_index('skeleton_id')
  491. connections = connections[connections['relation'] == 'downstream']
  492. total_synapses = pd.DataFrame(connections['total']).transpose()
  493. total_synapses.columns.names = ['targets']
  494. total_synapses.index = [n.id]
  495. total_synapses.index.names = ['sources']
  496. total_synapses_restricted = pymaid.filter_connectivity(total_synapses, restrict_to='MD1 04/06/21')
  498. #feedforward connections
  499. feedforward = pymaid.adjacency_matrix(n, targets=295)
  500. glom_connectivity=pymaid.filter_connectivity(feedforward, restrict_to='MD1 04/06/21')
  501. feedforward_value=int(glom_connectivity.iloc[0])
  503. #recurrent connections
  504. recurrent = pymaid.adjacency_matrix(n, targets=MD1neuronsALL)
  505. recurrent_glom_connectivity=pymaid.filter_connectivity(recurrent, restrict_to='MD1 04/06/21')
  506. print(recurrent_glom_connectivity)
  508. cable_overlap = nv.cable_overlap(n,MD1neuronsALL,dist= "2 microns")
  509. synapse_density_unrestricted = recurrent / (cable_overlap.values / 1000)
  510. synapse_density = recurrent_glom_connectivity / (cable_overlap.values / 1000)
  512. md1_total_synapses.append(sum(total_synapses.iloc[0]))
  513. md1_total_synapses_restricted.append(sum(total_synapses_restricted.iloc[0]))
  514. md1_synapse_count.append(recurrent.values)
  515. md1_synapse_count_feedforward.append(feedforward_value)
  516. md1_synapse_count_feedforward_unrestricted.append(int(feedforward.iloc[0]))
  517. md1_cable_overlap.append(cable_overlap.values)
  518. md1_feedforward_density.append((glom_connectivity/(cable_overlap_length/1000)).iloc[0].values[0])
  519. md1_feedforward_density_unrestricted.append((int(feedforward.iloc[0])/(cable_overlap_length/1000)).iloc[0].values[0])
  520. md1_recurrent_density.append(synapse_density.mean(axis=1).iloc[0])
  521. md1_recurrent_density_unrestricted.append(synapse_density_unrestricted.mean(axis=1).iloc[0])
  522. md1_cable_length.append(n.cable_length/1000)
  523. md1_cable_overlap_percentage.append(((cable_overlap.values /1000)/(n.cable_length/1000)).mean())
  524. sum_recurrent=recurrent.iloc[0].sum()
  525. md1_feedforward_fraction.append(feedforward_value / total_synapses_restricted.sum(axis=1)[0])
  526. md1_recurrent_synapses.append(recurrent.iloc[0].sum())
  527. md1_recurrent_synapses_restricted.append(recurrent_glom_connectivity.sum(axis=1)[0])
  528. md1_recurrent_fraction.append(sum_recurrent / total_synapses)
  529. md1_recurrent_fraction_restricted.append(recurrent_glom_connectivity.sum(axis=1)[0]/sum(total_synapses_restricted.iloc[0]))
  532. #append to dataframe
  533. flyneurons= pd.concat([flyneurons,
  534. pd.DataFrame({'id':n,
  535. 'total cable length':n.cable_length,
  536. 'cable length in glom': cable_overlap_length.values[0][0], #dlength,
  537. 'unilateral feedforward':feedforward_value,
  538. 'recurrent connections': sum_recurrent,
  539. 'glomerulus':'Glomerulus 1'}, index = [0])
  540. ],
  541. ignore_index=True)
  543. # %%
  544. #initialize empty dataframe
  545. flyneurons=pd.DataFrame(columns=['id', 'total cable length', 'cable length in glom', 'unilateral feedforward', 'bilateral left feedforward', 'bilateral right feedforward', 'recurrent connections', 'glomerulus'])
  547. md2_synapse_count = []
  548. md2_synapse_count_in_volume = []
  549. md2_total_synapses = []
  550. md2_total_synapses_restricted = []
  551. md2_synapse_count_feedforward = []
  552. md2_synapse_count_feedforward_unrestricted = []
  553. md2_cable_overlap = []
  554. md2_feedforward_density = []
  555. md2_feedforward_density_unrestricted = []
  556. md2_recurrent_density = []
  557. md2_recurrent_density_unrestricted = []
  558. md2_cable_length = []
  559. md2_cable_overlap_percentage = []
  560. md2_feedforward_fraction = []
  561. md2_recurrent_fraction = []
  562. md2_recurrent_fraction_restricted = []
  563. md2_recurrent_synapses = []
  564. md2_recurrent_synapses_restricted = []
  566. for n in MD2neurons:
  567. skeletonid=n.skeleton_id
  568. skelid=int(skeletonid)
  569. print(skelid)
  570. labels = pymaid.get_label_list()
  571. neuronlabels=labels[labels.skeleton_id==skelid]
  572. branchpoint=neuronlabels[neuronlabels.tag=='first branch point']
  574. if len(branchpoint) !=1:
  575. print('error, neuron skelid=%i does not have exactly 1 branchpoint tag'% skelid)
  576. else: bpnode=branchpoint.node_id.values[0]
  578. PN = pymaid.get_neuron('690')
  580. cable_overlap_length = nv.cable_overlap(n, PN, dist= "2 microns")
  582. connections = pymaid.get_partners(n).set_index('skeleton_id')
  583. connections = connections[connections['relation'] == 'downstream']
  584. total_synapses = pd.DataFrame(connections['total']).transpose()
  585. total_synapses.columns.names = ['targets']
  586. total_synapses.index = [n.id]
  587. total_synapses.index.names = ['sources']
  588. total_synapses_restricted = pymaid.filter_connectivity(total_synapses, restrict_to='MD2 2022')
  590. #feedforward connections
  591. feedforward = pymaid.adjacency_matrix(n, targets=690)
  592. glom_connectivity=pymaid.filter_connectivity(feedforward, restrict_to='MD2 2022')
  593. feedforward_value=int(glom_connectivity.iloc[0])
  595. #recurrent connections
  596. recurrent = pymaid.adjacency_matrix(n, targets=MD2neuronsALL)
  597. recurrent_glom_connectivity=pymaid.filter_connectivity(recurrent, restrict_to='MD2 2022')
  598. print(recurrent_glom_connectivity)
  600. cable_overlap = nv.cable_overlap(n,MD2neuronsALL,dist= "2 microns")
  601. synapse_density_unrestricted = recurrent / (cable_overlap.values / 1000)
  602. synapse_density = recurrent_glom_connectivity / (cable_overlap.values / 1000)
  604. md2_total_synapses.append(sum(total_synapses.iloc[0]))
  605. md2_total_synapses_restricted.append(sum(total_synapses_restricted.iloc[0]))
  606. md2_synapse_count.append(recurrent.values)
  607. md2_synapse_count_feedforward.append(feedforward_value)
  608. md2_synapse_count_feedforward_unrestricted.append(int(feedforward.iloc[0]))
  609. md2_cable_overlap.append(cable_overlap.values)
  610. md2_feedforward_density.append((glom_connectivity/(cable_overlap_length/1000)).iloc[0].values[0])
  611. md2_feedforward_density_unrestricted.append((int(feedforward.iloc[0])/(cable_overlap_length/1000)).iloc[0].values[0])
  612. md2_recurrent_density.append(synapse_density.mean(axis=1).iloc[0])
  613. md2_recurrent_density_unrestricted.append(synapse_density_unrestricted.mean(axis=1).iloc[0])
  614. md2_cable_length.append(n.cable_length/1000)
  615. md2_cable_overlap_percentage.append(((cable_overlap.values /1000)/(n.cable_length/1000)).mean())
  616. sum_recurrent=recurrent.iloc[0].sum()
  617. md2_feedforward_fraction.append(feedforward_value / total_synapses_restricted.sum(axis=1)[0])
  618. md2_recurrent_synapses.append(recurrent.iloc[0].sum())
  619. md2_recurrent_synapses_restricted.append(recurrent_glom_connectivity.sum(axis=1)[0])
  620. md2_recurrent_fraction.append(sum_recurrent / total_synapses)
  621. md2_recurrent_fraction_restricted.append(recurrent_glom_connectivity.sum(axis=1)[0]/sum(total_synapses_restricted.iloc[0]))
  624. #append to dataframe
  625. flyneurons= pd.concat([flyneurons,
  626. pd.DataFrame({'id':n,
  627. 'total cable length':n.cable_length,
  628. 'cable length in glom': cable_overlap_length.values[0][0], #dlength,
  629. 'unilateral feedforward':feedforward_value,
  630. 'recurrent connections': sum_recurrent,
  631. 'glomerulus':'Glomerulus 1'}, index = [0])
  632. ],
  633. ignore_index=True)
  635. # %%
  636. #initialize empty dataframe
  637. flyneurons=pd.DataFrame(columns=['id', 'total cable length', 'cable length in glom', 'unilateral feedforward', 'bilateral left feedforward', 'bilateral right feedforward', 'recurrent connections', 'glomerulus'])
  639. md3_synapse_count = []
  640. md3_synapse_count_in_volume = []
  641. md3_total_synapses = []
  642. md3_total_synapses_restricted = []
  643. md3_synapse_count_feedforward = []
  644. md3_synapse_count_feedforward_unrestricted = []
  645. md3_cable_overlap = []
  646. md3_feedforward_density = []
  647. md3_feedforward_density_unrestricted = []
  648. md3_recurrent_density = []
  649. md3_recurrent_density_unrestricted = []
  650. md3_cable_length = []
  651. md3_cable_overlap_percentage = []
  652. md3_feedforward_fraction = []
  653. md3_recurrent_fraction = []
  654. md3_recurrent_fraction_restricted = []
  655. md3_recurrent_synapses = []
  656. md3_recurrent_synapses_restricted = []
  658. for n in MD3neurons:
  659. skeletonid=n.skeleton_id
  660. skelid=int(skeletonid)
  661. print(skelid)
  662. labels = pymaid.get_label_list()
  663. neuronlabels=labels[labels.skeleton_id==skelid]
  664. branchpoint=neuronlabels[neuronlabels.tag=='first branch point']
  666. if len(branchpoint) !=1:
  667. print('error, neuron skelid=%i does not have exactly 1 branchpoint tag'% skelid)
  668. else: bpnode=branchpoint.node_id.values[0]
  670. PN = pymaid.get_neuron('11126')
  672. cable_overlap_length = nv.cable_overlap(n, PN, dist= "2 microns")
  674. connections = pymaid.get_partners(n).set_index('skeleton_id')
  675. connections = connections[connections['relation'] == 'downstream']
  676. total_synapses = pd.DataFrame(connections['total']).transpose()
  677. total_synapses.columns.names = ['targets']
  678. total_synapses.index = [n.id]
  679. total_synapses.index.names = ['sources']
  680. total_synapses_restricted = pymaid.filter_connectivity(total_synapses, restrict_to='MD3 2022')
  682. #feedforward connections
  683. feedforward = pymaid.adjacency_matrix(n, targets=11126)
  684. glom_connectivity=pymaid.filter_connectivity(feedforward, restrict_to='MD3 2022')
  685. feedforward_value=int(glom_connectivity.iloc[0])
  687. #recurrent connections
  688. recurrent = pymaid.adjacency_matrix(n, targets=MD3neuronsALL)
  689. recurrent_glom_connectivity=pymaid.filter_connectivity(recurrent, restrict_to='MD3 2022')
  690. print(recurrent_glom_connectivity)
  692. cable_overlap = nv.cable_overlap(n,MD3neuronsALL,dist= "2 microns")
  693. synapse_density_unrestricted = recurrent / (cable_overlap.values / 1000)
  694. synapse_density = recurrent_glom_connectivity / (cable_overlap.values / 1000)
  696. md3_total_synapses.append(sum(total_synapses.iloc[0]))
  697. md3_total_synapses_restricted.append(sum(total_synapses_restricted.iloc[0]))
  698. md3_synapse_count.append(recurrent.values)
  699. md3_synapse_count_feedforward.append(feedforward_value)
  700. md3_synapse_count_feedforward_unrestricted.append(int(feedforward.iloc[0]))
  701. md3_cable_overlap.append(cable_overlap.values)
  702. md3_feedforward_density.append((glom_connectivity/(cable_overlap_length/1000)).iloc[0].values[0])
  703. md3_feedforward_density_unrestricted.append((int(feedforward.iloc[0])/(cable_overlap_length/1000)).iloc[0].values[0])
  704. md3_recurrent_density.append(synapse_density.mean(axis=1).iloc[0])
  705. md3_recurrent_density_unrestricted.append(synapse_density_unrestricted.mean(axis=1).iloc[0])
  706. md3_cable_length.append(n.cable_length/1000)
  707. md3_cable_overlap_percentage.append(((cable_overlap.values /1000)/(n.cable_length/1000)).mean())
  708. sum_recurrent=recurrent.iloc[0].sum()
  709. md3_feedforward_fraction.append(feedforward_value / total_synapses_restricted.sum(axis=1)[0])
  710. md3_recurrent_synapses.append(recurrent.iloc[0].sum())
  711. md3_recurrent_synapses_restricted.append(recurrent_glom_connectivity.sum(axis=1)[0])
  712. md3_recurrent_fraction.append(sum_recurrent / total_synapses)
  713. md3_recurrent_fraction_restricted.append(recurrent_glom_connectivity.sum(axis=1)[0]/sum(total_synapses_restricted.iloc[0]))
  716. #append to dataframe
  717. flyneurons= pd.concat([flyneurons,
  718. pd.DataFrame({'id':n,
  719. 'total cable length':n.cable_length,
  720. 'cable length in glom': cable_overlap_length.values[0][0], #dlength,
  721. 'unilateral feedforward':feedforward_value,
  722. 'recurrent connections': sum_recurrent,
  723. 'glomerulus':'Glomerulus 1'}, index = [0])
  724. ],
  725. ignore_index=True)

retrieve_neuron_data.ipynb at commit 81b5589, under GPL-3.0 · at the source

Overview

Authors: Jialu Bao1,2, Wesley Alford3, Avinash Khandelwal1, Laurel Walsh1, George Lantz1, Santiago Poncio1, Laia Serratosa Capdevila1, Yervand Azatian1, Brian DePasquale4, David G C Hildebrand5,6, Meg A Younger3, Wei-Chung Allen Lee1,7
ORCID iDs: Wei-Chung Allen Lee
  1. Department of Neurobiology, Harvard Medical School, Boston, Massachusetts, United States of America
  2. Molecules, Cells, and Organisms Training Program, Department of Molecular and Cellular Biology, Harvard University, Cambridge, Massachusetts, United States of America
  3. Department of Biology, Boston University, Boston, Massachusetts, United States of America
  4. Department of Biomedical Engineering, Boston University, Boston, Massachusetts, United States of America
  5. Laboratory of Neural Systems, The Rockefeller University, New York, New York, United States of America
  6. Department of Vision Sciences, College of Optometry, University of Houston, Houston, Texas, United States of America
  7. F.M. Kirby Neurobiology Center, Boston Children’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
Institutions: Harvard University (United States); Boston University (United States); University of Houston (United States); Rockefeller University (United States); Boston Children's Hospital (United States)
Journal: PLoS biology, volume 24, issue 9, article e3003959
Dates: received 19 January 2026; accepted 5 August 2026; published online 10 September 2026
Type: Research article · Language: English
License: CC BY
Identifiers: DOI 10.1371/journal.pbio.3003959 · PMID 42721124 · PMCID PMC13561330 · OpenAlex W7212152300
Open access: gold, a free copy (OpenAlex)
Status: code verified
Categories: histology / microscopy (modality), drosophila (organism), cellular / molecular (subfield)
Methods: Statistics, fMRI & imaging
MeSH: Aedes*, Carbon Dioxide*, Olfactory Receptor Neurons*, Synapses*, Animals, Drosophila melanogaster, Female, Microscopy, Electron, Transmission, Odorants, Smell (* major topic)
Topic: Neurobiology and Insect Physiology Research (Cellular and Molecular Neuroscience, Neuroscience), according to OpenAlex
Funding: Alfred P. Sloan Foundation (n/a); National Institute of Neurological Disorders and Stroke (R01NS121874); NCRR NIH HHS (S10 RR028832); National Science Foundation Graduate Research Fellowship Program (n/a); Pew Charitable Trusts (n/a); Searle Scholars Program (n/a); NINDS NIH HHS (R01 NS121874); Wellcome Trust (316808/Z/24/Z); NHLBI NIH HHS (DP2 HL185099); National Heart, Lung, and Blood Institute of the National Institutes of Health (DP2HL185099); HMS Genise Goldenson Award (n/a); Leon Levy Foundation (n/a); Simons Foundation (n/a); Esther A. and Joseph Klingenstein Fund (n/a); Richard and Susan Smith Family Foundation (n/a)
Citations: cited by 1 paper (Europe PMC); 115 references in the paper
Notices: A comment on this paper has been published (42726699, from Europe PMC)

Abstract

The mosquito Aedes aegypti’s human host-seeking behavior depends on the integration of multiple sensory cues. One of these cues, carbon dioxide (CO2), gates odorant and heat pathways and activates host-seeking behavior. The neuronal circuits underlying processing of CO2 information remain unclear. We used automated serial-section transmission electron microscopy (EM) to image and reconstruct the circuitry of the glomeruli that are innervated by the Ae. aegypti maxillary palp, including the glomerulus that responds to CO2. Notably, CO2-sensitive olfactory sensory neurons (OSNs) make high levels of recurrent synaptic connections with one another, while making a low density of feedforward synapses. At some of these contacts between CO2 OSNs, we observe ribbon-like presynaptic structures, which may further enhance recurrent signaling. We compared both feedforward and recurrent connectivity with all olfactory glomeruli in Drosophila melanogaster, and we found more recurrent connections between the Ae. aegypti CO2-responsive OSNs than in any D. melanogaster glomeruli. We developed a computational circuit model that demonstrates recurrent synapses are necessary for robust CO2 detection under normal physiological conditions. Together, elevated levels of recurrent connectivity and ribbon-like structures may amplify sensory information detected by CO2-sensitive OSNs to support mosquito activation and sensitization by CO2, even in the presence of high levels of other odorants in the environment. We propose that this circuit organization supports the salience of CO2 as a mosquito host cue.

Reproduced under the paper's license (CC BY), from the paper cited above.

Repositories

Its files are read in the Code ↔ Paper reader above, with 4 matches between paragraphs and lines of code.

htem/aedes_public

License: GPL-3.0
State: the link answers, verified on 26 September 2026
Evidence: files inventoried
Commit: 81b5589a63e067c7e35601ee7ebceb342a082d7b, 16 July 2026
Languages: Jupyter (68), Python (1)
Size: 235 files, 69 scripts
Software Heritage: not archived
Found in: “Data Availability”
Holds: README, license file, environment (model/requirements.txt), 53 notebooks
Not found: CITATION.cff, tests, continuous integration, documentation
Tools: Matplotlib (67 files), seaborn (65 files), NumPy (58 files), pandas (57 files), SciPy (56 files), scikit-posthocs (51 files), statsmodels (37 files), scikit-learn (5 files), Scanpy (2 files)
Availability: 1 check, the latest on 26 September 2026: the link answers
  • 26 September 2026: the link answers
71 files
At the source: github.com/htem/aedes_public

Zenodo 21403889

License: GPL-3.0
State: the link answers, verified on 26 September 2026
Evidence: files inventoried
Size: 1 file
Software Heritage: not checked
Found in: “Data Availability”
Not found: README, license file, CITATION.cff, environment file, tests, continuous integration, documentation
Tools: Matplotlib (67 files), seaborn (65 files), NumPy (58 files), pandas (57 files), SciPy (56 files), scikit-posthocs (51 files), statsmodels (37 files), scikit-learn (5 files), Scanpy (2 files)
Availability: 1 check, the latest on 26 September 2026: the link answers (HTTP 200)
  • 26 September 2026: the link answers (HTTP 200)
71 files
At the source: doi.org/10.5281/zenodo.21403889

The paper's code and data availability statement is in the Data section.

Tracing map

Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.

What the map holds:

  • 2 repositories of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
  • 138 scripts, each with its path and the digest of its content;
  • 4 matches between paragraphs of the paper and lines of the code (method lexical-v1);
  • neither the text of the paper nor the code itself.

Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.

Data

No dataset and no data link were found in the paper.

Data Availability

Data are publicly available at: doi: https://doi.org/10.60533/boss-2025-gcyr. Data and code are also available at: https://github.com/htem/aedes_public and https://doi.org/10.5281/zenodo.21403889.

Reproduced under the paper's license (CC BY), from the paper cited above.

Versions

The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.

Version 1, 27 September 2026: the first record

Recorded: type, language, journal, volume, issue, pages, dates, 12 authors, 10 MeSH terms, 15 funders, 113 references, 1 integrity notice.

Cite

This paper

Bao, J., Alford, W., Khandelwal, A., Walsh, L., Lantz, G., Poncio, S., Capdevila, L. S., Azatian, Y., DePasquale, B., Hildebrand, D. G. C., Younger, M. A., & Lee, W.-C. A. (2026). Recurrent synapses between CO2-sensitive olfactory sensory neurons enable robust CO2 detection in Aedes aegypti mosquitoes. PLoS biology, 24(9), e3003959. https://doi.org/10.1371/journal.pbio.3003959

BibTeX

@article{bao2026recurrent,
author = {Bao, Jialu and Alford, Wesley and Khandelwal, Avinash and Walsh, Laurel and Lantz, George and Poncio, Santiago and Capdevila, Laia Serratosa and Azatian, Yervand and DePasquale, Brian and Hildebrand, David G C and Younger, Meg A and Lee, Wei-Chung Allen},
title = {{Recurrent synapses between CO2-sensitive olfactory sensory neurons enable robust CO2 detection in Aedes aegypti mosquitoes}},
journal = {PLoS biology},
year = {2026},
month = sep,
volume = {24},
number = {9},
pages = {e3003959},
publisher = {PLOS},
issn = {1544-9173},
doi = {10.1371/journal.pbio.3003959},
url = {https://doi.org/10.1371/journal.pbio.3003959},
pmid = {42721124},
pmcid = {PMC13561330}
}

RIS

TY - JOUR
AU - Bao, Jialu
AU - Alford, Wesley
AU - Khandelwal, Avinash
AU - Walsh, Laurel
AU - Lantz, George
AU - Poncio, Santiago
AU - Capdevila, Laia Serratosa
AU - Azatian, Yervand
AU - DePasquale, Brian
AU - Hildebrand, David G C
AU - Younger, Meg A
AU - Lee, Wei-Chung Allen
TI - Recurrent synapses between CO2-sensitive olfactory sensory neurons enable robust CO2 detection in Aedes aegypti mosquitoes
T2 - PLoS biology
J2 - PLoS Biol
PY - 2026
DA - 2026/09/10
VL - 24
IS - 9
SP - e3003959
SN - 1544-9173
PB - PLOS
DO - 10.1371/journal.pbio.3003959
UR - https://doi.org/10.1371/journal.pbio.3003959
LA - en
ER -

CSL-JSON

{
"id": "10.1371/journal.pbio.3003959",
"type": "article-journal",
"title": "Recurrent synapses between CO2-sensitive olfactory sensory neurons enable robust CO2 detection in Aedes aegypti mosquitoes",
"container-title": "PLoS biology",
"author": [
{
"family": "Bao",
"given": "Jialu"
},
{
"family": "Alford",
"given": "Wesley"
},
{
"family": "Khandelwal",
"given": "Avinash"
},
{
"family": "Walsh",
"given": "Laurel"
},
{
"family": "Lantz",
"given": "George"
},
{
"family": "Poncio",
"given": "Santiago"
},
{
"family": "Capdevila",
"given": "Laia Serratosa"
},
{
"family": "Azatian",
"given": "Yervand"
},
{
"family": "DePasquale",
"given": "Brian"
},
{
"family": "Hildebrand",
"given": "David G C"
},
{
"family": "Younger",
"given": "Meg A"
},
{
"family": "Lee",
"given": "Wei-Chung Allen"
}
],
"container-title-short": "PLoS Biol",
"volume": "24",
"issue": "9",
"page": "e3003959",
"DOI": "10.1371/journal.pbio.3003959",
"PMID": "42721124",
"PMCID": "PMC13561330",
"ISSN": "1544-9173",
"publisher": "PLOS",
"URL": "https://doi.org/10.1371/journal.pbio.3003959",
"language": "en",
"issued": {
"date-parts": [
[
2026,
9,
10
]
]
}
}

The tracing map gets a citation of its own once an author has validated it and it has a DOI.

Similar papers

The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.

[1] doi:10.3389/fnsys.2026.1822122 [code]
Convergence-divergence circuits for multimodal integration of innate and learned opponent valences.
Journal: Frontiers in systems neuroscience
In common: seaborn, scikit-learn, pandas, 3 other tools, 23 references
[2] doi:10.1038/s41586-026-10735-w [code]
Distributed control circuits across a brain-and-cord connectome.
Journal: Nature
In common: seaborn, scikit-learn, pandas, 3 other tools, drosophila, 17 references, author Wei-Chung Allen Lee
[3] doi:10.1038/s41467-026-72152-x [code]
Centralized brain networks controlling antennal grooming coordination.
Journal: Nature communications
In common: statsmodels, seaborn, scikit-learn, 4 other tools, drosophila, 8 references
[4] doi:10.1126/sciadv.aeg3535 [code]
Undoing of firing rate adaptation enables invariant population codes.
Journal: Science advances
In common: scikit-learn, SciPy, Matplotlib, 1 other tool, cellular / molecular, 9 references
[5] doi:10.1038/s41467-026-76682-2 [code]
Odorant receptor coexpression and multi-expression in the dengue mosquito.
Journal: Nature communications
In common: NumPy, drosophila, cellular / molecular, 9 references
[6] doi:10.1016/j.isci.2026.117359 [code]
Connectomic mapping of pharyngeal and gut sensory circuits in adult &lt;i&gt;Drosophila&lt;/i&gt;.
Journal: iScience
In common: statsmodels, seaborn, pandas, 3 other tools, drosophila, 7 references
[7] doi:10.1371/journal.pcbi.1014657 [code]
Population morphology implies a common developmental blueprint for Drosophila motion detectors.
Journal: PLoS computational biology
In common: statsmodels, seaborn, scikit-learn, 4 other tools, drosophila, 5 references
[8] doi:10.7554/elife.107276 [code]
Deployment of endocytic machinery to periactive zones of nerve terminals is independent of active zone assembly and evoked release.
Journal: eLife
In common: scikit-posthocs, seaborn, scikit-learn, 4 other tools, drosophila, cellular / molecular, 3 references
[9] doi:10.1371/journal.pbio.3003955 [code]
Beat- and Side-family cell-surface molecules are expressed combinatorially in the partner neurons of the olfactory circuit in Drosophila.
Journal: PLoS biology
In common: drosophila, cellular / molecular, 7 references
[10] doi:10.1038/s41593-026-02388-9 [code]
Hippocampal CA3 connectomics reveals a gradient of mossy fiber inputs and selective feedforward inhibition onto pyramidal cells.
Journal: Nature neuroscience
In common: seaborn, scikit-learn, pandas, 3 other tools, histology / microscopy, 4 references

Contribute

The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.

Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.

Request its removal

To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).

Discussion, reproductions, activity

Discussion: questions and error reports about this paper and its code, from signed-in readers and its authors. It opens with sign-in.

Reproductions: reports from readers who ran the authors' code: what they reproduced, with which environment, commit and data. It opens with sign-in.

Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.