The genetic content of chromosomal inversions across a wide latitudinal gradient

There is increasing evidence regarding the role of chromosomal inversions in relevant biological processes such as local adaptation and speciation. A classic example of the adaptive role of chromosomal polymorphisms is given by the clines of inversion frequencies in Drosophila subobscura, repeatable across continents. Nevertheless, not much is known about the molecular variation associated with these polymorphisms. We characterized the genetic content of ca. 600 individuals from nine European populations following a latitudinal gradient by analysing 19 microsatellite loci from two autosomes (J and U) and the sex chromosome (A), taking into account their chromosomal inversions. Our results clearly demonstrate the molecular genetic uniformity within a given chromosomal inversion across a large latitudinal gradient, particularly from Groningen (Netherlands) in the north to Málaga (Spain) in the south, experiencing highly diverse environmental conditions. This low genetic differentiation within the same gene arrangement across the nine European populations is consistent with the local adaptation hypothesis for th evolutionof chromosomal polymorphisms. We also show the effective role of chromosomal inversions in maintaining different genetic pools within these inverted genomic regions even in the presence of high gene flow. Inversions represent thus an important barrier to gene flux and can help maintain specific allelic combinations with positive effects on fitness. Consistent patterns of microsatellite allele-inversion linkage disequilibrium particularly in loci within inversions were also observed. Finally, we identified areas within inversions presenting clinal variation that might be under selection.

Consequences of lineage-specific gene loss on functional evolution of surviving paralogs: ALDH1A and retinoic acid signaling in vertebrate genomes

Genome duplications increase genetic diversity and may facilitate the evolution of gene subfunctions. Little attention, however, has focused on the evolutionary impact of lineage-specific gene loss. Here, we show that identifying lineage-specific gene loss after genome duplication is important for understanding the evolution of gene subfunctions in surviving paralogs and for improving functional connectivity among human and model organism genomes. We examine the general principles of gene loss following duplication, coupled with expression analysis of the retinaldehyde dehydrogenase Aldh1a gene family during retinoic acid signaling in eye development as a case study. Humans have three ALDH1A genes, but teleosts have just one or two. We used comparative genomics and conserved syntenies to identify loss of ohnologs (paralogs derived from genome duplication) and to clarify uncertain phylogenies. Analysis showed that Aldh1a1 and Aldh1a2 form a clade that is sister to Aldh1a3-related genes. Genome comparisons showed secondarily loss of aldh1a1 in teleosts, revealing that Aldh1a1 is not a tetrapod innovation and that aldh1a3 was recently lost in medaka, making it the first known vertebrate with a single aldh1a gene. Interestingly, results revealed asymmetric distribution of surviving ohnologs between co-orthologous teleost chromosome segments, suggesting that local genome architecture can influence ohnolog survival. We propose a model that reconstructs the chromosomal history of the Aldh1a family in the ancestral vertebrate genome, coupled with the evolution of gene functions in surviving Aldh1a ohnologs after R1, R2, and R3 genome duplications. Results provide evidence for early subfunctionalization and late subfunction-partitioning and suggest a mechanistic model based on altered regulation leading to heterochronic gene expression to explain the acquisition or modification of subfunctions by surviving ohnologs that preserve unaltered ancestral developmental programs in the face of gene loss.

GENETIC ARCHITECTURE, EVOLUTION AND MAINTENANCE OF A COLOR CLINE AT A CONTINENTAL SCALE

On a geological time scale the conditions on earth are very variable and biological patterns (for example the distributions of species) are very dynamic. Understanding large scale patterns of variation observed today thus requires a deep understanding of the historical factors that drove their evolution. In this thesis, we reevaluated the evolution and maintenance of a continental color cline observed in the European barn owl (Tyto alba) using population genetic tools. The colour cline spans from south-est Europe where most individual have pure white underparts to north and east Europe where most individuals have rufous-brown underparts. Our results globally showed that the old scenario, stipulating that the color cline evolved by secondary contact of two color morphs (white and rufous) that evolved in allopatry during the last ice age has to be revised. We collected samples of about 700 barn owls from the Western Palearctic to establish the first population genetic data set for this species. Individuals were genotyped at 22 microsatellites markers, at one mitochondrial gene, and at a candidate color gene. The color of each individuals was assessed and their sex determined by molecular methods. We first showed that the genetic variation in Western Europe is very limited compared to the heritable color variation. We found no evidences of different glacial lineages, and showed that selection must be involved in the maintenance of the color cline (chapter 1). Using computer simulations, we demonstrated that the post-glacial colonization of Europe occurred from the Iberian Peninsula and that the color cline could not have evolved by neutral demographic processes during this colonization (chapter 2). Finally we reevaluated the whole history of the establishment of the Western Palearctic variation of the barn owl (chapter 3): This study showed that all Western European barn owls descend from white barn owls phenotypes from the Middle East that colonized the Iberian Peninsula via North-Africa. Following the end of the last ice age (20'000 years ago), these white barn owls colonized Western Europe and under selection a novel rufous phenotype evolved (during or after the colonization). An important part of the color variation could be explained by a single mutation in the melanocortin-1-receptor (MC1R) gene that appeared during or after the colonization. The colonization of Europe reached until Greece, where the rufous birds encountered white ones (which reached Greece from the Middle East over the Bosporus) in a secondary contact zone. Our analyses show that white and rufous barn owls in Greece interbreed only to a limited extent. This suggests that barn owls are at the verge of becoming two species in Greece and demonstrates that European barn owls represent an incipient ring species around the Mediterranean. The revisited history of the establishment of the European barn owl color cline makes this model system remarkable for several aspects. It is a very clear example of strong local adaptation that can be achieved despite high gene flow (strong color and MC1R differentiation despite almost no neutral genetic differentiation). It also offers a wonderful model system to study the interactions between colonization processes and selection processes which have, for now, been remarkably understudied despite their potentially ubiquitous importance. Finally it represents a very interesting case in the speciation continuum and appeals for further studying the amount of gene flow that occurs between the color morphs in Greece. -- Sur l'échelle des temps géologiques, les conditions sur terre sont très variables et les patrons biologiques (telle que la distribution des espèces) sont très dynamiques. Si l'on veut comprendre des patrons que l'on peut observer à large échelle aujourd'hui, il est nécessaire de d'abord comprendre les facteurs historiques qui ont gouverné leur établissement. Dans cette thèse, nous allons réévaluer, grâce à des outils modernes de génétique des populations, l'évolution et la maintenance d'un cline de couleur continental observé chez l'effraie des clochers européenne (Tyto alba). Globalement, nos résultats montrent que le scenario accepté jusqu'à maintenant, qui stipule que le cline de couleur a évolué à partir du contact secondaire de deux morphes de couleur (blanches et rousses) ayant évolué en allopatrie durant les dernières glaciations, est à revoir. Afin de constituer le premier jeu de données de génétique des populations pour cette espèce, nous avons récolté des échantillons d'environ 700 effraies de l'ouest Paléarctique. Nous avons génotypé tous les individus à 22 loci microsatellites, sur un gène mitochondrial et sur un autre gène participant au déterminisme de la couleur. Nous avons aussi mesuré la couleur de tous les individus et déterminé leur sexe génétiquement. Nous avons tout d'abord pu montrer que la variation génétique neutre est négligeable en comparaison avec la variation héritable de couleur, qu'il n'existe qu'une seule lignée européenne et que de la sélection doit être impliquée dans le maintien du cline de couleur (chapitre 1). Grâce à des simulations informatiques, nous avons démontré que l'ensemble de l'Europe de l'ouest a été recolonisé depuis la Péninsule Ibérique après les dernières glaciations et que le cline de couleur ne peut pas avoir évolué par des processus neutre durant cette colonisation (chapitre 2). Finalement, nous avons réévalué l'ensemble de l'histoire postglaciaire de l'espèce dans l'ouest Paléarctique (chapitre 3): l'ensemble des effraies du Paléarctique descendent d'effraie claire du Moyen-Orient qui ont colonisé la péninsule ibérique en passant par l'Afrique du nord. Après la fin de la dernière glaciation (il y a 20'000 ans), ces effraies claires ont colonisé l'Europe de l'ouest et ont évolués par sélection le phénotype roux (durant ou après la colonisation). Une part importante de la variation de couleur peut être expliquée par une mutation sur le gène MC1R qui est apparue durant ou juste après la colonisation. Cette vague de colonisation s'est poursuivie jusqu'en Grèce où ces effraies rousses ont rencontré dans une zone de contact secondaire des effraies claires (qui sont remontées en Grèce depuis le Moyen-Orient via le Bosphore). Nos analyses montrent que le flux de gènes entre effraies blanches et rousses est limité en Grèce, ce qui suggère qu'elles sont en passe de former deux espèces et ce qui montre que les effraies constituent un exemple naissant de spéciation en anneaux autour de la Méditerranée. L'histoire revisitée des effraies des clochers de l'ouest Paléarctique en fait un système modèle remarquable pour plusieurs aspects. C'est un exemple très claire de forte adaptation locale maintenue malgré un fort flux de gènes (différenciation forte de couleur et sur le gène MC1R malgré presque aucune structure neutre). Il offre également un très bon système pour étudier l'interaction entre colonisation et sélection, un thème ayant été remarquablement peu étudié malgré son importance. Et il offre finalement un cas très intéressant dans le « continuum de spéciation » et il serait très intéressant d'étudier plus en détail l'importance du flux de gènes entre les morphes de couleur en Grèce.

Combining genetic and demographic data for the conservation of a Mediterranean habitat-forming species

The integration of ecological and evolutionary data is highly valuable for conservation planning. However, it has been rarely used in the marine realm, where the adequate design of marine protected areas (MPAs) is urgently needed. Here, we examined the interacting processes underlying the patterns of genetic structure and demographic strucuture of a highly vulnerable Mediterranean habitat-forming species (i.e. Paramuricea clavata (Risso, 1826)), with particular emphasis on the processes of contemporary dispersal, genetic drift, and colonization of a new population. Isolation by distance and genetic discontinuities were found, and three genetic clusters were detected; each submitted to variations in the relative impact of drift and gene flow. No founder effect was found in the new population. The interplay of ecology and evolution revealed that drift is strongly impacting the smallest, most isolated populations, where partial mortality of individuals was highest. Moreover, the eco-evolutionary analyses entailed important conservation implications for P. clavata. Our study supports the inclusion of habitat-forming organisms in the design of MPAs and highlights the need to account for genetic drift in the development of MPAs. Moreover, it reinforces the importance of integrating genetic and demographic data in marine conservation.

Conservation Genetics of Exploited Finnish Salmonid Fishes

Genetic diversity is one of the levels of biodiversity that the World Conservation Union (IUCN) has recognized as being important to preserve. This is because genetic diversity is fundamental to the future evolution and to the adaptive flexibility of a species to respond to the inherently dynamic nature of the natural world. Therefore, the key to maintaining biodiversity and healthy ecosystems is to identify, monitor and maintain locally-adapted populations, along with their unique gene pools, upon which future adaptation depends. Thus, conservation genetics deals with the genetic factors that affect extinction risk and the genetic management regimes required to minimize the risk. The conservation of exploited species, such as salmonid fishes, is particularly challenging due to the conflicts between different interest groups. In this thesis, I conduct a series of conservation genetic studies on primarily Finnish populations of two salmonid fish species (European grayling, Thymallus thymallus, and lake-run brown trout, Salmo trutta) which are popular recreational game fishes in Finland. The general aim of these studies was to apply and develop population genetic approaches to assist conservation and sustainable harvest of these populations. The approaches applied included: i) the characterization of population genetic structure at national and local scales; ii) the identification of management units and the prioritization of populations for conservation based on evolutionary forces shaping indigenous gene pools; iii) the detection of population declines and the testing of the assumptions underlying these tests; and iv) the evaluation of the contribution of natural populations to a mixed stock fishery. Based on microsatellite analyses, clear genetic structuring of exploited Finnish grayling and brown trout populations was detected at both national and local scales. Finnish grayling were clustered into three genetically distinct groups, corresponding to northern, Baltic and south-eastern geographic areas of Finland. The genetic differentiation among and within population groups of grayling ranged from moderate to high levels. Such strong genetic structuring combined with low genetic diversity strongly indicates that genetic drift plays a major role in the evolution of grayling populations. Further analyses of European grayling covering the majority of the species’ distribution range indicated a strong global footprint of population decline. Using a coalescent approach the beginning of population reduction was dated back to 1 000-10 000 years ago (ca. 200-2 000 generations). Forward simulations demonstrated that the bottleneck footprints measured using the M ratio can persist within small populations much longer than previously anticipated in the face of low levels of gene flow. In contrast to the M ratio, two alternative methods for genetic bottleneck detection identified recent bottlenecks in six grayling populations that warrant future monitoring. Consistent with the predominant role of random genetic drift, the effective population size (N_e) estimates of all grayling populations were very low with the majority of N_e estimates below 50. Taken together, highly structured local populations, limited gene flow and the small N_e of grayling populations indicates that grayling populations are vulnerable to overexploitation and, hence, monitoring and careful management using the precautionary principles is required not only in Finland but throughout Europe. Population genetic analyses of lake-run brown trout populations in the Inari basin (northernmost Finland) revealed hierarchical population structure where individual populations were clustered into three population groups largely corresponding to different geographic regions of the basin. Similar to my earlier work with European grayling, the genetic differentiation among and within population groups of lake-run brown trout was relatively high. Such strong differentiation indicated that the power to determine the relative contribution of populations in mixed fisheries should be relatively high. Consistent with these expectations, high accuracy and precision in mixed stock analysis (MSA) simulations were observed. Application of MSA to indigenous fish caught in the Inari basin identified altogether twelve populations that contributed significantly to mixed stock fisheries with the Ivalojoki river system being the major contributor (70%) to the total catch. When the contribution of wild trout populations to the fisheries was evaluated regionally, geographically nearby populations were the main contributors to the local catches. MSA also revealed a clear separation between the lower and upper reaches of Ivalojoki river system – in contrast to lower reaches of the Ivalojoki river that contributed considerably to the catch, populations from the upper reaches of the Ivalojoki river system (>140 km from the river mouth) did not contribute significantly to the fishery. This could be related to the available habitat size but also associated with a resident type life history and increased cost of migration. The studies in my thesis highlight the importance of dense sampling and wide population coverage at the scale being studied and also demonstrate the importance of critical evaluation of the underlying assumptions of the population genetic models and methods used. These results have important implications for conservation and sustainable fisheries management of Finnish populations of European grayling and brown trout in the Inari basin.

Genotypic diversity among brazilian isolates of Macrophomina phaseolina revealed by RAPD

Macrophomina phaseolina has been considered one of the most prevalent soybean (Glycine max) pathogens in Brazil. No genetic resistance has been determined in soybean and very little is known about the genetic diversity of this pathogen in tropical and sub-tropical regions. Fifty-five isolates from soybean roots were collected in different regions and analyzed through RAPD for genetic diversity. The UPGMA cluster analysis for 74 loci scored permitted identification of three divergent groups with an average similarity of 99%, 92% and 88%, respectively. The three groups corresponded to 5.45%, 59.95% and 34.6%, respectively of all isolates used. A single plant had three different haplotypes, while 10.9% of the analyzed plants had two different haplotypes. In another study the genetic similarity was evaluated among isolates from different hosts [soybean, sorghum (Sorghum bicolor), sunflower (Helianthus annuus), cowpea (Vigna unguiculata), corn (Zea mays) and wheat (Triticum aestivum)] as well as two soil samples from native areas. Results showed that more divergent isolates originated from areas with a single crop. Isolates from areas with crop rotation were less divergent, showing high similarity values and consequently formed the largest group. Amplification of the ITS region using primers ITS1 and ITS4 produced only one DNA fragment of 620 bp. None of the isolates were differentiated through PCR-RFLP. Our results demonstrated genetic variability among Brazilian isolates of M. phaseolina and showed that one single root can harbor more than one haplotype. Moreover, cultivation with crop rotation tends to induce less specialization of the pathogen isolates. Knowledge of this variation may be useful in screening soybean genotypes for resistance to charcoal rot.

REPRODUCTIVE SYSTEM AND POLLEN FLOW IN PROGENIES OF Qualea grandiflora Mart., A TYPICAL SPECIES OF THE BRAZILIAN CERRADO

This study analyzed the reproductive system and the pollen dispersion pattern of Qualea grandiflora progenies. This is a typical species from the Brazilian Cerrado about which there are not too many studies from the genetics point of view. The study was conducted in an area of 2.2 hectares located in the Conservation Unit managed by the Forest Institute of the state of São Paulo, Brazil (Assis State Forest). Total genomic DNA of 300 seeds from 25 plants (12 seeds from each plant) was extracted and amplified using specific primers to obtain microsatellite markers. Results showed that selfing is frequent among adults and progenies, and the species reproduces by outcrossing between related and unrelated individuals (0.913). The single-locus outcrossing rate was 0.632, which indicates that mating between unrelated individuals is more frequent than between related plants. The selfing rate was low (0.087), that is, the species is allogamous and self-fertilization is reduced. About 35% of the plants in the progenies were full-sibs, and about 57%, half-sibs. Besides, about 8% of the progenies were selfing siblings. The genetic differentiation coefficient within progenies was 0.139, whereas the fixation rate was about 27%. The estimate of the effective size revealed that the genetic representativeness of descent was lower than expected in random mating progenies: The analyzed samples corresponded to only 13.2 individuals of an ideal panmictic population. In environmental recovery programs, seeds, preferably from different fruits, should be collected from 95 trees to preserve the genetic diversity of the species.

Genotype characterization of Haematobia irritans from different Brazilian geographic regions based on randomly amplified polymorphic DNA (RAPD) analysis

Blood-sucking diptera are important parasites in bovine production systems, especially regarding confinement conditions. Haematobia irritans, the horn fly, is one of the most troublesome species within bovine production systems, due to the intense stress imposed to the animals. An important aspect while studying the variability within a species is the study of the geographic structure of its populations and, attempting to find out the genetic flow of Brazilian populations of horn fly, the RAPD technique, which is suited for this purpose, has been used. The use of molecular markers generated from RAPD made it possible to identify the geographic origin of samples from different Brazilian geographic regions, as well as to estimate the genotypic flow among the different Brazilian populations of the horn fly.

Sequence analysis of the 5′ third of glycoprotein C gene of South American bovine herpesviruses 1 and 5

Bovine herpesviruses 1 (BoHV-1) and 5 (BoHV-5) share high genetic and antigenic similarities, but exhibit marked differences in tissue tropism and neurovirulence. The amino-terminal region of glycoprotein C (gC), which is markedly different in each of the viruses, is involved in virus binding to cellular receptors and in interactions with the immune system. This study investigated the genetic and antigenic differences of the 5′ region of the gC (5′ gC) gene (amino-terminal) of South American BoHV-1 (n=19) and BoHV-5 (n=25) isolates. Sequence alignments of 374 nucleotides (104 amino acids) revealed mean similarity levels of 97.3 and 94.2% among BoHV-1 gC (gC1), respectively, 96.8 and 95.6% among BoHV-5 gC (gC5), and 62 and 53.3% between gC1 and gC5. Differences included the absence of 40 amino acid residues (27 encompassing predicted linear epitopes) scattered throughout 5′ gC1 compared to 5′ gC5. Virus neutralizing assays testing BoHV-1 and BoHV-5 antisera against each isolate revealed a high degree of cross-neutralization between the viruses, yet some isolates were neutralized at very low titers by heterologous sera, and a few BoHV-5 isolates reacted weakly with either sera. The virus neutralization differences observed within the same viral species, and more pronounced between BoHV-1 and BoHV-5, likely reflect sequence differences in neutralizing epitopes. These results demonstrate that the 5′ gC region is well conserved within each viral species but is divergent between BoHV-1 and BoHV-5, likely contributing to their biological and antigenic differences.

Disease-Gene Association Using Genetic Programming

As a result of mutation in genes, which is a simple change in our DNA, we will have undesirable phenotypes which are known as genetic diseases or disorders. These small changes, which happen frequently, can have extreme results. Understanding and identifying these changes and associating these mutated genes with genetic diseases can play an important role in our health, by making us able to find better diagnosis and therapeutic strategies for these genetic diseases. As a result of years of experiments, there is a vast amount of data regarding human genome and different genetic diseases that they still need to be processed properly to extract useful information. This work is an effort to analyze some useful datasets and to apply different techniques to associate genes with genetic diseases. Two genetic diseases were studied here: Parkinson’s disease and breast cancer. Using genetic programming, we analyzed the complex network around known disease genes of the aforementioned diseases, and based on that we generated a ranking for genes, based on their relevance to these diseases. In order to generate these rankings, centrality measures of all nodes in the complex network surrounding the known disease genes of the given genetic disease were calculated. Using genetic programming, all the nodes were assigned scores based on the similarity of their centrality measures to those of the known disease genes. Obtained results showed that this method is successful at finding these patterns in centrality measures and the highly ranked genes are worthy as good candidate disease genes for being studied. Using standard benchmark tests, we tested our approach against ENDEAVOUR and CIPHER - two well known disease gene ranking frameworks - and we obtained comparable results.

Reproduction et échanges génétiques horizontaux chez les champignons mycorhiziens à arbuscules

3 vidéos sont dans des fichiers complémentaires à ce mémoire

Sélection et polymorphisme chez des grenouilles mimétiques du Pérou (Dendrobatidae)

La diversification des signaux aposématiques dans un cadre de mimétisme müllérien est un phénomène intrigant. Alors que la théorie relative à l'aposématisme et au mimétisme suggère l'évolution vers un signal aposématique unique, d'impressionnantes variations peuvent être observées entre les populations, et cela à petite échelle spatiale. Il a été supposé que la variation spatiale des pressions de sélection engendrées par différents prédateurs puisse être à l'origine de ce phénomène. Afin de tester cette hypothèse, nous avons étudié la transition entre deux systèmes géographiques caractérisés par des patrons aposématiques distincts chez des grenouilles mimétiques et toxiques du nord du Pérou (Dendrobatidae) en combinant les outils de génétique des populations aux outils écologiques. Dans chacun de ces systèmes, Ranitomeya imitator vit en sympatrie avec R. ventrimaculata ou R. variabilis. Il s'agit du principal exemple empirique suggérant que dans un cadre de mimétisme müllérien, il n'y a pas convergence des signaux aposématiques des deux espèces, mais plutôt convergence unidirectionnelle où R. imitator, étant polymorphe, imite des espèces monomorphes avec lesquelles elle est sympatrique. Premièrement, les résultats réfutent les prémisses qui suggèrent que R. imitator converge vers le signal aposématique d’une autre espèce. La haute similarité génétique entre les espèces modèles suggère qu'elles ont divergé plus récemment que les populations de R. imitator ou qu'elles sont encore connectées par du flux génique. Ces résultats indiquent que ces espèces ont été identifiées à tort comme des espèces différentes. De fait, l'identification de l'espèce imitatrice basée sur la variabilité phénotypique est invalidée dans ce système puisque R. imitator et R. variabilis/ventrimaculata démontrent la même variabilité. Deuxièmement, nos résultats démontrent que la prédation varie spatialement, autant en intensité qu'en direction, créant ainsi un paysage hétérogène de pressions de sélection. Ainsi, de fortes pressions de prédation stabilisatrice permettent le maintien de l'organisation géographique de différents signaux aposématiques et expliquent l'uniformité de ces signaux ainsi que les relations mimétiques. Par contre, le relâchement temporaire des pressions de prédation permet l'apparition de nouveaux phénotypes aposématiques via les processus évolutifs neutres, conduisant à un haut polymorphisme au niveau de ces populations. L'interaction de ces modes sélectifs nous a permis de démontrer pour la première fois comment la théorie évolutive de Wright (shifting balance theory) permet la diversification adaptative dans un système naturel. Pour conclure, cette étude a permis de mettre en évidence à quel point les systèmes de mimétisme müllérien peuvent être dynamiques. L'alternance spatiale entre les processus évolutifs neutres et la sélection naturelle permet l'émergence de nouveaux phénotypes aposématiques à une échelle locale, ainsi que l'apparition d'une organisation géographique des signaux d'avertissement et des relations de mimétisme müllérien.

Structure génétique de populations montréalaises de salamandres cendrées (Plethodon cinereus) et de salamandres à points bleus (Ambystoma laterale)

Thèse numérisée par la Division de la gestion de documents et des archives de l'Université de Montréal

Studies on the Nature and Chemistry of Sediments and Water of Periyar and Chalakudy Rivers, Kerala, India

The present study is an attempt to address issues related to sediment properties like texture, mineralogy and geochemistry as well as water quality of two important rivers of central Kerala-the Periyar and the Chalakudy rivers. The main objectives of the study are to investigate the textural and mineralogical characteristics as well as transportation and depositional mechanisms of the sediments of Periyar and Chalakudy rivers, to find out the geochemical variability of organic carbon, phosphorus and certain major (Na,K,Ca and Mg) and minor/trace(Mn,Pb,Ni,Cr, and Zn) elements in the bulk sediments and mud fraction of these rivers, to evaluate the status of heavy metal pollution registered in the sediments of these rivers, to assess the physico-chemical characteristics and water quality of Periyar and Chalakudy rivers and to estimate the dissolved nutrient flux through the Periyar and Chalakudy rivers into the receiving coastal waters. The granulometric characteristics as well as statistical parameters of the sediments of Periyar and Chalakudy rivers depend on the flow pattern controlled by the gradient of the terrain. Compared to Periyar, fluctuations in the dispersal of particles are more in Chalakudy river. In Periyar river, the P and Fe in bulk sediments show a positive correlation with C-org, while in Chalakudy river, both the elements are related to THM concentration. In general, C-org, Fe and P Shows an increasing trend downstream. In Periyar river, the P and Fe in bulk sediments show a positive correlation with C-org, while in Chalakudy river, both the elements are related to THM concentration. Among these two rivers, the pollution of water is several fold higher in Periyar river due to influx due to influx of considerable quantity of liquid and solid wastes of industrial/domestic/urban origin. Nutrient analysis reveals 2-3 times increase in N and P during monsoon season whereas SiO2-Si shows a decreasing trend.