Ao longo de duas iniciativas separadas de conservação voltadas para a recuperação de Cyclargus thomasi bethunebakeri de fevereiro de 2003 a dezembro de 2010 e de novembro de 2016 até o presente, este protocolo foi usado para produzir com sucesso um excesso de 51.052 organismos viáveis. Com base no instantâneo resumo de um ano da produtividade geral da população cativa de junho de 2018 a junho de 2019, foram produzidos 10.166 organismos viáveis, representando 782,00 ± 118,93 organismos por mês em 13 gerações. Da mesma forma, a produção total total de ovos por fêmea em condições laboratoriais foi de 114,00 ± 26,12 (n = 12)31. A produtividade substancial do organismo está neste programa entre os maiores esforços ex situ nos EUA, juntamente com os de Euphydryas editha taylori, Speyeria zerene hippolyta, e Lycaeides melissa samuelis24. Parte dessa produtividade pode ser atribuída ao fato de que a borboleta é continuamente ninhada, produzindo uma geração aproximadamente a cada 4-6 semanas em cativeiro. A maioria dos outros programas de reprodução de conservação envolve taxas univoltinas ou bivoltinas. No entanto, mesmo para programas que envolvem taxas extremamente fecundas, como speyeria spp., o número total de organismos viáveis produzidos para translocação de conservação anualmente raramente excede alguns milharesde 32. Assim, nossa população cativa possibilitou a pesquisa direcionada e a coleta extensiva de dados sobre inúmeras lacunas de dados importantes para melhorar as melhores práticas de reprodução laboratorial e pecuária(Figura 1),bem como ajudar a informar as decisões de recuperação e gestão.
O tempo total de desenvolvimento médio da larva de neonato para adulto foi de 28,63 dias (Tabela 1). A maioria das larvas tinha quatro molts (Figura 2, Figura 3), embora duas tivessem cinco molts, e uma tinha seis molts. O comprimento médio total de todas as estrelas larvais foi de 5,97 mm, e as larvas foram as maiores na quarta fase de vida preparatória(Tabela 1). Quando só inclui variáveis com mais de 30 observações, o menor tempo foi gasto nas primeiras etapas instar e prepupal, e o mais longo foi gasto como pupae (Tabela 1, Figura 2). As fêmeas tipicamente se desenvolveram mais rapidamente em todas as etapas imaturas em comparação com os machos, embora este não tenha sido um efeito significativo (p = 0,625). Foram realizadas análises estatísticas utilizando rStudio versão 1.1.463 (R Core Team 2016)33. O comprimento médio do acorde da asa adulta foi de 12,64 mm(Tabela 2), e houve uma diferença significativa entre os sexos (p = 0,047). O teste t de dois lados foi executado para avaliar a diferença de acordes de asa entre os sexos. Modelo de regressão linear e regressão stepwise para o comprimento médio de cada fase de vida mostraram que o comprimento pupal foi o melhor preditor para o comprimento do acorde da asa adulta (Tabela 3, Tabela 4). Modelos de regressão para o tempo de desenvolvimento mostraram que o número de dias gastos na segunda e quarta estrelas e o número total de dias foram os melhores preditores para o comprimento do acorde da asa adulta, mas apenas o número de dias na quarta instar foi significativo (Tabela 5, Tabela 6). Como as variáveis eram contínuas, dois modelos de regressão linear eram executados para o tempo de desenvolvimento de cada fase de vida, bem como o comprimento de cada fase de vida, com o comprimento do acorde da asa adulta como variável dependente. Regressões stepwise foram executadas em ambos os modelos de regressão para determinar os melhores preditores do comprimento do acorde das asas adultas.

Figura Suplementar 1: Espécimes fixados do adulto Cyclargus thomasi bethunebackeri. (A)Homem adulto, dorsal (esquerda), ventral (direita). (B) Mulher adulta, dorsal (esquerda), ventral (direita). Clique aqui para ver uma versão maior deste valor.

Figura Suplementar 2: Gaiola de voo exibida abrigada em estufa controlada pela temperatura. (A) Interior mostra plantas de néctar adultoem em vasos e uma única planta de hospedeiro larval em vasos. (B)A prateleira metálica ajuda a elevar as plantas de néctar em vasos para que não haja mais de 30 cm de espaço desde o topo interno da gaiola até as flores mais altas florescendo. Clique aqui para ver uma versão maior deste valor.

Figura Suplementar 3: Procedimento para coleta de pares adultos na copula. (A)Par de acasalamento de Cyclargus thomasi bethunebakeri adulto dentro da gaiola de voo screened (feminino, direito e masculino, esquerda). (B) Pares de acasalamento coletados da gaiola de voo em frascos de tampa instantânea e trazidos para o laboratório. Clique aqui para ver uma versão maior deste valor.

Figura Suplementar 4: Procedimento para montagem da câmara de oviposição. (A)Dois sistemas de copos com material hospedeiro terminal e cotonetes de algodão. (B) Uma seringa sub-Q de 1 ml (0,45 mm x 16 mm) com bebidas esportivas aromatizadas saturando cotonetes de algodão no copo de papel. (C) Xícaras de habitação de fêmeas gravid garantidas com tule preto. Clique aqui para ver uma versão maior deste valor.

Figura Suplementar 5: Configuração laboratorial para maximizar a produção de ovos. (A)Câmaras de oviposição colocadas em um banco de laboratório uma luz de grampo com uma lâmpada incandescente de 40 W. (B)Um termômetro de monitoramento de memória rastreável é colocado ao lado das luzes com o sensor de temperatura descansando em cima de uma câmara de oviposição localizada diretamente uma luz de grampo. (C)Uma seringa sub-Q de 1 ml e um pequeno béquer segurando bebida esportiva aromatizada colocada ao lado das câmaras de oviposição para facilitar a refrescação dos cotonetes de algodão regularmente ao longo do dia. Clique aqui para ver uma versão maior deste valor.

Figura Suplementar 6: Configuração laboratorial para cuidados e manutenção larval. (A)Dois sistemas de xícara com cada um contendo material de hospedeiro terminal fresco e larvas. (B) A temperatura nos copos é mantida entre 25 °C-28 °C para atividade larval ideal e desenvolvimento por luzes de fixação aéreas com lâmpadas incandescentes de 40 W. (C)Um termômetro de monitoramento de memória rastreável com o sensor de temperatura colocado diretamente em um copo é usado para monitorar a temperatura. Clique aqui para ver uma versão maior deste valor.

Figura Suplementar 7: Câmaras de pupação preparadas. (A)Copos de porção de plástico individual alojados nas bandejas de copos de plástico transparente. (B)Um quadrado de papel ondulado é colocado em cada xícara de porção de plástico. (C) Uma única larva madura será colocada em cada copo de porção de plástico preparado para pupato. Clique aqui para ver uma versão maior deste valor.

Figura Suplementar 8: Preparando larvas para pupação e manutenção pupal. (A) Larva madura pronta para filhote em papel ondulado. É um uniforme maçante verde-marrom e perdeu qualquer chevrons. (B)Câmaras de pupação prontas para receber larvas maduras adjacentes a copos com larvas de alimentação. Todas as câmaras de pupação com larvas casa de tampas que estão se preparando para pupa. (C) Câmaras de pupação com pupas. (D) Bancos de câmaras de pupação com pupas organizadas por data e mantidos em condições laboratoriais. Clique aqui para ver uma versão maior deste valor.

Figura Suplementar 9: Gaiola de emergência laboratorial. (A)Uma gaiola de criação de malha dobrável que abriga as câmaras de pupação ocupadas. (B)As tampas de todas as câmaras de pupação são removidas para facilitar a eclosão adulta bem sucedida. (C) Todas as borboletas adultas viáveis resultantes serão lançadas na gaiola de voo exibida para garantir uma cópula bem sucedida. Clique aqui para ver uma versão maior deste valor.

Figura Suplementar 10: Borboleta macho adulta se aproximando com sucesso de pupa em um quadrado de papel ondulado. (A)Adulto se aproximando de pupa. (B)Adulto totalmente removido do invólucro pupal. (C) Adulto posicionado para expandir suas asas. (D)Adulto expandindo suas asas. Clique aqui para ver uma versão maior deste valor.

Figura Suplementar 11: Quinta larva instar marcada com tinta luminosa não tóxica. (A)Uma pequena gota de tinta luminosa vermelha e não tóxica contrastante é colocada no dorsum usando um pincel para marcar com sucesso a larva. Clique aqui para ver uma versão maior deste valor.

Figura Complementar 12: Criação de configuração para estudo de história da vida. (A)Exclusivamente rotulados copos de porção de plástico transparente de 2 onças. (B)Uma única larva é seqüerapada em cada xícara. (C) Todas as larvas são rastreadas individualmente em todos os estágios de desenvolvimento, do neonato a borboleta adulta. Clique aqui para ver uma versão maior deste valor.

Figura 1: Número de pares registrados em copula com base na temperatura (°C) dentro de uma gaiola de voo rastreada e exibida abrigada em uma estufa controlada pela temperatura. A temperatura foi registrada nos primeiros 2 minutos de um evento de pareamento bem sucedido (n = 411). Os dados resultantes foram usados para ajudar a refinar as condições ambientais controladas, a fim de maximizar o sucesso do acasalamento e, finalmente, a produtividade geral da propagação em cativeiro. Clique aqui para ver uma versão maior deste valor.

Figura 2: Tempo médio de desenvolvimento (número de dias) de cada fase de vida imatura. (A) As barras mostram a média de cada grupo, e as barras de erro representam os valores de desvio padrão superior e inferior para cada grupo. (B)As barras azuis escuras representam fêmeas, e azul claro representam machos. Clique aqui para ver uma versão maior deste valor.

Figura 3: Cápsulas de cabeça coletadas de #25 individuais usando protocolo de história da vida. Cápsulas de cabeça foram fotografadas por Johnathan Bremer usando um sistema de automontagem. Clique aqui para ver uma versão maior deste valor.
| Fase da vida | Comprimento médio do corpo (mm) | Erro de Std. (comprimento) | Tempo médio de desenvolvimento (num. dias) | Erro de Std. (tempo de v. ) |
| Instar I | 1.69478261 (n=23) | 0.02152643 | 2.90625 (n=32) | 0.08229783 |
| Instar II | 2.77248958 (n=32) | 0.04302826 | 3.375 (n=32) | 0.16649857 |
| Instar III | 5.45751042 (n=32) | 0.12120829 | 3.5 (n=32) | 0.20080483 |
| Instar IV | 10.2369688 (n=32) | 0.23653991 | 3.875 (n=32) | 0.18917265 |
| Instar V | 8,7625 (n=2) | 2.6125 | 1.5 (n=2) | 0.5 |
| Instar VI | 10.2666667 (n=1) | NA | 3 (n=1) | NA |
| Pré-pupa | 11.0858333 (n=24) | 0.23948251 | 2,9375 (n=32) | 0.21504641 |
| Pupa | 9.0316129 (n=31) | 0.12106792 | 11.6578947 (n=38) | 0.3272288 |
Tabela 1: Comprimento médio e tempo de desenvolvimento de cada fase de vida. Erro padrão incluído para cada variável e tamanho da amostra em parênteses.
| Fase da vida | Comprimento médio do acorde da asa (mm) | Erro de DST |
| Adulto | 12.63895 (n=38) | 0.1365516 |
| Fêmea | 12.960 (n=13) | 0.1465588 |
| Masculino | 12.472 (n=25) | 0.1863205 |
Tabela 2: Comprimento médio do acorde das asas para borboletas adultas. Inclui meios para mulheres, machos e todos os adultos (ambos os sexos combinados).
| Modelo LM 1 | Erro de estimativa do Std. | t valor | p-valor |
| Interceptar | 1.9179 | 3.128 | 0.0046 ** |
| Avg. comprimento segunda instar | 0.6822 | -1.11 | 0.278 |
| Avg. comprimento terceira instar | 0.2928 | 0.476 | 0.6381 |
| Avg. comprimento quarta instar | 0.1373 | -0.57 | 0.5739 |
| Pupae de comprimento avg. | 0.246 | 3.957 | 0.0005 *** |
| p < 0,001; ** p < 0,01; * p < 0,05. | | |
Tabela 3: Tabela de coeficientes para modelo de regressão linear (LM Modelo 1) para avaliar relação entre o comprimento médio de cada fase de vida (n > 30 incluídos na análise) e comprimento do acorde da asa adulta. Variável dependente: comprimento do acorde da asa adulta (mm).
| Coeficientes | Erro de estimativa do Std. | t valor | Pr (>|t|) |
| Interceptar | 1.7091 | 3.031 | 0.0053 ** |
| Pupae de comprimento avg. | 0.1878 | 4.414 | 0.0002 *** |
Tabela 4: Regressão stepwise (Stepwise 1). Variável dependente: comprimento do acorde da asa adulta (mm).
| Modelo LM 2 | Erro de estimativa do Std. | t valor | p-valor |
| Interceptar | 1.1888 | 12.643 | 4.21e-12 *** |
| Num. dias primeira instar | 0.3486 | 0.937 | 0.3583 |
| Num. dias segunda instar | 0.2603 | -0.686 | 0.4993 |
| Num. dias terceira instar | 0.2281 | 1.028 | 0.3141 |
| Num. dias quarta instar | 0.2048 | 2.378 | 0.0257 * |
| Num. dias pré-pupae | 0.222 | 1.133 | 0.2686 |
| Num. dias pupae | 0.2495 | 0.616 | 0.5435 |
| Total de num. dias | 0.1913 | -1.454 | 0.1589 |
| p < 0,001; ** p < 0,01; * p < 0,05. | | |
Tabela 5: Tabela coeficientes para modelo de regressão linear (LM Model 2) para avaliar relação entre tempo de desenvolvimento e comprimento do acorde da asa adulta. Variável dependente: comprimento do acorde da asa adulta (mm).
| Coeficientes | Erro de estimativa do Std. | t valor | p-valor |
| Interceptar | 0.89304 | 16.314 | 7.86e-16 *** |
| Num. dias segunda instar | 0.17974 | -1.809 | 0,0811 •
|
| Num. dias quarta instar | 0.16917 | 2.075 | 0.0473 * |
| Total de num. dias | 0.04184 | -1.787 | 0,0848 •
|
| p < 0,001; ** p < 0,01; * p < 0,05; • p < 0.1 | |
Tabela 6: Regressão stepwise (Stepwise 2) para o tempo de desenvolvimento. Variável dependente: comprimento do acorde da asa adulta (mm).