abstract: Gene-centric theories of evolution by natural selection have been popularized and remain generally accepted in both scientific and public paradigms. While gene-centrism is certainly parsimonious, its explanations fall short of describing two patterns of evolutionary and social phenomena: the evolution of sex and the evolution of social altruism. I review and analyze current theories on the evolution of sex. I then introduce the conflict presented to gene-centric evolution by social phenomena such as altruism and caste sterility in eusocial insects. I review gene-centric models of inclusive fitness and kin selection proposed by Hamilton and Maynard Smith. Based their assumptions, that relatedness should be equal between sterile workers and reproductives, I present several empirical examples that conflict with their models. Following that, I introduce a unique system of genetic caste determination (GCD) observed in hybrid populations of two sister-species of seed harvester ants, Pogonomyrmex rugosus and Pogonomyrmex barbatus. I review the evidence for GCD in those species, followed by a critique of the current gene-centric models used to explain it. In chapter two I present my own theoretical model that is both simple and extricable in nature to explain the origin, evolution, and maintenance of GCD in Pogonomyrmex. Furthermore, I use that model to fill in the gaps left behind by the contributing authors of the other GCD models. As both populations in my study system formed from inter-specific hybridization, I review modern discussions of heterosis (also called hybrid vigor) and use those to help explain the ecological competitiveness of GCD. I empirically address the inbreeding depression the lineages of GCD must overcome in order to remain ecologically stable, demonstrating that as a result of their unique system of caste determination, GCD lineages have elevated recombination frequencies. I summarize and conclude with an argument for why GCD evolved under selective mechanisms which cannot be considered gene-centric, providing evidence that natural selection can effectively operate on non-heritable genotypes appearing in groups and other social contexts. / Dissertation/Thesis / M.S. Biology 2012

Vários estudos encontraram maiores taxas de mortalidade por batidas no trânsito entre homens do que entre mulheres e entre homens jovens do que entre homens mais velhos. No entanto, há uma lacuna quanto a explicações para essas diferenças. Esse estudo, diferentemente dos antecessores, parte de hipóteses postuladas a priori, sobre como devem variar as taxas de mortalidade por batidas no trânsito, de acordo com sexo e idade. Espera-se encontrar a presença e ausência de certas associações e interações, com base na teoria moderna da seleção sexual aplicada ao estudo do comportamento humano, especificamente, a síndrome do macho jovem. O objetivo do presente estudo é comparar tendências nas taxas de mortalidade por batida de carro e moto segundo sexo e idade. Trata-se de um estudo ecológico com Estado como unidade geográfica de análise e o ano calendário como unidade temporal de análise. A população do estudo foi composta por homens e mulheres de 18 a 60 anos residentes no Estado do Rio de Janeiro, entre 2004 e 2010, e no Estado do Rio Grande do Sul, entre os anos de 2001 e 2010. Os resultados mostraram que o número de mortes por batida de carro e de moto, considerando o número de habilitados para guiar tais veículos, foi maior em homens do que em mulheres e em indivíduos mais jovens do que mais velhos. Além dessa interessante congruência entre os resultados encontrados e as hipóteses postuladas, descobrimos uma intrigante exceção: no Estado do Rio de Janeiro, a taxa de mortalidade por batida de moto foi consideravelmente maior em mulheres do que em homens. As tendências nas taxas de mortalidade por batida de carro e moto no Estado do Rio de Janeiro vêm apresentando uma queda nos últimos anos. No Rio Grande do Sul, as taxas de mortalidade por batida de moto também vêm caindo ao longo dos anos, enquanto que as taxas de mortalidade por batida de carro vêm apresentando aumento. Tais resultados levam a concluir que a Síndrome do Macho Jovem parece ser uma potente e promissora teoria para ajudar a postular hipóteses sobre o risco de morte no trânsito. / Several studies have found higher rates of mortality in traffic crash among men than among women and among young men than among older men. However, there is a gap concerning explanations for these differences. In this study, unlike predecessors, hypotheses were postulated a priori, on how rates of mortality due to crashes in traffic should vary, according to sex and age. We expected to find the presence and absence of certain associations and interactions, based on the modern theory of sexual selection applied to the study of human behavior, specifically, the "young male syndrome." The aim of this study is to compare trends in mortality rates from car and motorcycle crashes according to sex and age. It is an ecological study with the state as the geographical unit of analysis and the calendar year as the temporal unit of analysis. The study population was composed of 18-60 years old men and women residing in the State of Rio de Janeiro, between 2004 and 2010, and in the state of Rio Grande do Sul, between the years 2001 and 2010. The results have shown that the number of deaths from car and motorcycle crashes considering the number of qualified people to drive such vehicles was higher in men than in women and in younger than older people. In addition to the interesting congruence between the findings and hypotheses postulated, we discovered an intriguing exception: at the state of Rio de Janeiro, the mortality rate from motorcycle crash was considerably higher in women than in men. Trends in mortality rates from car and motorcycle crashes in the state of Rio de Janeiro have shown a decline in recent years. In Rio Grande do Sul, the mortality rates from motorcycle crashes have also been falling over the years, while mortality rates from car crashes have increased. These results lead to conclude that the "Young Male Syndrome" seems to be a powerful and promising theory to help postulating hypotheses about the risk of death in traffic.

Made available in DSpace on 2015-03-26T13:42:28Z (GMT). No. of bitstreams: 1 texto completo.pdf: 2717003 bytes, checksum: 430ed9ce97311b4450b054a78245fcf2 (MD5) Previous issue date: 2013-02-01 / Conselho Nacional de Desenvolvimento Científico e Tecnológico / Cobalamin-independent methionine synthase (MetE) is a key enzyme for de novo synthesis of methionine in plants and it is also responsible for the regeneration of S-adenosylmethionine, which is a metabolite associated with the process of interaction with pathogens, plant growth and stresses response. In this study, 95 complete sequences of the gene for methionine synthase were obtained from 26 specimens of Carapichea ipecacuanha (species of Rubiaceae family, which is popularly known as ipecac) with aims to understand and model evolutionary events that occurred on MetE locus and on the specimens from three disjuncts ranges of occurrence of this species (Atlantic, Amazon and Central America ranges). Sequences of ITS and 5.8S from nrDNA and trnT-trnL spacers from cpDNA were also used for a multiloci population analysis. The sequenced region of MetE covered 3,890pb, where 2,298pb correspond to the coding region (composed of 11 exons) and 1,592pb to the noncoding region (comprising 10 introns and a portion of 3'UTR). In the entire sample were obtained 75 haplotypes and a nucleotide diversity of 0.00899. FST and AMOVA analyzes of MetE indicated no structure in three regions, opposite to detected for other loci analyzed here and in previous studies in which high divergence among the specimens from three ranges was observed. However, in the network haplotype was possible to observe the sympatric distribution of the haplotypes. Possibly, excess of tip haplotypes and the small number of mutational steps among neighboring haplotypes from the three ranges do not allow that the differentiation among these was significant. The Bayesian phylogenetic tree rooted by Carapichea lucida (outgroup) using the coding region sequences showed the formation of clades by geographic regions and provided evidence that sequences from "Atlantic" range is closer to the ancestral state of the gene MetE than sequences from others ranges. The heterogeneity of results for neutrality tests and for values of haplotype and nucleotide diversity at different loci suggests that these are evolving independently and also suggests that demographic events do not explain the pattern found. Thereby, tests for the presence of selection site-by-site were performed for the coding region of the gene MetE and protein alignment against others MetE from diferents plant species available in the Genbank was develop too. Four sites under positive selection were found, but the results show that negative selection (purifying selection) is the main force acting on the locus, favoring the conservation of protein structure and function in all gene copies from C. ipecacuanha. / Metionina sintase independente de cobalamina (MetE) é uma enzima chave na via de síntese de novo de metionina em plantas e também é responsável pela regeneração de S-adenosilmetionina, metabólito associado à vias de interação com patógenos, crescimento vegetal e resposta à estresses. Neste estudo, 95 sequências do gene para metionina sintase foram obtidas a partir de 26 indivíduos de Carapichea ipecacuanha (espécie da família Rubiaceae, popularmente conhecida com ipeca) com objetivo de entender e modelar os eventos evolutivos que ocorreram sobre o locus MetE e sobre os indivíduos de três regiões disjuntas de ocorrência desta espécie (Mata Atlântica, Amazônia e América Central). Sequências de ITS do nrDNA e dos espaçadores trnT-trnL do cpDNA também foram utilizadas para uma análise populacional multiloci. A região sequenciada de MetE abrangeu 3.890pb, sendo 2.298pb correspondentes à região codificante (composta por 11 exons) e 1.592pb à região não codificante (composta por 10 introns e uma parte da região 3 UTR). No total, foram obtidos 75 haplótipos e uma diversidade nucleotídica de 0,00899. As análises FST e AMOVA de MetE indicaram ausência de estruturação em relação às três regiões, diferentemente do observado para os demais loci aqui analisados e em estudos prévios em que elevada divergência entre as três regiões foi observada. Por meio rede de haplótipos foi possível observar a distribuição simpátrica destes; no entanto, o excesso de haplótipos de ponta e o pequeno número de passos mutacionais entre haplótipos vizinhos das três regiões não permitem que a diferenciação destes seja significativa. A árvore filogenética Bayesiana enraizada por Carapichea lucida (grupo externo) utilizando a região codificante também mostrou a formação de clados por regiões geográficas e forneceu evidências de que sequências da Mata Atlântica sejam mais próximas do estado ancestral do gene MetE que as sequências das demais regiões. A heterogeneidade dos resultados encontrados para os testes de neutralidade e para os valores de diversidade nucleotídica e haplotípica nos diferentes loci sugerem que estes estão evoluindo independentemente e que eventos demográficos não explicam o padrão obtido. Assim, análises para a presença de seleção foram conduzidas para a região codificante do gene MetE e também foi realizado alinhamento de proteínas preditas para esta espécie contra proteínas MetE de outras espécies vegetais disponíveis no GenBank. Em quatro sítios houve evidência de seleção positiva, mas os resultados sugerem que a seleção negativa (purificadora) é a principal força atuando sobre o locus, favorecendo a conservação da estrutura e função proteica em todas as cópias gênicas de C. ipecacuanha.

Submitted by PRISCILA SANTOS CASADO null (pri.casado@hotmail.com) on 2017-10-22T18:46:39Z No. of bitstreams: 1 Dissertação Final - Priscila Casado.pdf: 1866790 bytes, checksum: 8d4b39875813bd43bc592acd4b71faef (MD5) / Approved for entry into archive by Luiz Galeffi (luizgaleffi@gmail.com) on 2017-10-26T15:50:29Z (GMT) No. of bitstreams: 1 casado_ps_me_ilha.pdf: 1866790 bytes, checksum: 8d4b39875813bd43bc592acd4b71faef (MD5) / Made available in DSpace on 2017-10-26T15:50:29Z (GMT). No. of bitstreams: 1 casado_ps_me_ilha.pdf: 1866790 bytes, checksum: 8d4b39875813bd43bc592acd4b71faef (MD5) Previous issue date: 2017-08-31 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) / Na primeira etapa da nossa pesquisa propusemos e testamos a utilização de um método baseado em leitor automatizado de microplacas para detectar resistência a fungicidas em populações do fungo hemibiotrófico fitopatogênico Pyricularia graminis-tritici (Pygt) associado à brusone do trigo no Brasil. A disponibilidade de um método mais rápido, acurado e eficiente para detecção de resistência a fungicidas facilitaria a tomada de decisão sobre o manejo químico da brusone do trigo no país. Como há indícios da ineficácia de fungicidas triazóis (DMI) no controle da brusone do trigo e evidências da ocorrência generalizada de resistência a estrobilurinas (QoI) em populações de Pygt no Brasil, inicialmente, testou-se a aplicação do método de microplacas para detecção de resistência a estes dois grupos de fungicidas. O método de microplaca foi acurado em discriminar a variação fenotípica na sensibilidade entre isolados do patógeno aos fungicidas DMI tebuconazol e epoxiconazol e ao fungicida QoI azoxistrobina. Economicamente, ambos os métodos representam um alto investimento inicial. Portanto, é necessário analisar, de forma integral, as necessidades de cada laboratório ou grupo de pesquisa quanto à aquisição dos equipamentos necessários para a fenotipagem visando detecção de resistência a fungicidas em atividades de rotina. Na segunda etapa do nosso estudo buscou-se evidências de resistência ao fungicida SDHI fluxapiroxade em populações de Pygt. Para se estabelecer o baseline de populações do patógeno da brusone do trigo ao fungicida fluxapiroxade, determinou-se a distribuição de frequência de padrões de sensibilidade, com base na EC50 (concentração suficiente para inibir 50% do crescimento micelial fúngico) de cada isolado. Fenótipos extremos de sensibilidade e resistência a fluxapiroxade foram detectados nas diferentes populações do patógeno. Considerando que os fungicidas da segunda geração de SDHI são de alto risco para desenvolvimento de resistência em populações de fipatógenos, é necessária a avaliação do registro e recomendação de uso contínuo do fluxapiroxade, como único princípio ativo numa formulação para o manejo da brusone e de outras doenças do trigo no Brasil. Como estratégia anti-emergência para evitar resistência a SDHI, esse fungicida só poderia ser registrado em co-formulações com outro princípio ativo de baixo risco para resistência. / In the first part of our study, we proposed and tested the use of an automated microplate reader to detect resistance to fungicides in populations of the plant pathogenic hemibiotrophic fungus Pyricularia graminis-tritici (Pygt) associated with the wheat blast disease in Brazil. The availability of a faster, more accurate and efficient method for detecting resistance to fungicies would facilitate the decision making on the chemical management of wheat blast in the country. Because there is indication of inefficacy of triazole fungicides (DMI) for wheat blast control and evidence of the widespread occurrence of strobilurin (QoI) resistance in pathogen populations in Brazil, initially, we tested the application of the microplate method for detecting resistance to these two fungicide groups. The microplate method was accurate in discriminating the phenotypic variation in sensitivity between isolates of the pathogen to the fungicides DMI tebuconazole and epoxiconazole and the fungicide QoI azoxystrobin. Economically, both methods require a high initial investment. Therefore, it is necessary to fully analyze the needs of each laboratory or research group regarding the acquisition of the equipments necessary for phenotyping to detect resistance to fungicides in routine activities. In the second part of our study, using the microplate method, we searched for evidence of resistance to the SDHI fungicide fluxapiroxade in populations of Pygt. To establish the baseline of populations of the wheat blast pathogen to fluxapiroxade, the frequency distribution of sensitivity patterns was determined based on the EC50 (concentration sufficient to inhibit 50% mycelial fungal growth) of each isolate. Extreme phenotypes of sensitivity and resistance to fluxapiroxade were detected in the different pathogen populations. Whereas the second generation of SDHI fungicides are at high risk for development of resistance in populations of plant pathogens, it is necessary to review the labeling and recommendation of continuous application of fluxapyroxad, as the sole active ingredient in any formulation for chemical management of wheat blast and other wheat diseases in Brazil. As an antiemergence strategy to avoid the resistance to SDHI, this fungicide could only be labeled in co-formulations with another active ingredient of low risk for resistance.

Populações nativas da América do Sul apresentam diversidade genética reduzida em relação às demais populações do mundo e alta diferenciação interpopulacional dentro do continente. As pressões seletivas sobre a região cis-reguladora do CCR5 geram uma assinatura de diversidade oposta àquela esperada com base na história demográfica, reduzindo a variação interpopulacional e contribuindo para uma elevação das taxas de polimorfismo. No presente estudo investigamos a interação entre esses processos micro-evolutivos, analisando o polimorfismo da região cis-reguladora de ameríndios da América do Sul e comparando os resultados com aqueles obtidos em outras regiões do mundo. Sequenciamos 927 pares de base (pb) da região controladora do CCR5 em 7 tribos indígenas da região amazônica. Em ameríndios, nenhum haplótipo exclusivo foi encontrado em ameríndios e os dois haplogrupos mais comuns foram de diferentes clusters filogenéticos. A diversidade nucleotídica (&#960;) para a amostra total, foi a mais alta dentre todas as regiões do mundo já estudadas (&#960;=0,0027). A estimativa da diversidade genética populacional para ameríndios, baseada no número de sítios polimórficos (&#952;s), foi similar aos valores encontrados para as populações não-africanas, sendo que em africanos o valor foi praticamente o dobro. Conseqüentemente, foi observado um valor de D de Tajima positivo e estatisticamente significativo (D=2,82, p<0,01), maior do que aquele visto para outras regiões do mundo. Através de comparações empíricas e de simulações coalescentes verificamos que o alto valor de D de Tajima encontrado para ameríndios não é explicável por um recente modelo demográfico de evolução humana. O teste Ewens-Watterson indicou para a amostra ameríndia uma taxa de heterozigose maior do que esperada para um população neutra. O valor de Fst observado para comparação envolvendo asiáticos e ameríndios foi mais baixo que 1000 valores de Fst calculados a partir de 783 microssatélites genotipados em amostras de mesmo tamanho de regiões geográficas similares. Essas caraterísticas corroboram a hipótese de seleção balanceadora na região cis-reguladora do CCR5 de ameríndios da região amazônica brasileira. No presente estudo demonstramos também a existência de uma associação estatisticamente significativa (p<0,01) entre os alelos mantidos em freqüências intermediárias na população humana e regiões ativadoras de splicing localizadas em exons (ESEs). Finalmente, nós acreditamos que um modelo demográfico de evolução humana, complementar ao utilizado neste trabalho, deve ser testado para se refutar a hipótese de que a forte assinatura de seleção balanceadora observada em ameríndios não foi causada por seleção natural que atuou em população ancestral à de indígenas americanos / Native american populations show lower genetic diversity and higher interpopulational genetic differences than populations from other continents. Within South America groups from the non-andean geographic region shows extremely high genetic differentiation. The strong signature of balancing selection observed for 5\' cis-regulatory region of the CCR5 gene at the worldwide scale represents an opposite pattern of genetic variation to that expected by demographic process which acted in Amerindians. To evaluate the impact of complex demographic history on the 5\' cis-regulatory region of CCR5 we resequenced 927 bp of this locus in 62 individuals from different native groups located in the Brazilian Amazonian region. No new haplotype was detected and the two most common haplogroups observed were from different phylogenetic clusters, according with the pattern observed for other human populations. The level of heterozigosity of the total sample, measured by nucleotide diversity (&#960;) was the highest yet described for human populations (&#960;=0,0027) while values based on polymorphic sites (&#952;s) were similar among Amerindians and non-Africans populations. Consequently we observed a positive and significant Tajima\'s D value (D=2,82, p <0,01). This value was higher than all D values observed for the other human populations. Observed summary statistics (&#960; and D) were significantly higher than same statistics estimated from 2000 simulated samples assuming a recently proposed demographic model of human evolution. We also found significant deviations in Ewens-Watterson homozygosity test toward an excess of heterozygosity for Amerindian sample. Another striking feature of CCR5 cis-regulatory region was the low Fst value among populations. The Fst among Asians and Amerindians was an outlier among 1000 Fst values estimated by 783 microssatelites of samples from similar geographic regions and with the same sample sizes as those of the CCR5 5\' cis-regulatory data. These features corroborates the strong signature of balancing selection described for CCR5 5\'cis regulatory region. We also contributed to discussion about functional aspects of CCR5 cis-regulatory region demonstrating a significant association between intermediate frequency SNPs and exonic splicing enhancers (ESEs) regions (p<0,01). Finally, we believe that improved models of demographic human evolution must be tested to refute the hypothesis that the strong signature of balancing selection observed in Amerindians was not caused by a selection pressure which occurred in an ancestral population.

Flexibilidade comportamental é a capacidade de ajustes comportamentais que permitem aos indivíduos responderem a um novo desafio ambiental podendo ser, ou pela modificação de suas respostas, ou por desenvolver novas estratégias diante de um novo estímulo. Esse novo comportamento pode ser alvo da seleção natural, apesar de não ser hereditário. O que é transferida para as próximas gerações, é a capacidade de ajustar o comportamento. Assim, é plausível assumir a existência de uma associação entre capacidades cognitivas e o desafio do ambiente. Neste contexto, espera-se que indivíduos territoriais exibam capacidades cognitivas associadas ao seu habitat específico. Estudos têm frequentemente focado em diferenças nas capacidades de resolução de problemas em animais de espécies e populações distintas. No entanto, o presente estudo compara capacidades de resolução de problemas de lagartos Tropidurus torquatus de uma mesma população que utilizam dois tipos de territórios, os quais diferem no grau de desafio: muros, ambientes desafiadores, e árvores, ambientes permissíveis. Essa espécie exibe um comportamento territorial em que machos dominantes e fêmeas associadas possuem uma fidelidade evidente aos seus territórios. Lagartos foram observados no campo para a caracterização do território e, posteriormente, realizaram testes de resolução de problemas no laboratório. Os resultados demonstraram que lagartos de ambientes desafiadores obtiveram um desempenho superior aos de ambientes permissíveis nos testes de resolução de problemas, e foram mais propensos em exibir uma inovação comportamental durante a resolução dos testes. As interpretações para esses resultados consistem em duas possíveis explicações. Lagartos de muros podem se beneficiar de capacidades cognitivas avançadas, pois a seleção natural teria agido na flexibilidade comportamental e selecionaria indivíduos que poderiam lidar com problemas de uma ambiente desafiador. Alternativamente, é possível que a maioria dos lagartos Tropidurus torquatus apresentem melhores capacidades de resolução de problema, porém esta característica estaria críptica, e apenas a interação entre o indivíduo e a vivência em ambientes desafiadores (muros) induziria o surgimento de melhores capacidades cognitivas. Os resultados desse estudo provavelmente indicam que diferenças nas capacidades de resolução de problemas e habilidades cognitivas entre indivíduos poderiam estar sob efeitos de regimes seletivos antes mesmo de ocorrer alguma segregação populacional. / Behavioral flexibility is the ability of behavioral adjustments that allow individuals to respond to a new environmental challenge by either modifying their responses or developing new strategies when facing novel stimuli. This new behavior may be target of natural selection, despite not being heritable. What is transferred to the next generation is the capacity to adjust behavior. Thus, it is plausible to assume the existence of an association between cognitive capabilities and environmental severity. In this context, it is expected that territorial individuals exhibit cognitive capabilities associated with their specific habitat. Studies have often been focused on differences of problem-solving capabilities in animals from different species and populations. Nevertheless, the present study compares problem-solving abilities of Tropidurus torquatus lizards from the same population that use two types of territories that differ in the level of harshness: house walls, a harsh environment, and trees, a mild environment. This species exhibits a territorial system where dominant males and associated females have an evident fidelity to their territories. Lizards were observed in the field for territory characterization and then performed problem-solving tests in the laboratory. The results show that lizards from harsh environments outperformed those from mild environments in the problem-solving tasks and were more likely to show a behavioral innovation when solving the task. Interpretations for these results resides on two possible explanations. Lizards from walls might benefit from enhanced cognitive capabilities because natural selection acted on behavioral flexibility and selected individuals that can deal with problems of a harsh environment. Alternatively, it is possible that the majority of Tropidurus torquatus lizards exhibit enhanced problem-solving capability, but this characteristic might be cryptic, and only the interaction between the individual and living in a harsh environment (walls) induces the emergence of enhanced cognitive capabilities. These findings probably indicate that differences in problem-solving and cognitive capabilities might be exposed to specific selective regimes even before population segregation occurs.

No description available.

In this thesis I dealt with on the relationship between religious belief and acceptance or rejection the theories about origin of life and evolution-it means evolution by natural selection and theories such as creationism and Intelligent design. Respondents of the survey were students of universities and high-schools focused on natural sciences. The review of literature consists of chapters dedicated to evolution, creationism, ID, religion and surveys on this topic. In the research part of the thesis there is own analysis on the sample of students. The research was realized in the form of survey on students of universities, mostly and the Faculty of Science of Charles University in Prague and on the high schools, which are the teaching schools of this faculty. The benefit of this thesis is finding the relationship between their religious belief and their opinion about theories origin of life, evolution by natural selection, creationism and ID. Keywords: Evolution by natural selection, Intelligent design, Creationism, Religion, Belief

In this thesis I dealt with on the relationship between religious belief and acceptance or rejection the theories about origin of life and evolution-it means evolution by natural selection and theories such as creationism and Intelligent design. Respondents of the survey were students of universities and high-schools focused on natural sciences. The review of literature consists of chapters dedicated to evolution, creationism, Intelligent design, religion and surveys on this topic. In the research part of the thesis there is own analysis on the sample of students. The research was realized in the form of survey on students of universities, mostly and the Faculty of Science of Charles University and on the high schools, which are the teaching schools of this faculty. The benefit of this thesis is finding the relationship between their religious belief and their opinion about theories origin of life, evolution by natural selection, creationism and Intelligent design. Keywords: Evolution by natural selection, Intelligent design, Creationism, Religion, Belief

No description available.