Sex-linked genetic markers in the shrews of the Sorex araneus group : from free gene flow to speciation

SUMMARY: The shrews of the Sorex araneus group are morphologically .very similar, but have undergone a spectacular chromosomal evolution. Altogether, the shrews of this group present a complete array of every possible level of chromosomal and genetic differentiation. In South-Western Europe, four species are recognised: S. antiriorii, S. araneus, S. coronatus and S. granarius, which differ essentially by the amount and the composition of Robertsonian metacentric chromosomés. Additionally, several chromosome races of S. araneus are also present in the same region (i.e. Bretolet, Carlit, Cordon, Jura and Vaud). The objective of this thesis was to examine the genetic relationships between populations, races and /or species of the Sorex araneus group with a special emphasis onsex-specific markers (mtDNA and Y chromosome). We first investigate the evolutionary history of the shrews of the Sorex araneus group distributed in the South-Western Europe. The results of. these analyses confirmed the difficulty to draw a single dichotomic tree within this group. Incongruent mtDNA and Y chromosome phylogenies suggest further that genetic and chromosomal evolution are in this group partially independent processes and that the evolutionary history of the south-western European populations of the S. araneus group can only be understood if we consider secondary contacts between taxa, after their divergence (with genetic exchanges by means of hybridization and / or introgression). Using one male-inherited, one female inherited and eight biparentally inherited markers, we investigate the population genetic structure of the Valais shrew (Sorex antinorii). Overall there results suggest that two already well-differentiated genetic lineages colonized the Swiss Alps after the last glacial period and came into contact in the Rhône Valley. After the Valais shrew (Sorex antinorii) reached the Swiss Alps, it came into contact with the common shrew (Sorex araneus). When two species come into contact and hybridize, endogenous counter-selection of hybrids is usually first expressed as a reduced fertility or viability in hybrids of the heterogametic sex, a mechanism know as Haldane's rule (Haldane 1922). We first evaluated the extent of introgression for Y chromosome, mtDNA and autosomal markers in a hybrid zone between S. antinoriii and S. araneus. The overall level of genetic and karyotypic differentiation between the two species must be strong .enough to allow the detection asymmetric introgression. Secondly, we compared the levels of gene flow between chromosome common to both species and chromosome differently rearranged in each of them. We detected a significantly stronger genetic structure in rearranged chromosomes. Over a 10-year period, we even observed a decrease of genetic structure for common chromosomes. These results strongly support the role of chromosomal rearrangements in the reproductive barrier between S. araneus and S. anfinorii. Overall, this thesis underlines the need to use different inherited (paternally, maternally and / or biparentally) and chromosomally located (on common vs. on rearranged chromosomes) markers to obtain more accurate pictures of genetic relationships between populations or species. RÉSUMÉ: Les musaraignes du groupe Sorex araneus sont morphologiquement très proches, mais ont connu une spectaculaire évolution chromosomique. Prises dans leur ensemble, les musaraignes de ce groupe présentent tous les nivaux possibles de différenciation génétique et chromosomique. Dans le sud-ouest de l'Europe, quatre espèces appartenant à ce groupe sont présentes : S. antinorii, S. araneus, S. coronatus et S. granarius. Celles-ci diffèrent essentiellement par leur caryotype dont la variabilité est principalement due à des fusions Robertsoniennes. De plus, plusieurs races chromosomiques appartenant à S. araneus sont aussi présentes dans la même région (i.e. les races Bretolet, Carlit, Cordon, Jura et Vaud). L'objectif de cette thèse était d'examiner les relations génétiques entre populations, races et/ou espèces du groupe S. araneus, en utilisant particulièrement des marqueurs liés aux sexes (ADN mitochondrial et Chromosome Y). Nous avons dans un premier temps retracé l'histoire évolutive des musaraignes de ce groupe dans le sud-ouest de l'Europe. Les résultats dé ces analyses confirment qu'il est difficile de tracer un simple arbre dichotomique au sein de ce groupe. Les arbres phylogénétiques obtenus sur l'ADN mitochondrial et le chromosome Y sont incongruents et suggèrent de plus que l'évolution génétique et chromosomique sont des processus indépendants. L'histoire évolutive -des populations de ce groupe ne peut. être comprise qu'en considérant des contacts secondaires entre taxa postérieure à leur divergence et induisant des échanges génétiques par hybridation et/ou introgression. Par la suite, nous avons examiné la structure génétique des populations de la musaraigne du Valais, S. antinorii, en utilisant un marqueur transmis par les mâles, un marqueur transmis par les femelles et huit marqueurs transmis par les 2 sexes. Nos résultats suggèrent que deux lignées génétiquement bien différenciées aient colonisé les Alpes Suisses, après les dernières glaciations et entrent en contact dans là Vallée du Rhône. Après avoir franchi les Alpes Suisses, la musaraigne du Valais est entrée en contact avec là musaraigne commune (S. araneus). Lorsque deux espèces entrent en contact et s'hybrident, la sélection contre les hybrides implique habituellement une baisse de fertilité ou de viabilité des hybrides du sexe hétérogamétique (i.e. les mâles XY chez les mammifères). Ce mécanisme est connu sous le nom de règle de Haldane (Haldane 1922) et implique une plus forte structuration génétique de marqueurs males - spécifiques que des marqueurs femelles spécifiques. Nous avons donc évalué le degré d'introgression des marqueurs situés sur le chromosome Y, sur l'ADN mitochondrial et sur des autosomes dans une zone hybride entre S. araneus et S. antinorii. Le niveau de différenciation chromosomique et génétique entre les 2 espèces doit être suffisamment fort pour ne pas permettre la détection d'une introgression asymétrique entre les sexes. Dans un second temps, nous avons comparé les niveaux de flux de gênes mesurés à l'échelle du chromosome, pour des chromosomes communs aux deux espèces et pour des chromosomes différemment arrangées dans chacune des deux espèces. Nous avons détecté une structure génétique significativement plus forte sur les chromosomes réarrangés et comme la zone hybride a été étudiée à dix années d'intervalle, nous observons même une diminution de la structure génétique pour les chromosomes communs au cours du temps.. Ces résultats soutiennent fortement l'hypothèse d'un rôle des réarrangements chromosomiques dans l'établissement d'une barrière reproductive entre S. araneus et S. antinorii. Ainsi cette thèse souligne l'utilité d'utiliser des marqueurs génétiques avec différents modes de transmission. (par les mâles, par les femelles et/ou par les 2 sexes) ou localisés au niveau du chromosome (chromosomes communs vs chromosomes réarrangés) afin d'obtenir une image plus juste ou du moins plus complète des relations génétiques entre populations ou espèces.

Genetic isolation of insular populations of the Maghrebian bat, Myotis punicus, in the Mediterranean Basin

We investigate the population genetic structure of the Maghrebian bat, Myotis punicus, between the mainland and islands to assess the island colonization pattern and current gene flow between nearby islands and within the mainland. Location North Africa and the Mediterranean islands of Corsica and Sardinia. Methods We sequenced part of the control region (HVII) of 79 bats across 11 colonies. The phylogeographical pattern was assessed by analysing molecular diversity indices, examining differentiation among populations and estimating divergence time. In addition, we genotyped 182 bats across 10 colonies at seven microsatellite loci. We used analysis of molecular variance and a Bayesian approach to infer nuclear population structure. Finally, we estimated sex-specific dispersal between Corsica and Sardinia. Results Mitochondrial analyses indicated that colonies between Corsica, Sardinia and North Africa are highly differentiated. Within islands there was no difference between colonies, while at the continental level Moroccan and Tunisian populations were highly differentiated. Analyses with seven microsatellite loci showed a similar pattern. The sole difference was the lack of nuclear differentiation between populations in North Africa, suggesting a male-biased dispersal over the continental area. The divergence time of Sardinian and Corsican populations was estimated to date back to the early and mid-Pleistocene. Main conclusions Island colonization by the Maghrebian bats seems to have occurred in a stepping-stone manner and certainly pre-dated human colonization. Currently, open water seems to prevent exchange of bats between the two islands, despite their ability to fly and the narrowness of the strait of Bonifacio. Corsican and Sardinian populations are thus currently isolated from any continental gene pool and must therefore be considered as different evolutionarily significant units (ESU).

Scale-specific sex-biased dispersal in the Valais shrew unveiled by genetic variation on the Y chromosome, autosomes, and mitochondrial DNA.

We investigated sex specificities in the evolutionary processes shaping Y chromosome, autosomes, and mitochondrial DNA patterns of genetic structure in the Valais shrew (Sorex antinorii), a mountain dwelling species with a hierarchical distribution. Both hierarchical analyses of variance and isolation-by-distance analyses revealed patterns of population structure that were not consistent across maternal, paternal, and biparentally inherited markers. Differentiation on a Y microsatellite was lower than expected from the comparison with autosomal microsatellites and mtDNA, and it was mostly due to genetic variance among populations within valleys, whereas the opposite was observed on other markers. In addition, there was no pattern of isolation by distance for the Y, whereas there was strong isolation by distance on mtDNA and autosomes. We use a hierarchical island model of coancestry dynamics to discuss the relative roles of the microevolutionary forces that may induce such patterns. We conclude that sex-biased dispersal is the most important driver of the observed genetic structure, but with an intriguing twist: it seems that dispersal is strongly male biased at large spatial scale, whereas it is mildly biased in favor of females at local scale. These results add to recent reports of scale-specific sex-biased dispersal patterns, and emphasize the usefulness of the Y chromosome in conjunction with mtDNA and autosomes to infer sex specificities.

Genetic Association of Recovery from Eating Disorders: The Role of GABA Receptor SNPs.

Follow-up studies of eating disorders (EDs) suggest outcomes ranging from recovery to chronic illness or death, but predictors of outcome have not been consistently identified. We tested 5151 single-nucleotide polymorphisms (SNPs) in approximately 350 candidate genes for association with recovery from ED in 1878 women. Initial analyses focused on a strictly defined discovery cohort of women who were over age 25 years, carried a lifetime diagnosis of an ED, and for whom data were available regarding the presence (n=361 ongoing symptoms in the past year, ie, 'ill') or absence (n=115 no symptoms in the past year, ie, 'recovered') of ED symptoms. An intronic SNP (rs17536211) in GABRG1 showed the strongest statistical evidence of association (p=4.63 × 10(-6), false discovery rate (FDR)=0.021, odds ratio (OR)=0.46). We replicated these findings in a more liberally defined cohort of women age 25 years or younger (n=464 ill, n=107 recovered; p=0.0336, OR=0.68; combined sample p=4.57 × 10(-6), FDR=0.0049, OR=0.55). Enrichment analyses revealed that GABA (γ-aminobutyric acid) SNPs were over-represented among SNPs associated at p<0.05 in both the discovery (Z=3.64, p=0.0003) and combined cohorts (Z=2.07, p=0.0388). In follow-up phenomic association analyses with a third independent cohort (n=154 ED cases, n=677 controls), rs17536211 was associated with trait anxiety (p=0.049), suggesting a possible mechanism through which this variant may influence ED outcome. These findings could provide new insights into the development of more effective interventions for the most treatment-resistant patients.

Clear detection of ADIPOQ locus as the major gene for plasma adiponectin: results of genome-wide association analyses including 4659 European individuals.

OBJECTIVE: Plasma adiponectin is strongly associated with various components of metabolic syndrome, type 2 diabetes and cardiovascular outcomes. Concentrations are highly heritable and differ between men and women. We therefore aimed to investigate the genetics of plasma adiponectin in men and women. METHODS: We combined genome-wide association scans of three population-based studies including 4659 persons. For the replication stage in 13795 subjects, we selected the 20 top signals of the combined analysis, as well as the 10 top signals with p-values less than 1.0 x 10(-4) for each the men- and the women-specific analyses. We further selected 73 SNPs that were consistently associated with metabolic syndrome parameters in previous genome-wide association studies to check for their association with plasma adiponectin. RESULTS: The ADIPOQ locus showed genome-wide significant p-values in the combined (p=4.3 x 10(-24)) as well as in both women- and men-specific analyses (p=8.7 x 10(-17) and p=2.5 x 10(-11), respectively). None of the other 39 top signal SNPs showed evidence for association in the replication analysis. None of 73 SNPs from metabolic syndrome loci exhibited association with plasma adiponectin (p>0.01). CONCLUSIONS: We demonstrated the ADIPOQ gene as the only major gene for plasma adiponectin, which explains 6.7% of the phenotypic variance. We further found that neither this gene nor any of the metabolic syndrome loci explained the sex differences observed for plasma adiponectin. Larger studies are needed to identify more moderate genetic determinants of plasma adiponectin.

Genetic polymorphisms of the main transcription factors for adiponectin gene promoter in regulation of adiponectin levels: association analysis in three European cohorts.

Adiponectin serum concentrations are an important biomarker in cardiovascular epidemiology with heritability etimates of 30-70%. However, known genetic variants in the adiponectin gene locus (ADIPOQ) account for only 2%-8% of its variance. As transcription factors are thought to play an under-acknowledged role in carrying functional variants, we hypothesized that genetic polymorphisms in genes coding for the main transcription factors for the ADIPOQ promoter influence adiponectin levels. Single nucleotide polymorphisms (SNPs) at these genes were selected based on the haplotype block structure and previously published evidence to be associated with adiponectin levels. We performed association analyses of the 24 selected SNPs at forkhead box O1 (FOXO1), sterol-regulatory-element-binding transcription factor 1 (SREBF1), sirtuin 1 (SIRT1), peroxisome-proliferator-activated receptor gamma (PPARG) and transcription factor activating enhancer binding protein 2 beta (TFAP2B) gene loci with adiponectin levels in three different European cohorts: SAPHIR (n = 1742), KORA F3 (n = 1636) and CoLaus (n = 5355). In each study population, the association of SNPs with adiponectin levels on log-scale was tested using linear regression adjusted for age, sex and body mass index, applying both an additive and a recessive genetic model. A pooled effect size was obtained by meta-analysis assuming a fixed effects model. We applied a significance threshold of 0.0033 accounting for the multiple testing situation. A significant association was only found for variants within SREBF1 applying an additive genetic model (smallest p-value for rs1889018 on log(adiponectin) = 0.002, β on original scale = -0.217 µg/ml), explaining ∼0.4% of variation of adiponectin levels. Recessive genetic models or haplotype analyses of the FOXO1, SREBF1, SIRT1, TFAPB2B genes or sex-stratified analyses did not reveal additional information on the regulation of adiponectin levels. The role of genetic variations at the SREBF1 gene in regulating adiponectin needs further investigation by functional studies.

Taxonomic status of Hylomys parvus and Hylomys suillus (Insectivora: Erinaceidae): biochemical and morphological analyses

The genus Hylomys was thought to be represented by a single widespread species. Biochemical and morphometric analyses of several Southeast Asian populations reveal that Sumatra is inhabited by two distinct species, the dwarf gymnure (H. parvus) and the lesser gymnure (H. suillus). The absence of interbreeding between these two groups along with their relatively ancient common origins are documented by several diagnostic loci and a large Nei's genetic distance (D = 0.353 +/- 0.035). The dwarf gymnure has been reported only from the slopes of the Mt. Kerinci volcano in Sumatra, where the species lives at higher elevations than its potential competitor, the lesser gymnure. Other populations of Hylomys from Java, Borneo, and Malaysia are more closely related to the Sumatran sample of H. suillus, but they exhibit strong interpopulational genetic differentiation (D = 0.165 +/- 0.040) that may be accounted for by their isolated montane habitat. In addition, a principal-components analysis based on 16 measurements of the skull clearly separates adult specimens of both species. There is little overlap in the measurements between H. suillus (which is larger) and H. parvus. On Sumatra where both species may be sympatric, the notched space between premaxillary tips, soft texture of the fur, and more delicate skull and dentition are diagnostic of H. parvus.

Genetical and environmental analyses of yield in six biparental soybean crosses

The yearly genetic progress obtained by breeding for increased soybean yield has been considered acceptable worldwide. It is common sense, however, that this progress can be improved further if refined breeding techniques, developed from the knowledge of the genetic mechanisms controlling soybean yield, are used. In this paper, data from four cultivars and/or lines and their derived sets of F2, F3, F7, F8, F9 and F10 generations assayed in 17 environments were analyzed to allow an insight of the genetic control of soybean yield under different environmental conditions. The general picture was of a complex polygene system controlling yield in soybeans. Additive genetic effects predominated although dominance was often found to be significant. Complications such as epistasis, linkage and macro and micro genotype x environment (G x E) interactions were also commonly detected. The overall heritability was 0.29. The relative magnitude of the additive effects and the complicating factors allowed the inference that the latter are not a serious problem to the breeder. The low heritability values and the considerable magnitude of G x E interactions for yield, however, indicated that careful evaluation through experiments designed to allow for the presence of these effects is necessary for successful selection.

Characterization of mandarin citrus germplasm from Southern Brazil by morphological and molecular analyses

The objective of this work was to characterize mandarin (Citrus spp.) germplasm from Southern Brazil by morphological and molecular analyses. Thirty seven cultivars from 34 distinct mandarin varieties were evaluated by morphological and agronomic traits of leaves, flowers and fruits, and by microsatellite markers. The morphological and agronomic characteristics suggested that almost all varieties can be produced for commercial use, and some, as the Satsuma variety, are recommended for breeding programs. Pooled DNA samples from 1-5 plants belonging to each cultivar were tested. Eight of the nine primers detected polymorphisms. Specific markers were found for some accessions. The dendrogram constructed with the morphological results divided the 37 cultivars into four groups, while that obtained with the microsatellites clustered 35 of the 37 cultivars into three groups only. Generally, intervarietal differences are not high, and this lack of agreement in the two multifactorial analyses indicates that diverse evolutionary factors are acting at these two levels of investigation.

Genetic variability of the Brazilian hair sheep breeds

The objectives of this work were to investigate the genetic structure of the Brazilian hair sheep breeds and to determine the origin of the Santa Inês breed. Molecular similarity was determined using Randomly Amplified Polymorphic DNA - Polymerase Chain Reaction markers in 238 individuals from five naturalized sheep breeds: Santa Inês (48 animals), Rabo Largo (48), Somali (48), Morada Nova (48) and Bergamasca (46), collected in Goiás, Sergipe, Bahia, and Ceará States as well as in the Federal District. Fifty-four loci were selected from 19 primers, after a pilot test using 140 primers. Qualitative analyses indicate diagnostic markers for all breeds. All breeds were significantly different from each other. Interbreed differences were explained by 14.92% of the total variation. Santa Inês clustered with Bergamasca (97% bootstrap) and with Rabo Largo, composing the third member of the group (81% bootstrap) while Morada Nova and Somali breeds clustered separately. Each breed should be considered as a separate management and conservation unit, and special care should be taken with Rabo Largo, Morada Nova and Somali breeds, represented by small herds in Brazil.

Genetic abnormalities and survival in multiple myeloma: the experience of the Intergroupe Francophone du Myélome.

Acquired genomic aberrations have been shown to significantly impact survival in several hematologic malignancies. We analyzed the prognostic value of the most frequent chromosomal changes in a large series of patients with newly diagnosed symptomatic myeloma prospectively enrolled in homogeneous therapeutic trials. All the 1064 patients enrolled in the IFM99 trials conducted by the Intergroupe Francophone du Myélome benefited from an interphase fluorescence in situ hybridization analysis performed on purified bone marrow plasma cells. They were systematically screened for the following genomic aberrations: del(13), t(11;14), t(4;14), hyperdiploidy, MYC translocations, and del(17p). Chromosomal changes were observed in 90% of the patients. The del(13), t(11;14), t(4;14), hyperdiploidy, MYC translocations, and del(17p) were present in 48%, 21%, 14%, 39%, 13%, and 11% of the patients, respectively. After a median follow-up of 41 months, univariate statistical analyses revealed that del(13), t(4;14), nonhyperdiploidy, and del(17p) negatively impacted both the event-free survival and the overall survival, whereas t(11;14) and MYC translocations did not influence the prognosis. Multivariate analyses on 513 patients annotated for all the parameters showed that only t(4;14) and del(17p) retained prognostic value for both the event-free and overall survivals. When compared with the currently used International Staging System, this prognostic model compares favorably. In myeloma, the genomic aberrations t(4;14) and del(17p), together with beta2-microglobulin level, are important independent predictors of survival. These findings have implications for the design of risk-adapted treatment strategies.

Genetic distance estimated by RAPD markers and its relationship with hybrid performance in maize

The objective of this work was to evaluate the genetic diversity of 16 maize inbred lines, and to determine the correlation between genetic distance and hybrid performance, using random amplified polymorphic DNA (RAPD) molecular markers. Twenty-two different random primers were used, which resulted in the amplification of 265 fragments, 237 (84.44%) of them being polymorphic. A genetic similarity matrix was created from the RAPD data, using Jaccard coefficient, and a dendrogram was constructed. Hybrid analyses were carried out using random block design and Griffing method VI for diallel crossings. The genetic associations showed five distinct heterotic groups. Correlations between genetic divergences detected by RAPD, as well as the means observed in the diallel crossings were positive and significant for plant height, ear height, prolificacy, and grain weight. The correlation of genetic divergences, detected by RAPD, and the specific combining ability between heterotic group associations, showed significance in all characteristics under study, except prolificacy. A direct relationship between genetic divergence and productivity was found in 79.2% of the 120 hybrids confirming the hypothesis that genetic divergence is directly related to the performance of hybrids and is efficient in predicting it.

GWAS of 126,559 individuals identifies genetic variants associated with educational attainment.

A genome-wide association study (GWAS) of educational attainment was conducted in a discovery sample of 101,069 individuals and a replication sample of 25,490. Three independent single-nucleotide polymorphisms (SNPs) are genome-wide significant (rs9320913, rs11584700, rs4851266), and all three replicate. Estimated effects sizes are small (coefficient of determination R(2) ≈ 0.02%), approximately 1 month of schooling per allele. A linear polygenic score from all measured SNPs accounts for ≈2% of the variance in both educational attainment and cognitive function. Genes in the region of the loci have previously been associated with health, cognitive, and central nervous system phenotypes, and bioinformatics analyses suggest the involvement of the anterior caudate nucleus. These findings provide promising candidate SNPs for follow-up work, and our effect size estimates can anchor power analyses in social-science genetics.

Genetic diversity in caribou linked to past and future climate change

Climate-driven range fluctuations during the Pleistocene have continuously reshaped species distribution leading to populations of contrasting genetic diversity. Contemporary climate change is similarly influencing species distribution and population structure, with important consequences for patterns of genetic diversity and species' evolutionary potential1. Yet few studies assess the impacts of global climatic changes on intraspecific genetic variation2, 3, 4, 5. Here, combining analyses of molecular data with time series of predicted species distributions and a model of diffusion through time over the past 21 kyr, we unravel caribou response to past and future climate changes across its entire Holarctic distribution. We found that genetic diversity is geographically structured with two main caribou lineages, one originating from and confined to Northeastern America, the other originating from Euro-Beringia but also currently distributed in western North America. Regions that remained climatically stable over the past 21 kyr maintained a high genetic diversity and are also predicted to experience higher climatic stability under future climate change scenarios. Our interdisciplinary approach, combining genetic data and spatial analyses of climatic stability (applicable to virtually any taxon), represents a significant advance in inferring how climate shapes genetic diversity and impacts genetic structure.

Multivariate analyses of genotype x environment interaction of popcorn

The objectives of this work were to evaluate the genotype x environment (GxE) interaction for popcorn and to compare two multivariate analyses methods. Nine popcorn cultivars were sown on four dates one month apart during each of the agricultural years 1998/1999 and 1999/2000. The experiments were carried out using randomized block designs, with four replicates. The cv. Zélia contributed the least to the GxE interaction. The cv. Viçosa performed similarly to cv. Rosa-claro. Optimization of GxE was obtained for cv. CMS 42 for a favorable mega-environment, and for cv. CMS 43 for an unfavorable environment. Multivariate analysis supported the results from the method of Eberhart & Russell. The graphic analysis of the Additive Main effects and Multiplicative Interaction (AMMI) model was simple, allowing conclusions to be made about stability, genotypic performance, genetic divergence between cultivars, and the environments that optimize cultivar performance. The graphic analysis of the Genotype main effects and Genotype x Environment interaction (GGE) method added to AMMI information on environmental stratification, defining mega-environments and the cultivars that optimized performance in those mega-environments. Both methods are adequate to explain the genotype x environment interactions.