High-density linkage maps fail to detect any genetic component to sex determination in a Rana temporaria family.

Sex chromosome differentiation in Rana temporaria varies strikingly among populations or families: whereas some males display well-differentiated Y haplotypes at microsatellite markers on linkage group 2 (LG2 ), others are genetically undistinguishable from females. We analysed with RADseq markers one family from a Swiss lowland population with no differentiated sex chromosomes, and where sibship analyses had failed to detect any association between the phenotypic sex of progeny and parental haplotypes. Offspring were reared in a common tank in outdoor conditions and sexed at the froglet stage. We could map a total of 2177 SNPs (1123 in the mother, 1054 in the father), recovering in both adults 13 linkage groups (= chromosome pairs) that were strongly syntenic to Xenopus tropicalis despite > 200 My divergence. Sexes differed strikingly in the localization of crossovers, which were uniformly distributed in the female but limited to chromosome ends in the male. None of the 2177 markers showed significant association with offspring sex. Considering the very high power of our analysis, we conclude that sex determination was not genetic in this family; which factors determined sex remain to be investigated.

Genetic associations at 53 loci highlight cell types and biological pathways relevant for kidney function.

Reduced glomerular filtration rate defines chronic kidney disease and is associated with cardiovascular and all-cause mortality. We conducted a meta-analysis of genome-wide association studies for estimated glomerular filtration rate (eGFR), combining data across 133,413 individuals with replication in up to 42,166 individuals. We identify 24 new and confirm 29 previously identified loci. Of these 53 loci, 19 associate with eGFR among individuals with diabetes. Using bioinformatics, we show that identified genes at eGFR loci are enriched for expression in kidney tissues and in pathways relevant for kidney development and transmembrane transporter activity, kidney structure, and regulation of glucose metabolism. Chromatin state mapping and DNase I hypersensitivity analyses across adult tissues demonstrate preferential mapping of associated variants to regulatory regions in kidney but not extra-renal tissues. These findings suggest that genetic determinants of eGFR are mediated largely through direct effects within the kidney and highlight important cell types and biological pathways.

Genome-wide Association Studies Identify Genetic Loci Associated With Albuminuria in Diabetes.

Elevated concentrations of albumin in the urine, albuminuria, are a hallmark of diabetic kidney disease and are associated with an increased risk for end-stage renal disease and cardiovascular events. To gain insight into the pathophysiological mechanisms underlying albuminuria, we conducted meta-analyses of genome-wide association studies and independent replication in up to 5,825 individuals of European ancestry with diabetes and up to 46,061 without diabetes, followed by functional studies. Known associations of variants in CUBN, encoding cubilin, with the urinary albumin-to-creatinine ratio (UACR) were confirmed in the overall sample (P = 2.4 × 10(-10)). Gene-by-diabetes interactions were detected and confirmed for variants in HS6ST1 and near RAB38/CTSC. Single nucleotide polymorphisms at these loci demonstrated a genetic effect on UACR in individuals with but not without diabetes. The change in the average UACR per minor allele was 21% for HS6ST1 (P = 6.3 × 10(-7)) and 13% for RAB38/CTSC (P = 5.8 × 10(-7)). Experiments using streptozotocin-induced diabetic Rab38 knockout and control rats showed higher urinary albumin concentrations and reduced amounts of megalin and cubilin at the proximal tubule cell surface in Rab38 knockout versus control rats. Relative expression of RAB38 was higher in tubuli of patients with diabetic kidney disease compared with control subjects. The loci identified here confirm known pathways and highlight novel pathways influencing albuminuria.

Hybridization Capture Using RAD Probes (hyRAD), a New Tool for Performing Genomic Analyses on Collection Specimens.

In the recent years, many protocols aimed at reproducibly sequencing reduced-genome subsets in non-model organisms have been published. Among them, RAD-sequencing is one of the most widely used. It relies on digesting DNA with specific restriction enzymes and performing size selection on the resulting fragments. Despite its acknowledged utility, this method is of limited use with degraded DNA samples, such as those isolated from museum specimens, as these samples are less likely to harbor fragments long enough to comprise two restriction sites making possible ligation of the adapter sequences (in the case of double-digest RAD) or performing size selection of the resulting fragments (in the case of single-digest RAD). Here, we address these limitations by presenting a novel method called hybridization RAD (hyRAD). In this approach, biotinylated RAD fragments, covering a random fraction of the genome, are used as baits for capturing homologous fragments from genomic shotgun sequencing libraries. This simple and cost-effective approach allows sequencing of orthologous loci even from highly degraded DNA samples, opening new avenues of research in the field of museum genomics. Not relying on the restriction site presence, it improves among-sample loci coverage. In a trial study, hyRAD allowed us to obtain a large set of orthologous loci from fresh and museum samples from a non-model butterfly species, with a high proportion of single nucleotide polymorphisms present in all eight analyzed specimens, including 58-year-old museum samples. The utility of the method was further validated using 49 museum and fresh samples of a Palearctic grasshopper species for which the spatial genetic structure was previously assessed using mtDNA amplicons. The application of the method is eventually discussed in a wider context. As it does not rely on the restriction site presence, it is therefore not sensitive to among-sample loci polymorphisms in the restriction sites that usually causes loci dropout. This should enable the application of hyRAD to analyses at broader evolutionary scales.

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.

Distribution and population genetic variation of cryptic species of the Alpine mayfly Baetis alpinus (Ephemeroptera: Baetidae) in the Central Alps.

BACKGROUND: Many species contain evolutionarily distinct groups that are genetically highly differentiated but morphologically difficult to distinguish (i.e., cryptic species). The presence of cryptic species poses significant challenges for the accurate assessment of biodiversity and, if unrecognized, may lead to erroneous inferences in many fields of biological research and conservation. RESULTS: We tested for cryptic genetic variation within the broadly distributed alpine mayfly Baetis alpinus across several major European drainages in the central Alps. Bayesian clustering and multivariate analyses of nuclear microsatellite loci, combined with phylogenetic analyses of mitochondrial DNA, were used to assess population genetic structure and diversity. We identified two genetically highly differentiated lineages (A and B) that had no obvious differences in regional distribution patterns, and occurred in local sympatry. Furthermore, the two lineages differed in relative abundance, overall levels of genetic diversity as well as patterns of population structure: lineage A was abundant, widely distributed and had a higher level of genetic variation, whereas lineage B was less abundant, more prevalent in spring-fed tributaries than glacier-fed streams and restricted to high elevations. Subsequent morphological analyses revealed that traits previously acknowledged as intraspecific variation of B. alpinus in fact segregated these two lineages. CONCLUSIONS: Taken together, our findings indicate that even common and apparently ecologically well-studied species may consist of reproductively isolated units, with distinct evolutionary histories and likely different ecology and evolutionary potential. These findings emphasize the need to investigate hidden diversity even in well-known species to allow for appropriate assessment of biological diversity and conservation measures.

Population connectivity buffers genetic diversity loss in a seabird

Background Ancient DNA has revolutionized conservation genetic studies as it allows monitoring of the genetic variability of species through time and predicting the impact of ecosystems" threats on future population dynamics and viability. Meanwhile, the consequences of anthropogenic activities and climate change to island faunas, particularly seabirds, remain largely unknown. In this study, we examined temporal changes in the genetic diversity of a threatened seabird, the Cory"s shearwater (Calonectris borealis). Findings We analysed the mitochondrial DNA control region of ancient bone samples from the late-Holocene retrieved from the Canary archipelago (NE Atlantic) together with modern DNA sequences representative of the entire breeding range of the species. Our results show high levels of ancient genetic diversity in the Canaries comparable to that of the extant population. The temporal haplotype network further revealed rare but recurrent long-distance dispersal between ocean basins. The Bayesian demographic analyses reveal both regional and local population size expansion events, and this is in spite of the demographic decline experienced by the species over the last millennia. Conclusions Our findings suggest that population connectivity of the species has acted as a buffer of genetic losses and illustrate the use of ancient DNA to uncover such cryptic genetic events.

Genetic diversity and species richness patterns in Baetidae (Ephemeroptera) in the Montseny Mountain range (North-East Iberian Peninsula)

This study aimed to describe patterns of diversity of Baetidae (Ephemeroptera) at the ommunity and population levels within the Montseny Mountain range (North-East Iberian Peninsula). We studied both the distribution of 4 species of baetids in 20 sites among three catchments along the altitudinal gradient (350-1700 masl); and the genetic diversity of the mtDNA cytochrome c oxidase subunit I (cox1) gene of the two common species Baetis alpinus and Baetis rhodani. We found a gradual replacement of the dominant species along the altitudinal gradient. Baetis alpinus inhabited sites at high-altitudes, and this species was replaced by B. rhodani when the altitude decreased. Baetis melanonyx and Alainites muticus attained low abundance at all river sections, and no clear altitudinal trend appeared. Our hypothesis at the population level was that genetic structuring is associated with geographic distance and limited by drainage boundaries among the three studied catchments because of the short-time dispersion of adults. Unexpectedly, analyses of molecular variance (AMOVA) and isolation-bydistance (IBD) showed genetic diversity was unstructured by distance for both species, which may be explained by the relatively short spatial scale studied and small topographic barriers among the three catchments. The Generalized Mixed Yule-Coalescent (GMYC) model showed that B. rhodani had two differentiated genetic lineages that co-occurred in all sites. Overall, diversity of baetids was structured at the community level along the altitudinal gradient, whereas it was unstructured at the population level within the Montseny Mountain range.

Genetic diversity and species richness patterns in Baetidae (Ephemeroptera) in the Montseny Mountain range (North-East Iberian Peninsula)

This study aimed to describe patterns of diversity of Baetidae (Ephemeroptera) at the ommunity and population levels within the Montseny Mountain range (North-East Iberian Peninsula). We studied both the distribution of 4 species of baetids in 20 sites among three catchments along the altitudinal gradient (350-1700 masl); and the genetic diversity of the mtDNA cytochrome c oxidase subunit I (cox1) gene of the two common species Baetis alpinus and Baetis rhodani. We found a gradual replacement of the dominant species along the altitudinal gradient. Baetis alpinus inhabited sites at high-altitudes, and this species was replaced by B. rhodani when the altitude decreased. Baetis melanonyx and Alainites muticus attained low abundance at all river sections, and no clear altitudinal trend appeared. Our hypothesis at the population level was that genetic structuring is associated with geographic distance and limited by drainage boundaries among the three studied catchments because of the short-time dispersion of adults. Unexpectedly, analyses of molecular variance (AMOVA) and isolation-bydistance (IBD) showed genetic diversity was unstructured by distance for both species, which may be explained by the relatively short spatial scale studied and small topographic barriers among the three catchments. The Generalized Mixed Yule-Coalescent (GMYC) model showed that B. rhodani had two differentiated genetic lineages that co-occurred in all sites. Overall, diversity of baetids was structured at the community level along the altitudinal gradient, whereas it was unstructured at the population level within the Montseny Mountain range.

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.

PCR amplification and sequence analyses of reverse transcriptase-like genes in Crinipellis perniciosa isolates

Reverse transcriptase (RT) sequence analysis is an important technique used to detect the presence of transposable elements in a genome. Putative RT sequences were analyzed in the genome of the pathogenic fungus C. perniciosa, the causal agent of witches' broom disease of cocoa. A 394 bp fragment was amplified from genomic DNA of different isolates of C. perniciosa belonging to C-, L-, and S-biotypes and collected from various geographical areas. The cleavage of PCR products with restriction enzymes and the sequencing of various RT fragments indicated the presence of several sequences showing transition events (G:C to A:T). Southern blot analysis revealed high copy numbers of RT signals, forming different patterns among C-, S-, and L-biotype isolates. Sequence comparisons of the predicted RT peptide indicate a close relationship with the RT protein from thegypsy family of LTR-retrotransposons. The possible role of these retrotransposons in generating genetic variability in the homothallic C. perniciosa is discussed.

Genetic variability of Puccinia triticina Eriks. in Brazil

Studies on the genetic variability of Puccinia triticina in inoculum collected in Brazil started in 1941 with Vallega (20). The pioneering work in Brazil dates from 1949 (16) at "Instituto Agronômico do Sul", Ministry of Agriculture (MA), in Pelotas, Rio Grande do Sul State (RS), and continued after 1975 at Embrapa Wheat in Passo Fundo, RS. In 2002, analyses for the identification of P. triticina races continued at OR Seed breeding, simultaneously to Embrapa's program, both in Passo Fundo. The investigators involved in the identification of races in Brazil were Ady Raul da Silva in Pelotas (MA), Eliza Coelho in Pelotas (MA) and in Passo Fundo (Embrapa), Amarilis Labes Barcellos in Pelotas (MA) and in Passo Fundo (Embrapa and OR), Camila Turra in Passo Fundo (OR) and Marcia Chaves in Passo Fundo (Embrapa). From 1979 to 2010 growing season, 59 races were determined, according to the differentiation based on the expression of each Lr resistance gene. On average, one to three new races are detected per year. Research has focused on the use of vertical resistance; however, lately some institutes have searched more durable resistance, of the adult-plant type (horizontal, less race-specific). The uninterrupted monitoring of the wheat rust pathogenic population in Brazil during so many decades allowed the understanding of the evolution and virulence of races. The use of international nomenclature adopted by some programs has allowed the comparison of the fungus variability in Brazil with that in other countries, especially where frontiers are not barriers for spore transportation, confirmed by the occurrence of the same races all over one region.

Genetic structure and diversity of Copaifera langsdorffii Desf. in Cerrado fragments of the São Paulo State, Brazil

The loss of large areas of Cerrado (Brazilian savanna) in Brazil can lead to reduced biodiversity and to the extinction of species. Therefore, the present study aimed to investigate the genetic fragility of populations of Copaifera langsdorffii Desf exposed to different anthropic conditions in fragments of Cerrado in the state of São Paulo. The study was carried out in two Experimental Stations operated by the Forest Institute (Assis and Itirapina), in one fully protected conservation unit (Pedregulho) and in one private property (Brotas). Analyses were conducted using leaf samples from 353 adult specimens and eight pairs of microsatellite loci. The number of alleles per locus ranged from 13 to 15 in all populations, but the mean number of effective alleles was approximately half this value (7.2 to 9-1). Observed heterozygosity was significant and lower than the expected in all populations. Consequently, all populations deviated from Hardy-Weinberg expected frequencies. Fixation indexes were significant for all populations, with the Pedregulho population having the lowest value (0.189) and Itirapina having the highest (0.283). The analysis of spatial genetic structure detected family structures at distance classes of 20 to 65 m in the populations studied. No clones were detected in the populations. Estimates of effective population size were low, but the area occupied by each population studied was large enough for conservation, medium and long term. Recent reductions or bottlenecks were detected in all four populations. Mean Gst’ (genetic divergence) indicated that most of the variation was within populations. Cluster structure analysis based on the genotypes detected K= 4 clusters with distinct allele frequencies patterns. The genetic differentiation observed among populations is consistent with the hypothesis of genetic and geographic isolation. Therefore, it is essential to adopt conservation strategies that raise the gene flow between fragments.

In silico phylogenetic and virulence gene profile analyses of avian pathogenic Escherichia coli genome sequences

Avian pathogenic Escherichia coli (APEC) infections are responsible for significant losses in the poultry industry worldwide. A zoonotic risk has been attributed to APEC strains because they present similarities to extraintestinal pathogenic E. coli (ExPEC) associated with illness in humans, mainly urinary tract infections and neonatal meningitis. Here, we present in silico analyses with pathogenic E. coli genome sequences, including recently available APEC genomes. The phylogenetic tree, based on multi-locus sequence typing (MLST) of seven housekeeping genes, revealed high diversity in the allelic composition. Nevertheless, despite this diversity, the phylogenetic tree was able to cluster the different pathotypes together. An in silico virulence gene profile was also determined for each of these strains, through the presence or absence of 83 well-known virulence genes/traits described in pathogenic E. coli strains. The MLST phylogeny and the virulence gene profiles demonstrated a certain genetic similarity between Brazilian APEC strains, APEC isolated in the United States, UPEC (uropathogenic E. coli) and diarrheagenic strains isolated from humans. This correlation corroborates and reinforces the zoonotic potential hypothesis proposed to APEC.

Holstein white coat color and performance: phenotypic, genetic and environmental correlations

Correlations of measures of percentages of white coat color, five measures of production and two measures of reproduction were obtained from 4293 first lactation Holsteins from eight Florida dairy farms. Percentages of white coat color were analyzed as recorded and transformed by an extension of Box-Cox procedures. Statistical analyses were by derivative-free restricted maximum likelihood (DFREML) with an animal model. Phenotypic and genetic correlations of white percentage (not transformed) were with milk yield, 0.047 and 0.097; fat yield, 0.002 and 0.004; fat percentage, -0.047 and -0.090; protein yield, 0.024 and 0.048; protein percentage, -0.070 and -0.116; days open, -0.012 and -0.065; and calving interval, -0.007 and -0.029. Changes in magnitude of correlations were very small for all variables except days open. Genetic and phenotypic correlations of transformed values with days open were -0.027 and -0.140. Modest positive correlated responses would be expected for white coat color percentage following direct selection for milk, fat, and protein yields, but selection for fat and protein percentages, days open, or calving interval would lead to small decreases.