A genome-wide library of CB4856/N2 introgression lines of Caenorhabditis elegans

Recombinant inbred lines (RILs) derived from Caenorhabditis elegans wild-type N2 and CB4856 are increasingly being used for mapping genes underlying complex traits. To speed up mapping and gene discovery, introgression lines (ILs) offer a powerful tool for more efficient QTL identification. We constructed a library of 90 ILs, each carrying a single homozygous CB4856 genomic segment introgressed into the genetic background of N2. The ILs were genotyped by 123 single-nucleotide polymorphism (SNP) markers. The proportion of the CB4856 segments in most lines does not exceed 3%, and together the introgressions cover 96% of the CB4856 genome. The value of the IL library was demonstrated by identifying novel loci underlying natural variation in two ageing-related traits, i.e. lifespan and pharyngeal pumping rate. Bin mapping of lifespan resulted in six QTLs, which all have a lifespan-shortening effect on the CB4856 allele. We found five QTLs for the decrease in pumping rate, of which four colocated with QTLs found for average lifespan. This suggests pleiotropic or closely linked QTL associated with lifespan and pumping rate. Overall, the presented IL library provides a versatile resource toward easier and efficient fine mapping and functional analyses of loci and genes underlying complex traits in C. elegans.     

Mos1-Mediated Transgenesis to Probe Consequences of Single Gene Mutations in Variation-Rich Isolates of Caenorhabditis elegans

Caenorhabditis elegans, especially the N2 isolate, is an invaluable biological model system. Numerous additional natural C. elegans isolates have been shown to have unexpected genotypic and phenotypic variations which has encouraged researchers to use next generation sequencing methodology to develop a more complete picture of genotypic variations among the isolates. To understand the phenotypic effects of a genomic variation (GV) on a single gene, in a variation-rich genetic background, one should analyze that particular GV in a well understood genetic background. In C. elegans, the analysis is usually done in N2, which requires extensive crossing to bring in the GV. This can be a very time consuming procedure thus it is important to establish a fast and efficient approach to test the effect of GVs from different isolates in N2. Here we use a Mos1-mediated single-copy insertion (MosSCI) method for phenotypic assessments of GVs from the variation-rich Hawaiian strain CB4856 in N2. Specifically, we investigate effects of variations identified in the CB4856 strain on tac-1 which is an essential gene that is necessary for mitotic spindle elongation and pronuclear migration. We show the usefulness of the MosSCI method by using EU1004 tac-1(or402) as a control. or402 is a temperature sensitive lethal allele within a well-conserved TACC domain (transforming acidic coiled-coil) that results in a leucine to phenylalanine change at amino acid 229. CB4856 contains a variation that affects the second exon of tac-1 causing a cysteine to tryptophan change at amino acid 94 also within the TACC domain. Using the MosSCI method, we analyze tac-1 from CB4856 in the N2 background and demonstrate that the C94W change, albeit significant, does not cause any obvious decrease in viability. This MosSCI method has proven to be a rapid and efficient way to analyze GVs.     

Quantitative PCR approach to SNP detection and linkage mapping in Caenorhabditis elegans

I report a method for single nucleotide polymorphism (SNP) detection and linkage mapping in Caenorhabditis elegans using automated oligonucleotide design and fluorescence-based quantitative PCR detection. Nine hundred twenty-three oligonucleotide pairs were designed to produce small products of <150 bp for efficient amplification in a PCR, with one oligonucleotide of each pair overlapping a SNP site at the 3'-most nucleotide. A subset of the pairs were tested, and efficient allelic discrimination was obtained for SNPs between N2, the canonical laboratory strain, and CB4856, a strain isolated from Hawaii commonly used for mapping studies. Linkage mapping is demonstrated using the unc-119 locus of C. elegans. This quantitative PCR method provides an inexpensive, uniform, and automatable detection alternative for genetic mapping strategies in C. elegans or other organisms.     

Mitochondrial DNA Variant in COX1 Subunit Significantly Alters Energy Metabolism of Geographically Divergent Wild Isolates in Caenorhabditis elegans

Mitochondrial DNA (mtDNA) sequence variation can influence the penetrance of complex diseases and climatic adaptation. While studies in geographically defined human populations suggest that mtDNA mutations become fixed when they have conferred metabolic capabilities optimally suited for a specific environment, it has been challenging to definitively assign adaptive functions to specific mtDNA sequence variants in mammals. We investigated whether mtDNA genome variation functionally influences Caenorhabditis elegans wild isolates of distinct mtDNA lineages and geographic origins. We found that, relative to N2 (England) wild-type nematodes, CB4856 wild isolates from a warmer native climate (Hawaii) had a unique p.A12S amino acid substitution in the mtDNA-encoded COX1 core catalytic subunit of mitochondrial complex IV (CIV). Relative to N2, CB4856 worms grown at 20 °C had significantly increased CIV enzyme activity, mitochondrial matrix oxidant burden, and sensitivity to oxidative stress but had significantly reduced lifespan and mitochondrial membrane potential. Interestingly, mitochondrial membrane potential was significantly increased in CB4856 grown at its native temperature of 25 °C. A transmitochondrial cybrid worm strain, chpIR (M, CB4856 > N2), was bred as homoplasmic for the CB4856 mtDNA genome in the N2 nuclear background. The cybrid strain also displayed significantly increased CIV activity, demonstrating that this difference results from the mtDNA-encoded p.A12S variant. However, chpIR (M, CB4856 > N2) worms had significantly reduced median and maximal lifespan relative to CB4856, which may relate to their nuclear–mtDNA genome mismatch. Overall, these data suggest that C. elegans wild isolates of varying geographic origins may adapt to environmental challenges through mtDNA variation to modulate critical aspects of mitochondrial energy metabolism.     

A Caenorhabditis elegans Wild Type Defies the Temperature–Size Rule Owing to a Single Nucleotide Polymorphism in tra-3

Ectotherms rely for their body heat on surrounding temperatures. A key question in biology is why most ectotherms mature at a larger size at lower temperatures, a phenomenon known as the temperature–size rule. Since temperature affects virtually all processes in a living organism, current theories to explain this phenomenon are diverse and complex and assert often from opposing assumptions. Although widely studied, the molecular genetic control of the temperature–size rule is unknown. We found that the Caenorhabditis elegans wild-type N2 complied with the temperature–size rule, whereas wild-type CB4856 defied it. Using a candidate gene approach based on an N2 × CB4856 recombinant inbred panel in combination with mutant analysis, complementation, and transgenic studies, we show that a single nucleotide polymorphism in tra-3 leads to mutation F96L in the encoded calpain-like protease. This mutation attenuates the ability of CB4856 to grow larger at low temperature. Homology modelling predicts that F96L reduces TRA-3 activity by destabilizing the DII-A domain. The data show that size adaptation of ectotherms to temperature changes may be less complex than previously thought because a subtle wild-type polymorphism modulates the temperature responsiveness of body size. These findings provide a novel step toward the molecular understanding of the temperature–size rule, which has puzzled biologists for decades.     

Development and evaluation of a 9K SNP array for peach by internationally coordinated SNP detection and validation in breeding germplasm

Although a large number of single nucleotide polymorphism (SNP) markers covering the entire genome are needed to enable molecular breeding efforts such as genome wide association studies, fine mapping, genomic selection and marker-assisted selection in peach [Prunus persica (L.) Batsch] and related Prunus species, only a limited number of genetic markers, including simple sequence repeats (SSRs), have been available to date. To address this need, an international consortium (The International Peach SNP Consortium; IPSC) has pursued a coordinated effort to perform genome-scale SNP discovery in peach using next generation sequencing platforms to develop and characterize a high-throughput Illumina Infinium® SNP genotyping array platform. We performed whole genome re-sequencing of 56 peach breeding accessions using the Illumina and Roche/454 sequencing technologies. Polymorphism detection algorithms identified a total of 1,022,354 SNPs. Validation with the Illumina GoldenGate® assay was performed on a subset of the predicted SNPs, verifying ∼75% of genic (exonic and intronic) SNPs, whereas only about a third of intergenic SNPs were verified. Conservative filtering was applied to arrive at a set of 8,144 SNPs that were included on the IPSC peach SNP array v1, distributed over all eight peach chromosomes with an average spacing of 26.7 kb between SNPs. Use of this platform to screen a total of 709 accessions of peach in two separate evaluation panels identified a total of 6,869 (84.3%) polymorphic SNPs.The almost 7,000 SNPs verified as polymorphic through extensive empirical evaluation represent an excellent source of markers for future studies in genetic relatedness, genetic mapping, and dissecting the genetic architecture of complex agricultural traits. The IPSC peach SNP array v1 is commercially available and we expect that it will be used worldwide for genetic studies in peach and related stone fruit and nut species.     

Genetic mapping of variation in dauer larvae development in growing populations of Caenorhabditis elegans

In the nematode Caenorhabditis elegans, the appropriate induction of dauer larvae development within growing populations is likely to be a primary determinant of genotypic fitness. The underlying genetic architecture of natural genetic variation in dauer formation has, however, not been thoroughly investigated. Here, we report extensive natural genetic variation in dauer larvae development within growing populations across multiple wild isolates. Moreover, bin mapping of introgression lines (ILs) derived from the genetically divergent isolates N2 and CB4856 reveals 10 quantitative trait loci (QTLs) affecting dauer formation. Comparison of individual ILs to N2 identifies an additional eight QTLs, and sequential IL analysis reveals six more QTLs. Our results also show that a behavioural, laboratory-derived, mutation controlled by the neuropeptide Y receptor homolog npr-1 can affect dauer larvae development in growing populations. These findings illustrate the complex genetic architecture of variation in dauer larvae formation in C. elegans and may help to understand how the control of variation in dauer larvae development has evolved.     

Genetic Dissection of Drought and Heat Tolerance in Chickpea through Genome-Wide and Candidate Gene-Based Association Mapping Approaches

To understand the genetic basis of tolerance to drought and heat stresses in chickpea, a comprehensive association mapping approach has been undertaken. Phenotypic data were generated on the reference set (300 accessions, including 211 mini-core collection accessions) for drought tolerance related root traits, heat tolerance, yield and yield component traits from 1–7 seasons and 1–3 locations in India (Patancheru, Kanpur, Bangalore) and three locations in Africa (Nairobi, Egerton in Kenya and Debre Zeit in Ethiopia). Diversity Array Technology (DArT) markers equally distributed across chickpea genome were used to determine population structure and three sub-populations were identified using admixture model in STRUCTURE. The pairwise linkage disequilibrium (LD) estimated using the squared-allele frequency correlations (r²; when r²<0.20) was found to decay rapidly with the genetic distance of 5 cM. For establishing marker-trait associations (MTAs), both genome-wide and candidate gene-sequencing based association mapping approaches were conducted using 1,872 markers (1,072 DArTs, 651 single nucleotide polymorphisms [SNPs], 113 gene-based SNPs and 36 simple sequence repeats [SSRs]) and phenotyping data mentioned above employing mixed linear model (MLM) analysis with optimum compression with P3D method and kinship matrix. As a result, 312 significant MTAs were identified and a maximum number of MTAs (70) was identified for 100-seed weight. A total of 18 SNPs from 5 genes (ERECTA, 11 SNPs; ASR, 4 SNPs; DREB, 1 SNP; CAP2 promoter, 1 SNP and AMDH, 1SNP) were significantly associated with different traits. This study provides significant MTAs for drought and heat tolerance in chickpea that can be used, after validation, in molecular breeding for developing superior varieties with enhanced drought and heat tolerance.     

A pipeline for high throughput detection and mapping of SNPs from EST databases

Single nucleotide polymorphisms (SNPs) represent the most abundant type of genetic variation that can be used as molecular markers. The SNPs that are hidden in sequence databases can be unlocked using bioinformatic tools. For efficient application of these SNPs, the sequence set should be error-free as much as possible, targeting single loci and suitable for the SNP scoring platform of choice. We have developed a pipeline to effectively mine SNPs from public EST databases with or without quality information using QualitySNP software, select reliable SNP and prepare the loci for analysis on the Illumina GoldenGate genotyping platform. The applicability of the pipeline was demonstrated using publicly available potato EST data, genotyping individuals from two diploid mapping populations and subsequently mapping the SNP markers (putative genes) in both populations. Over 7000 reliable SNPs were identified that met the criteria for genotyping on the GoldenGate platform. Of the 384 SNPs on the SNP array approximately 12% dropped out. For the two potato mapping populations 165 and 185 SNPs segregating SNP loci could be mapped on the respective genetic maps, illustrating the effectiveness of our pipeline for SNP selection and validation.     

A Multi-Trait,Meta-analysis for Detecting Pleiotropic Polymorphisms for Stature,Fatness and Reproduction in Beef Cattle

Polymorphisms that affect complex traits or quantitative trait loci (QTL) often affect multiple traits. We describe two novel methods (1) for finding single nucleotide polymorphisms (SNPs) significantly associated with one or more traits using a multi-trait, meta-analysis, and (2) for distinguishing between a single pleiotropic QTL and multiple linked QTL. The meta-analysis uses the effect of each SNP on each of n traits, estimated in single trait genome wide association studies (GWAS). These effects are expressed as a vector of signed t-values (t) and the error covariance matrix of these t values is approximated by the correlation matrix of t-values among the traits calculated across the SNP (V). Consequently, t''V⁻¹t is approximately distributed as a chi-squared with n degrees of freedom. An attractive feature of the meta-analysis is that it uses estimated effects of SNPs from single trait GWAS, so it can be applied to published data where individual records are not available. We demonstrate that the multi-trait method can be used to increase the power (numbers of SNPs validated in an independent population) of GWAS in a beef cattle data set including 10,191 animals genotyped for 729,068 SNPs with 32 traits recorded, including growth and reproduction traits. We can distinguish between a single pleiotropic QTL and multiple linked QTL because multiple SNPs tagging the same QTL show the same pattern of effects across traits. We confirm this finding by demonstrating that when one SNP is included in the statistical model the other SNPs have a non-significant effect. In the beef cattle data set, cluster analysis yielded four groups of QTL with similar patterns of effects across traits within a group. A linear index was used to validate SNPs having effects on multiple traits and to identify additional SNPs belonging to these four groups.     

Recombinational Landscape and Population Genomics of Caenorhabditis elegans

Recombination rate and linkage disequilibrium, the latter a function of population genomic processes, are the critical parameters for mapping by linkage and association, and their patterns in Caenorhabditis elegans are poorly understood. We performed high-density SNP genotyping on a large panel of recombinant inbred advanced intercross lines (RIAILs) of C. elegans to characterize the landscape of recombination and, on a panel of wild strains, to characterize population genomic patterns. We confirmed that C. elegans autosomes exhibit discrete domains of nearly constant recombination rate, and we show, for the first time, that the pattern holds for the X chromosome as well. The terminal domains of each chromosome, spanning about 7% of the genome, exhibit effectively no recombination. The RIAILs exhibit a 5.3-fold expansion of the genetic map. With median marker spacing of 61 kb, they are a powerful resource for mapping quantitative trait loci in C. elegans. Among 125 wild isolates, we identified only 41 distinct haplotypes. The patterns of genotypic similarity suggest that some presumed wild strains are laboratory contaminants. The Hawaiian strain, CB4856, exhibits genetic isolation from the remainder of the global population, whose members exhibit ample evidence of intercrossing and recombining. The population effective recombination rate, estimated from the pattern of linkage disequilibrium, is correlated with the estimated meiotic recombination rate, but its magnitude implies that the effective rate of outcrossing is extremely low, corroborating reports of selection against recombinant genotypes. Despite the low population, effective recombination rate and extensive linkage disequilibrium among chromosomes, which are techniques that account for background levels of genomic similarity, permit association mapping in wild C. elegans strains.     

Genome-Wide Association Study Dissects the Genetic Architecture of Seed Weight and Seed Quality in Rapeseed (Brassica napus L.)

Association mapping can quickly and efficiently dissect complex agronomic traits. Rapeseed is one of the most economically important polyploid oil crops, although its genome sequence is not yet published. In this study, a recently developed 60K Brassica Infinium^® SNP array was used to analyse an association panel with 472 accessions. The single-nucleotide polymorphisms (SNPs) of the array were in silico mapped using ‘pseudomolecules’ representative of the genome of rapeseed to establish their hypothetical order and to perform association mapping of seed weight and seed quality. As a result, two significant associations on A8 and C3 of Brassica napus were detected for erucic acid content, and the peak SNPs were found to be only 233 and 128 kb away from the key genes BnaA.FAE1 and BnaC.FAE1. BnaA.FAE1 was also identified to be significantly associated with the oil content. Orthologues of Arabidopsis thaliana HAG1 were identified close to four clusters of SNPs associated with glucosinolate content on A9, C2, C7 and C9. For seed weight, we detected two association signals on A7 and A9, which were consistent with previous studies of quantitative trait loci mapping. The results indicate that our association mapping approach is suitable for fine mapping of the complex traits in rapeseed.     

Fluorescence detection of single nucleotide polymorphisms using nucleic acid probe containing tricyclic base-linked acyclonucleoside

A reliable and simple method for detecting nucleobase mutations is very important clinically because sequence variations in human DNA cause genetic diseases and genetically influenced traits. A majority of sequence variations are attributed to single nucleotide polymorphisms (SNPs). Here, we developed a method for SNP detection using DNA probes that contained a fluorescent tricyclic base-linked acyclonucleoside N. The type of nucleobases involved in the SNP sites in an RNA target could be determined using four DNA probes containing N. Further, we found that the SNP in the RNA target could be detected by a visible color. Thus, this system would provide a novel and simple method for detecting SNPs in an RNA target.     

A Caenorhabditis elegans wild type defies the temperature-size rule owing to a single nucleotide polymorphism in tra-3

Ectotherms rely for their body heat on surrounding temperatures. A key question in biology is why most ectotherms mature at a larger size at lower temperatures, a phenomenon known as the temperature-size rule. Since temperature affects virtually all processes in a living organism, current theories to explain this phenomenon are diverse and complex and assert often from opposing assumptions. Although widely studied, the molecular genetic control of the temperature-size rule is unknown. We found that the Caenorhabditis elegans wild-type N2 complied with the temperature-size rule, whereas wild-type CB4856 defied it. Using a candidate gene approach based on an N2 x CB4856 recombinant inbred panel in combination with mutant analysis, complementation, and transgenic studies, we show that a single nucleotide polymorphism in tra-3 leads to mutation F96L in the encoded calpain-like protease. This mutation attenuates the ability of CB4856 to grow larger at low temperature. Homology modelling predicts that F96L reduces TRA-3 activity by destabilizing the DII-A domain. The data show that size adaptation of ectotherms to temperature changes may be less complex than previously thought because a subtle wild-type polymorphism modulates the temperature responsiveness of body size. These findings provide a novel step toward the molecular understanding of the temperature-size rule, which has puzzled biologists for decades.     

Development and Validation of a 20K Single Nucleotide Polymorphism (SNP) Whole Genome Genotyping Array for Apple (Malus × domestica Borkh)

High-density SNP arrays for genome-wide assessment of allelic variation have made high resolution genetic characterization of crop germplasm feasible. A medium density array for apple, the IRSC 8K SNP array, has been successfully developed and used for screens of bi-parental populations. However, the number of robust and well-distributed markers contained on this array was not sufficient to perform genome-wide association analyses in wider germplasm sets, or Pedigree-Based Analysis at high precision, because of rapid decay of linkage disequilibrium. We describe the development of an Illumina Infinium array targeting 20K SNPs. The SNPs were predicted from re-sequencing data derived from the genomes of 13 Malus × domestica apple cultivars and one accession belonging to a crab apple species (M. micromalus). A pipeline for SNP selection was devised that avoided the pitfalls associated with the inclusion of paralogous sequence variants, supported the construction of robust multi-allelic SNP haploblocks and selected up to 11 entries within narrow genomic regions of ±5 kb, termed focal points (FPs). Broad genome coverage was attained by placing FPs at 1 cM intervals on a consensus genetic map, complementing them with FPs to enrich the ends of each of the chromosomes, and by bridging physical intervals greater than 400 Kbps. The selection also included ∼3.7K validated SNPs from the IRSC 8K array. The array has already been used in other studies where ∼15.8K SNP markers were mapped with an average of ∼6.8K SNPs per full-sib family. The newly developed array with its high density of polymorphic validated SNPs is expected to be of great utility for Pedigree-Based Analysis and Genomic Selection. It will also be a valuable tool to help dissect the genetic mechanisms controlling important fruit quality traits, and to aid the identification of marker-trait associations suitable for the application of Marker Assisted Selection in apple breeding programs.     

High density SNP mapping and QTL analysis for fruit quality characteristics in peach (Prunus persica L.)

Single nucleotide polymorphisms (SNPs) were used to construct an integrated SNP linkage map of peach (Prunus persica (L.) Batsch). A set of 1,536 SNPs were evaluated with the GoldenGate® Genotyping assay in two mapping populations, Pop-DF, and Pop-DG. After genotyping and filtering, a final set of 1,400 high quality SNPs in Pop-DF and 962 in Pop-DG with full map coverage were selected and used to construct two linkage maps with JoinMap®4.0. The Pop-DF map covered 422 cM of the peach genome and included 1,037 SNP markers, and Pop-DG map covered 369 cM and included 738 SNPs. A consensus map was constructed with 588 SNP markers placed in eight linkage groups (n?=?8 for peach), with map coverage of 454 cM and an average distance of 0.81 cM/marker site. Placements of SNPs on the “peach v1.0” physical map were compared to placement on the linkage maps and several differences were observed. Using the SNP linkage map of Pop-DG and phenotypic data collected for three harvest seasons, a QTL analysis for fruit quality traits and chilling injury symptoms was carried out with the mapped SNPs. Significant QTL effects were detected for mealiness (M) and flesh bleeding (FBL) QTLs on linkage group 4 and flesh browning (FBr) on linkage group 5. This study represents one of the first examples of QTL detection for quality traits and chilling injury symptoms using a high-density SNP map in a single peach F1 family.