Functionally associated molecular genetic marker map construction in perennial ryegrass (Lolium perenne L.)

A molecular marker-based map of perennial ryegrass (Lolium perenne L.) has been constructed through the use of polymorphisms associated with expressed sequence tags (ESTs). A pair-cross between genotypes from a North African ecotype and the cultivar Aurora was used to generate a two-way pseudo-testcross population. A selection of 157 cDNAs assigned to eight different functional categories associated with agronomically important biological processes was used to detect polymorphic EST–RFLP loci in the F1(NA6 × AU6) population. A comprehensive set of EST–SSR markers was developed from the analysis of 14,767 unigenes, with 310 primer pairs showing efficient amplification and detecting 113 polymorphic loci. Two parental genetic maps were produced: the NA6 genetic map contains 88 EST–RFLP and 71 EST–SSR loci with a total map length of 963 cM, while the AU6 genetic map contains 67 EST–RFLP and 58 EST–SSR loci with a total map length of 757 cM. Bridging loci permitted the alignment of homologous chromosomes between the parental maps, and a sub-set of genomic DNA-derived SSRs was used to relate linkage groups to the perennial ryegrass reference map. Regions of segregation distortion were identified, in some instances in common with other perennial ryegrass maps. The EST-derived marker-based map provides the basis for in silico comparative genetic mapping, as well as the evaluation of co-location between QTLs and functionally associated genetic loci.

Construction of microsatellite linkage maps for Corymbia

The genus Corymbia is closely related to the genus Eucalyptus, and like Eucalyptus contains tree species that are important for sub-tropical forestry. Corymbia's close relationship with Eucalyptus suggests genetic studies in Corymbia should benefit from transfer of genetic information from its more intensively studied relatives. Here we report a genetic map for Corymbia spp. based on microsatellite markers identified de novo in Corymbia sp or transferred from Eucalyptus. A framework consensus map was generated from an outbred F 2 population (n = 90) created by crossing two unrelated Corymbia torelliana x C. citriodora subsp. variegata F1 trees. The map had a total length of 367 cM (Kosambi) and was composed of 46 microsatellite markers distributed across 13 linkage groups (LOD 3). A high proportion of Eucalyptus microsatellites (90%) transferred to Corymbia. Comparative analysis between the Corymbia map and a published Eucalyptus map identified eight homeologous linkage groups in Corymbia with 13 markers mapping on one or both maps. Further comparative analysis was limited by low power to detect linkage due to low genome coverage in Corymbia, however, there was no convincing evidence for chromosomal structural differences because instances of non-synteny were associated with large distances on the Eucalyptus map. Segregation distortion was primarily restricted to a single linkage group and due to a deficit of hybrid genotypes, suggesting that hybrid inviability was one factor shaping the genetic composition of the F2 population in this inter-subgeneric hybrid. The conservation of microsatellite loci and synteny between Corymbia and Eucalyptus suggests there will be substantial value in exchanging information between the two groups.

A consensus genetic map of sorghum that integrates multiple component maps and high-throughput Diversity Array Technology (DArT) markers

Background: Sorghum genome mapping based on DNA markers began in the early 1990s and numerous genetic linkage maps of sorghum have been published in the last decade, based initially on RFLP markers with more recent maps including AFLPs and SSRs and very recently, Diversity Array Technology (DArT) markers. It is essential to integrate the rapidly growing body of genetic linkage data produced through DArT with the multiple genetic linkage maps for sorghum generated through other marker technologies. Here, we report on the colinearity of six independent sorghum component maps and on the integration of these component maps into a single reference resource that contains commonly utilized SSRs, AFLPs, and high-throughput DArT markers. Results: The six component maps were constructed using the MultiPoint software. The lengths of the resulting maps varied between 910 and 1528 cM. The order of the 498 markers that segregated in more than one population was highly consistent between the six individual mapping data sets. The framework consensus map was constructed using a "Neighbours" approach and contained 251 integrated bridge markers on the 10 sorghum chromosomes spanning 1355.4 cM with an average density of one marker every 5.4 cM, and were used for the projection of the remaining markers. In total, the sorghum consensus map consisted of a total of 1997 markers mapped to 2029 unique loci ( 1190 DArT loci and 839 other loci) spanning 1603.5 cM and with an average marker density of 1 marker/0.79 cM. In addition, 35 multicopy markers were identified. On average, each chromosome on the consensus map contained 203 markers of which 58.6% were DArT markers. Non-random patterns of DNA marker distribution were observed, with some clear marker-dense regions and some marker-rare regions. Conclusion: The final consensus map has allowed us to map a larger number of markers than possible in any individual map, to obtain a more complete coverage of the sorghum genome and to fill a number of gaps on individual maps. In addition to overall general consistency of marker order across individual component maps, good agreement in overall distances between common marker pairs across the component maps used in this study was determined, using a difference ratio calculation. The obtained consensus map can be used as a reference resource for genetic studies in different genetic backgrounds, in addition to providing a framework for transferring genetic information between different marker technologies and for integrating DArT markers with other genomic resources. DArT markers represent an affordable, high throughput marker system with great utility in molecular breeding programs, especially in crops such as sorghum where SNP arrays are not publicly available.

AFLP-based genetic linkage maps of the Pacific oyster Crassostrea gigas Thunberg

Amplified fragment length polymorphisms (AFLPs) were used for genome mapping in the Pacific Oyster Crassostrea gigas Thunberg. Seventeen selected primer combinations produced 1106 peaks, of which 384 (34.7%) were polymorphic in a backcross family. Among the polymorphic markers, 349 were segregating through either the female or the male parent. Chi-square analysis indicated that 255 (73.1%) of the markers segregated in a Mendelian ratio, and 94 (26.9%) showed significant (P < 0.05) segregation distortion. Separate genetic linkage maps were constructed for the female and male parents. The female framework map consisted of 119 markers in 11 linkage groups, spanning 1030.7 cM, with an average interval of 9.5 cM per marker. The male map contained 96 markers in 10 linkage groups, covering 758.4 cM, with 8.8 cM per marker. The estimated genome length of the Pacific oyster was 1258 cM for the female and 933 cM for the male, and the observed coverage was 82.0% for the female map and 81.3% for the male map. Most distorted markers were deficient for homozygotes and closely linked to each other on the genetic map, suggesting the presence of major recessive deleterious genes in the Pacific oyster.

Genome-wide linkage analysis for human longevity: Genetics of healthy aging study

Clear evidence exists for heritability of humanlongevity, and much interest is focused on identifying genes associated with longer lives. To identify such longevity alleles, we performed the largest genome-wide linkage scan thus far reported. Linkage analyses included 2118nonagenarian Caucasian sibling pairs that have been enrolled in 15 study centers of 11 European countries as part of the Genetics of Healthy Aging (GEHA) project. In the joint linkage analyses, we observed four regions that show linkage with longevity; chromosome 14q11.2 (LOD = 3.47), chromosome 17q12-q22 (LOD = 2.95), chromosome 19p13.3-p13.11 (LOD = 3.76), and chromosome 19q13.11-q13.32 (LOD = 3.57). To fine map these regions linked to longevity, we performed association analysis using GWAS data in a subgroup of 1228 unrelated nonagenarian and 1907 geographically matched controls. Using a fixed-effect meta-analysis approach, rs4420638 at the TOMM40/ APOE/APOC1 gene locus showed significant association with longevity (P-value = 9.6 × 10). By combined modeling of linkage and association, we showed that association of longevity with APOEe4 and APOEe2 alleles explain the linkage at 19q13.11-q13.32 with P-value = 0.02 and P-value = 1.0 × 10, respectively. In the largest linkage scan thus far performed for human familial longevity, we confirm that the APOE locus is a longevity gene and that additional longevity loci may be identified at 14q11.2, 17q12-q22, and 19p13.3-p13.11. As the latter linkage results are not explained by common variants, we suggest that rare variants play an important role in human familial longevity.

Physical map and cosmid contig encompassing a new interstitial deletion of the X-linked lymphoproliferative syndrome region

The X-linked lymphoproliferative syndrome (XLP) is an inherited immuno-deficiency to Epstein-Barr virus infection that has been mapped to chromosome Xq25. Molecular analysis of XLP patients from ten different families identified a small interstitial constitutional deletion in 1 patient (XLP-D). This deletion, initially defined by a single marker, DF83, known to map to interval Xq24-q26.1, is nested within a previously reported and much larger deletion in another XLP patient (XLP-739). A cosmid minilibrary was constructed from a single mega-YAC and used to establish a contig encompassing the whole XLP-D deletion and a portion of the XLP-739 deletion. Based on this contig, the size of the XLP-D deletion can be estimated at 130 kb. The identification of this minimal deletion, within which at least a portion of the XLP gene is likely to reside, should greatly facilitate efforts in isolating the gene.

Trans-Ancestral Studies Fine Map the SLE-Susceptibility Locus TNFSF4

We previously established an 80 kb haplotype upstream of TNFSF4 as a susceptibility locus in the autoimmune disease SLE. SLE-associated alleles at this locus are associated with inflammatory disorders, including atherosclerosis and ischaemic stroke. In Europeans, the TNFSF4 causal variants have remained elusive due to strong linkage disequilibrium exhibited by alleles spanning the region. Using a trans-ancestral approach to fine-map the locus, utilising 17,900 SLE and control subjects including Amerindian/Hispanics (1348 cases, 717 controls), African-Americans (AA) (1529, 2048) and better powered cohorts of Europeans and East Asians, we find strong association of risk alleles in all ethnicities; the AA association replicates in African-American Gullah (152,122). The best evidence of association comes from two adjacent markers: rs2205960-T (P = 1.71×10-34, OR = 1.43[1.26-1.60]) and rs1234317-T (P = 1.16×10-28, OR = 1.38[1.24-1.54]). Inference of fine-scale recombination rates for all populations tested finds the 80 kb risk and non-risk haplotypes in all except African-Americans. In this population the decay of recombination equates to an 11 kb risk haplotype, anchored in the 5′ region proximal to TNFSF4 and tagged by rs2205960-T after 1000 Genomes phase 1 (v3) imputation. Conditional regression analyses delineate the 5′ risk signal to rs2205960-T and the independent non-risk signal to rs1234314-C. Our case-only and SLE-control cohorts demonstrate robust association of rs2205960-T with autoantibody production. The rs2205960-T is predicted to form part of a decameric motif which binds NF-κBp65 with increased affinity compared to rs2205960-G. ChIP-seq data also indicate NF-κB interaction with the DNA sequence at this position in LCL cells. Our research suggests association of rs2205960-T with SLE across multiple groups and an independent non-risk signal at rs1234314-C. rs2205960-T is associated with autoantibody production and lymphopenia. Our data confirm a global signal at TNFSF4 and a role for the expressed product at multiple stages of lymphocyte dysregulation during SLE pathogenesis. We confirm the validity of trans-ancestral mapping in a complex trait. © 2013 Manku et al.

Genetic structure and linkage disequilibrium in landrace populations of barley in Sardinia

Multilocus digenic linkage disequilibria (LD) and their population structure were investigated in eleven landrace populations of barley (Hordeum vulgare ssp. vulgare L.) in Sardinia, using 134 dominant simple-sequence amplified polymorphism markers. The analysis of molecular variance for these markers indicated that the populations were partially differentiated (F ST = 0.18), and clustered into three geographic areas. Consistent with this population pattern, STRUCTURE analysis allocated individuals from a bulk of all populations into four genetic groups, and these groups also showed geographic patterns. In agreement with other molecular studies in barley, the general level of LD was low (13 % of locus pairs, with P < 0.01) in the bulk of 337 lines, and decayed steeply with map distance between markers. The partitioning of multilocus associations into various components indicated that genetic drift and founder effects played a major role in determining the overall genetic makeup of the diversity in these landrace populations, but that epistatic homogenising or diversifying selection was also present. Notably, the variance of the disequilibrium component was relatively high, which implies caution in the pooling of barley lines for association studies. Finally, we compared the analyses of multilocus structure in barley landrace populations with parallel analyses in both composite crosses of barley on the one hand and in natural populations of wild barley on the other. Neither of these serves as suitable mimics of landraces in barley, which require their own study. Overall, the results suggest that these populations can be exploited for LD mapping if population structure is controlled.

A SNP-based consensus genetic map for synteny-based trait targeting in faba bean (Vicia faba L.)

Faba bean (Vicia faba L.) is a globally important nitrogen-fixing legume, which is widely grown in a diverse range of environments. In this work, we mine and validate a set of 845 SNPs from the aligned transcriptomes of two contrasting inbred lines. Each V. faba SNP is assigned by BLAST analysis to a single Medicago orthologue. This set of syntenically anchored polymorphisms were then validated as individual KASP assays, classified according to their informativeness and performance on a panel of 37 inbred lines, and the best performing 757 markers used to genotype six mapping populations. The six resulting linkage maps were merged into a single consensus map on which 687 SNPs were placed on six linkage groups, each presumed to correspond to one of the six V. faba chromosomes. This sequence-based consensus map was used to explore synteny with the most closely-related crop species, lentil, and the most closely related fully sequenced genome, Medicago. Large tracts of uninterrupted colinearity were found between faba bean and Medicago, making it relatively straightforward to predict gene content and order in mapped genetic interval. As a demonstration of this, we mapped a flower colour gene to a 2 cM interval of Vf chromosome 2 which was highly collinear with Mt3. The obvious candidate gene from 77 gene models in the collinear Medicago chromosome segment was the previously characterized MtWD40-1 gene (Mt3g092830, Mt3g092840) controlling anthocyanin production in Medicago and re-sequencing of the Vf orthologue showed a putative causative deletion of the entire 5’ end of the gene.

A radiation hybrid map of river buffalo (Bubalus bubalis) chromosome 7 and comparative mapping to the cattle and human genomes

A preliminary radiation hybrid (RH) map containing 50 loci on chromosome 7 of the domestic river buffalo Bubalus bubalis (BBU; 2n = 50) was constructed based on a comparative mapping approach. The RH map of BBU7 includes thirty-seven gene markers and thirteen microsatellites. All loci have been previously assigned to Bos taurus (BTA) chromosome BTA6, which is known for its association with several economically important milk production traits in cattle. The map consists of two linkage groups spanning a total length of 627.9 cR(5,000). Comparative analysis of the BBU7 RH 5,000 map with BTA6 in cattle gave new evidence for strong similarity between the two chromosomes over their entire length and exposed minor differences in locus order. Comparison of the BBU7 RH 5,000 map with the Homo sapiens (HSA) genome revealed similarity with a large chromosome segment of HSA4. Comparative analysis of loci in both species revealed more variability than previously known in gene order and several chromosome rearrangements including centromere relocation. The data obtained in our study define the evolutionarily conserved segment on BBU7 and HSA4 to be between 3.5 megabases (Mb) and 115.8 Mb in the HSA4 (genome build 36) DNA sequence. Copyright (c) 2008 S. Karger AG, Basel.

A radiation hybrid map of river buffalo (Bubalus bubalis) chromosome 1 (BBU1)

The largest chromosome in the river buffalo karyotype, BBU1, is a submetacentric chromosome with reported homology between BBU1q and bovine chromosome 1 and between BBU1p and BTA27. We present the first radiation hybrid map of this chromosome containing 69 cattle derived markers including 48 coding genes, 17 microsatellites and four ESTs distributed in two linkage groups spanning a total length of 1330.1 cR(5000). The RH map was constructed based on analysis of a recently developed river buffalo-hamster whole genome radiation hybrid (BBURH5000) panel. The retention frequency of individual markers across the panel ranged from 17.8 to 52.2%. With few exceptions, the order of markers within linkage groups is identical to the order established for corresponding cattle RH maps. The BBU1 map provides a starting point for comparison of gene order rearrangements between river buffalo chromosome 1 and its bovine homologs. Copyright (C) 2007 S. Karger AG, Basel.

First radiation hybrid map of the river buffalo X chromosome (BBUX) and comparison with BTAX

We report the first radiation hybrid map of the river buffalo X chromosome generated from a recently constructed river buffalo (Bubalus bubalis) whole-genome radiation hybrid panel (BBURH5000). This map contains a total of 33 cattle-derived markers, including 10 genes, four ESTs and 19 microsatellites. The markers are distributed in two linkage groups: LG1 contains eight markers spanning 125.6 cR, and LG2 contains 25 markers spanning 366.3 cR. LG1 contains six markers in common with bovine sequence assembly BUILD 3.1. With the exception of BMS2152, the order of these markers on our BBUX map is shuffled when compared to the cow X chromosome (Bos taurus; BTAX). From LG2, two markers (AMELX and BL22) map to a more distal portion of BTAX compared to BBUX. In addition, two pairs of LG2 markers exhibit inversions compared to BTAX (ILSTS017 and ATRX; XBM38 and PPEF1). Alternatively, when compared to the most recent bovine RH map (Bov-Gen 3000rads), BL1098 and BMS2227 from LG1 as well as PLS3 and BMS1820 from LG2 showed inverted positions on the BBUX map. These discrepancies in buffalo and cattle maps may reflect evolutionary divergence of the chromosomes or mapping errors in one of the two species. Although the set of mapped markers does not cover the entire X chromosome, this map is a starting point for the construction of a high-resolution map, which is necessary for characterization of small rearrangements that might have occurred between the Bubalus bubalis and Bos taurus X chromosomes.

A first generation whole genome RH map of the river buffalo with comparison to domestic cattle

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Preliminary radiation hybrid map for river buffalo chromosome 6 and comparison to bovine chromosome 3

We present the first radiation hybrid (RH) map of river buffalo (Bubalus bubalis) chromosome 6 (BBU6) developed with a recently constructed river buffalo whole-genome RH panel (BBURH5000). The preliminary map contains 33 cattle-derived markers, including 12 microsatellites, 19 coding genes and two ESTs, distributed across two linkage groups. Retention frequencies for markers ranged from 14.4% to 40.0%. Most of the marker orders within the linkage groups on BBU6 were consistent with the cattle genome sequence and RH maps. This preliminary RH map is the starting point for comparing gene order between river buffalo and cattle, presenting an opportunity for the examination of micro-rearrangements of these chromosomes. Also, resources for positional candidate cloning in river buffalo are enhanced.

A high-resolution radiation hybrid map of the river buffalo major histocompatibility complex and comparison with BoLA

The major histocompatibility complex (MHC) in mammals codes for antigen-presenting proteins. For this reason, the MHC is of great importance for immune function and animal health. Previous studies revealed this gene-dense and polymorphic region in river buffalo to be on the short arm of chromosome 2, which is homologous to cattle chromosome 23. Using cattle-derived STS markers and a river buffalo radiation hybrid (RH) panel (BBURH5000), we generated a high-resolution RH map of the river buffalo MHC region. The buffalo MHC RH map (cR5000) was aligned with the cattle MHC RH map (cR 12000) to compare gene order. The buffalo MHC had similar organization to the cattle MHC, with class II genes distributed in two segments, class IIa and class IIb. Class IIa was closely associated with the class I and class III regions, and class IIb was a separate cluster. A total of 53 markers were distributed into two linkage groups based on a two-point LOD score threshold of ≥8. The first linkage group included 32 markers from class IIa, class I and class III. The second linkage group included 21 markers from class IIb. Bacterial artificial chromosome clones for seven loci were mapped by fluorescence in situ hybridization on metaphase chromosomes using single- and double-color hybridizations. The order of cytogenetically mapped markers in the region corroborated the physical order of markers obtained from the RH map and served as anchor points to align and orient the linkage groups. © 2012 The Authors, Animal Genetics © 2012 Stichting International Foundation for Animal Genetics.