Genetic architecture is more complex for resistance to Septoria tritici blotch than to Fusarium head blight in Central European winter wheat

Background

Fusarium head blight (FHB) and Septoria tritici blotch (STB) severely impair wheat production. With the aim to further elucidate the genetic architecture underlying FHB and STB resistance, we phenotyped 1604 European wheat hybrids and their 135 parental lines for FHB and STB disease severities and determined genotypes at 17,372 single-nucleotide polymorphic loci.

Results

Cross-validated association mapping revealed the absence of large effect QTL for both traits. Genomic selection showed a three times higher prediction accuracy for FHB than STB disease severity for test sets largely unrelated to the training sets.

Conclusions

Our findings suggest that the genetic architecture is less complex and, hence, can be more properly tackled to perform accurate prediction for FHB than STB disease severity. Consequently, FHB disease severity is an interesting model trait to fine-tune genomic selection models exploiting beyond relatedness also knowledge of the genetic architecture.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1628-8) contains supplementary material, which is available to authorized users.     

Mapping in the era of sequencing: high density genotyping and its application for mapping TYLCV resistance in Solanum pimpinellifolium

Background

A RIL population between Solanum lycopersicum cv. Moneymaker and S. pimpinellifolium G1.1554 was genotyped with a custom made SNP array. Additionally, a subset of the lines was genotyped by sequencing (GBS).

Results

A total of 1974 polymorphic SNPs were selected to develop a linkage map of 715 unique genetic loci. We generated plots for visualizing the recombination patterns of the population relating physical and genetic positions along the genome.This linkage map was used to identify two QTLs for TYLCV resistance which contained favourable alleles derived from S. pimpinellifolium. Further GBS was used to saturate regions of interest, and the mapping resolution of the two QTLs was improved. The analysis showed highest significance on Chromosome 11 close to the region of 51.3 Mb (qTy-p11) and another on Chromosome 3 near 46.5 Mb (qTy-p3). Furthermore, we explored the population using untargeted metabolic profiling, and the most significant differences between susceptible and resistant plants were mainly associated with sucrose and flavonoid glycosides.

Conclusions

The SNP information obtained from an array allowed a first QTL screening of our RIL population. With additional SNP data of a RILs subset, obtained through GBS, we were able to perform an in silico mapping improvement to further confirm regions associated with our trait of interest. With the combination of different ~ omics platforms we provide valuable insight into the genetics of S. pimpinellifolium-derived TYLCV resistance.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1152) contains supplementary material, which is available to authorized users.     

Efficient linkage mapping using exome capture and extreme QTL in schistosome parasites

Background

Identification of parasite genes that underlie traits such as drug resistance and host specificity is challenging using classical linkage mapping approaches. Extreme QTL (X-QTL) methods, originally developed by rodent malaria and yeast researchers, promise to increase the power and simplify logistics of linkage mapping in experimental crosses of schistosomes (or other helminth parasites), because many 1000s of progeny can be analysed, phenotyping is not required, and progeny pools rather than individuals are genotyped. We explored the utility of this method for mapping a drug resistance gene in the human parasitic fluke Schistosoma mansoni.

Results

We staged a genetic cross between oxamniquine sensitive and resistant parasites, then between two F1 progeny, to generate multiple F2 progeny. One group of F2s infecting hamsters was treated with oxamniquine, while a second group was left untreated. We used exome capture to reduce the size of the genome (from 363 Mb to 15 Mb) and exomes from pooled F2 progeny (treated males, untreated males, treated females, untreated females) and the two parent parasites were sequenced to high read depth (mean = 95-366×) and allele frequencies at 14,489 variants compared. We observed dramatic enrichment of alleles from the resistant parent in a small region of chromosome 6 in drug-treated male and female pools (combined analysis: = 11.07, p = 8.74 × 10^-29). This region contains Smp_089320 a gene encoding a sulfotransferase recently implicated in oxamniquine resistance using classical linkage mapping methods.

Conclusions

These results (a) demonstrate the utility of exome capture for generating reduced representation libraries in Schistosoma mansoni, and (b) provide proof-of-principle that X-QTL methods can be successfully applied to an important human helminth. The combination of these methods will simplify linkage analysis of biomedically or biologically important traits in this parasite.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-617) contains supplementary material, which is available to authorized users.     

Construction of a high-density genetic map using specific length amplified fragment markers and identification of a quantitative trait locus for anthracnose resistance in walnut (Juglans regia L.)

Background

Walnut (Juglans regia, 2n = 32, approximately 606 Mb per 1C genome) is an economically important tree crop. Resistance to anthracnose, caused by Colletotrichum gloeosporioides, is a major objective of walnut genetic improvement in China. The recently developed specific length amplified fragment sequencing (SLAF-seq) is an efficient strategy that can obtain large numbers of markers with sufficient sequence information to construct high-density genetic maps and permits detection of quantitative trait loci (QTLs) for molecular breeding.

Results

SLAF-seq generated 161.64 M paired-end reads. 153,820 SLAF markers were obtained, of which 49,174 were polymorphic. 13,635 polymorphic markers were sorted into five segregation types and 2,577 markers of them were used to construct genetic linkage maps: 2,395 of these fell into 16 linkage groups (LGs) for the female map, 448 markers for the male map, and 2,577 markers for the integrated map. Taking into account the size of all LGs, the marker coverage was 2,664.36 cM for the female map, 1,305.58 cM for the male map, and 2,457.82 cM for the integrated map. The average intervals between two adjacent mapped markers were 1.11 cM, 2.91 cM and 0.95 cM for three maps, respectively. ‘SNP_only’ markers accounted for 89.25 % of the markers on the integrated map. Mapping markers contained 5,043 single nucleotide polymorphisms (SNPs) loci, which corresponded to two SNP loci per SLAF marker. According to the integrated map, we used interval mapping (Logarithm of odds, LOD > 3.0) to detect our quantitative trait. One QTL was detected for anthracnose resistance. The interval of this QTL ranged from 165.51 cM to 176.33 cM on LG14, and ten markers in this interval that were above the threshold value were considered to be linked markers to the anthracnose resistance trait. The phenotypic variance explained by each marker ranged from 16.2 to 19.9 %, and their LOD scores varied from 3.22 to 4.04.

Conclusions

High-density genetic maps for walnut containing 16 LGs were constructed using the SLAF-seq method with an F1 population. One QTL for walnut anthracnose resistance was identified based on the map. The results will aid molecular marker-assisted breeding and walnut resistance genes identification.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1822-8) contains supplementary material, which is available to authorized users.     

A bi-filtering method for processing single nucleotide polymorphism array data improves the quality of genetic map and accuracy of quantitative trait locus mapping in doubled haploid populations of polyploid Brassica napus

Background

Single nucleotide polymorphism (SNP) markers have a wide range of applications in crop genetics and genomics. Due to their polyploidy nature, many important crops, such as wheat, cotton and rapeseed contain a large amount of repeat and homoeologous sequences in their genomes, which imposes a huge challenge in high-throughput genotyping with sequencing and/or array technologies. Allotetraploid Brassica napus (AACC, 2n = 4x = 38) comprises of two highly homoeologous sub-genomes derived from its progenitor species B. rapa (AA, 2n = 2x = 20) and B. oleracea (CC, 2n = 2x = 18), and is an ideal species to exploit methods for reducing the interference of extensive inter-homoeologue polymorphisms (mHemi-SNPs and Pseudo-simple SNPs) between closely related sub-genomes.

Results

Based on a recent B. napus 6K SNP array, we developed a bi-filtering procedure to identify unauthentic lines in a DH population, and mHemi-SNPs and Pseudo-simple SNPs in an array data matrix. The procedure utilized both monomorphic and polymorphic SNPs in the DH population and could effectively distinguish the mHemi-SNPs and Pseudo-simple SNPs that resulted from superposition of the signals from multiple SNPs. Compared with conventional procedure for array data processing, the bi-filtering method could minimize the pseudo linkage relationship caused by the mHemi-SNPs and Pseudo-simple SNPs, thus improving the quality of SNP genetic map. Furthermore, the improved genetic map could increase the accuracies of mapping of QTLs as demonstrated by the ability to eliminate non-real QTLs in the mapping population.

Conclusions

The bi-filtering analysis of the SNP array data represents a novel approach to effectively assigning the multi-loci SNP genotypes in polyploid B. napus and may find wide applications to SNP analyses in polyploid crops.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1559-4) contains supplementary material, which is available to authorized users.     

DNAH5 is associated with total lung capacity in chronic obstructive pulmonary disease

Background

Chronic obstructive pulmonary disease (COPD) is characterized by expiratory flow limitation, causing air trapping and lung hyperinflation. Hyperinflation leads to reduced exercise tolerance and poor quality of life in COPD patients. Total lung capacity (TLC) is an indicator of hyperinflation particularly in subjects with moderate-to-severe airflow obstruction. The aim of our study was to identify genetic variants associated with TLC in COPD.

Methods

We performed genome-wide association studies (GWASs) in white subjects from three cohorts: the COPDGene Study; the Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints (ECLIPSE); and GenKOLS (Bergen, Norway). All subjects were current or ex-smokers with at least moderate airflow obstruction, defined by a ratio of forced expiratory volume in 1 second to forced vital capacity (FEV₁/FVC) <0.7 and FEV₁ < 80% predicted on post-bronchodilator spirometry. TLC was calculated by using volumetric computed tomography scans at full inspiration (TLC_CT). Genotyping in each cohort was completed, with statistical imputation of additional markers. To find genetic variants associated with TLC_CT, linear regression models were used, with adjustment for age, sex, pack-years of smoking, height, and principal components for genetic ancestry. Results were summarized using fixed-effect meta-analysis.

Results

Analysis of a total of 4,543 COPD subjects identified one genome-wide significant locus on chromosome 5p15.2 (rs114929486, β = 0.42L, P = 4.66 × 10⁻⁸).

Conclusions

In COPD, TLC_CT was associated with a SNP in dynein, axonemal, heavy chain 5 (DNAH5), a gene in which genetic variants can cause primary ciliary dyskinesia. DNAH5 could have an effect on hyperinflation in COPD.

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-014-0097-y) contains supplementary material, which is available to authorized users.     

Identification of male-specific amh duplication,sexually differentially expressed genes and microRNAs at early embryonic development of Nile tilapia (Oreochromis niloticus)

Background

The probable influence of genes and the environment on sex determination in Nile tilapia suggests that it should be regarded as a complex trait. Detection of sex determination genes in tilapia has both scientific and commercial importance. The main objective was to detect genes and microRNAs that were differentially expressed by gender in early embryonic development.

Results

Artificial fertilization of Oreochromis niloticus XX females with either sex-reversed ΔXX males or genetically-modified YY ‘supermales’ resulted in all-female and all-male embryos, respectively. RNA of pools of all-female and all-male embryos at 2, 5 and 9 dpf were used as template for a custom Agilent eArray hybridization and next generation sequencing. Fifty-nine genes differentially expressed between genders were identified by a false discovery rate of p < 0.05. The most overexpressed genes were amh and tspan8 in males, and cr/20β-hsd, gpa33, rtn4ipl and zp3 in females (p < 1 × 10⁻⁹). Validation of gene expression using qPCR in embryos and gonads indicated copy number variation in tspan8, gpa33, cr/20β-hsd and amh. Sequencing of amh identified a male-specific duplication of this gene, denoted amhy, differing from the sequence of amh by a 233 bp deletion on exonVII, hence lacking the capability to encode the protein motif that binds to the transforming growth factor beta receptor (TGF-β domain). amh and amhy segregated in the mapping family in full concordance with SD-linked marker on LG23 signifying the QTL for SD. We discovered 831 microRNAs in tilapia embryos of which nine had sexually dimorphic expression patterns by a false discovery rate of p < 0.05. An up-regulated microRNA in males, pma-mir-4585, was characterized with all six predicted target genes including cr/20β-hsd, down-regulated in males.

Conclusions

This study reports the first discovery of sexually differentially expressed genes and microRNAs at a very early stage of tilapia embryonic development, i.e. from 2 dpf. Genes with sexually differential expression patterns are enriched for copy number variation. A novel male-specific duplication of amh, denoted amhy, lacking the TGF-β domain was identified and mapped to the QTL region on LG23 for SD, thus indicating its potential role in SD.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-774) contains supplementary material, which is available to authorized users.     

Identification of a common risk haplotype for canine idiopathic epilepsy in the ADAM23 gene

Background

Idiopathic epilepsy is a common neurological disease in human and domestic dogs but relatively few risk genes have been identified to date. The seizure characteristics, including focal and generalised seizures, are similar between the two species, with gene discovery facilitated by the reduced genetic heterogeneity of purebred dogs. We have recently identified a risk locus for idiopathic epilepsy in the Belgian Shepherd breed on a 4.4 megabase region on CFA37.

Results

We have expanded a previous study replicating the association with a combined analysis of 157 cases and 179 controls in three additional breeds: Schipperke, Finnish Spitz and Beagle (p_c = 2.9e–07, p_GWAS = 1.74E-02). A targeted resequencing of the 4.4 megabase region in twelve Belgian Shepherd cases and twelve controls with opposite haplotypes identified 37 case-specific variants within the ADAM23 gene. Twenty-seven variants were validated in 285 cases and 355 controls from four breeds, resulting in a strong replication of the ADAM23 locus (p_raw = 2.76e–15) and the identification of a common 28 kb-risk haplotype in all four breeds. Risk haplotype was present in frequencies of 0.49–0.7 in the breeds, suggesting that ADAM23 is a low penetrance risk gene for canine epilepsy.

Conclusions

These results implicate ADAM23 in common canine idiopathic epilepsy, although the causative variant remains yet to be identified. ADAM23 plays a role in synaptic transmission and interacts with known epilepsy genes, LGI1 and LGI2, and should be considered as a candidate gene for human epilepsies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1651-9) contains supplementary material, which is available to authorized users.     

Impaired macrophage phagocytosis of bacteria in severe asthma

Background

Bacteria are frequently cultured from sputum samples of severe asthma patients suggesting a defect in bacterial clearance from the airway. We measured the capacity of macrophages from patients with asthma to phagocytose bacteria.

Methods

Phagocytosis of fluorescently-labelled polystyrene beads, Haemophilus influenzae or Staphylococcus aureus by broncholaveolar lavage alveolar macrophages (AM) and by monocyte-derived macrophages (MDM) from non-asthmatics, mild-moderate and severe asthmatic patients was assessed using fluorimetry.

Results

There were no differences in phagocytosis of polystyrene beads by AMs or MDMs from any of the subject groups. There was reduced phagocytosis of Haemophilus influenzae and Staphylococcus aureus in MDMs from patients with severe asthma compared to non-severe asthma (p < 0.05 and p < 0.01, respectively) and healthy subjects (p < 0.01and p < 0.001, respectively). Phagocytosis of Haemophilus influenzae and Staphylococcus aureus by AM was also reduced in severe asthma compared to normal subjects (p < 0.05). Dexamethasone and formoterol did not suppress phagocytosis of bacteria by MDMs from any of the groups.

Conclusions

Persistence of bacteria in the lower airways may result partly from a reduced phagocytic capacity of macrophages for bacteria. This may contribute to increased exacerbations, airway colonization and persistence of inflammation.     

Secretome analysis reveals effector candidates associated with broad host range necrotrophy in the fungal plant pathogen Sclerotinia sclerotiorum

Background

The white mold fungus Sclerotinia sclerotiorum is a devastating necrotrophic plant pathogen with a remarkably broad host range. The interaction of necrotrophs with their hosts is more complex than initially thought, and still poorly understood.

Results

We combined bioinformatics approaches to determine the repertoire of S. sclerotiorum effector candidates and conducted detailed sequence and expression analyses on selected candidates. We identified 486 S. sclerotiorum secreted protein genes expressed in planta, many of which have no predicted enzymatic activity and may be involved in the interaction between the fungus and its hosts. We focused on those showing (i) protein domains and motifs found in known fungal effectors, (ii) signatures of positive selection, (iii) recent gene duplication, or (iv) being S. sclerotiorum-specific. We identified 78 effector candidates based on these properties. We analyzed the expression pattern of 16 representative effector candidate genes on four host plants and revealed diverse expression patterns.

Conclusions

These results reveal diverse predicted functions and expression patterns in the repertoire of S. sclerotiorum effector candidates. They will facilitate the functional analysis of fungal pathogenicity determinants and should prove useful in the search for plant quantitative disease resistance components active against the white mold.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-336) contains supplementary material, which is available to authorized users.     

Genome-wide association mapping of quantitative resistance to sudden death syndrome in soybean

Background

Sudden death syndrome (SDS) is a serious threat to soybean production that can be managed with host plant resistance. To dissect the genetic architecture of quantitative resistance to the disease in soybean, two independent association panels of elite soybean cultivars, consisting of 392 and 300 unique accessions, respectively, were evaluated for SDS resistance in multiple environments and years. The two association panels were genotyped with 52,041 and 5,361 single nucleotide polymorphisms (SNPs), respectively. Genome-wide association mapping was carried out using a mixed linear model that accounted for population structure and cryptic relatedness.

Result

A total of 20 loci underlying SDS resistance were identified in the two independent studies, including 7 loci localized in previously mapped QTL intervals and 13 novel loci. One strong peak of association on chromosome 18, associated with all disease assessment criteria across the two panels, spanned a physical region of 1.2 Mb around a previously cloned SDS resistance gene (GmRLK18-1) in locus Rfs2. An additional variant independently associated with SDS resistance was also found in this genomic region. Other peaks were within, or close to, sequences annotated as homologous to genes previously shown to be involved in plant disease resistance. The identified loci explained an average of 54.5% of the phenotypic variance measured by different disease assessment criteria.

Conclusions

This study identified multiple novel loci and refined the map locations of known loci related to SDS resistance. These insights into the genetic basis of SDS resistance can now be used to further enhance durable resistance to SDS in soybean. Additionally, the associations identified here provide a basis for further efforts to pinpoint causal variants and to clarify how the implicated genes affect SDS resistance in soybean.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-809) contains supplementary material, which is available to authorized users.     

Association of 25-Hydroxyvitamin D status and genetic variation in the vitamin D metabolic pathway with FEV1 in the Framingham Heart Study

Background

Vitamin D is associated with lung function in cross-sectional studies, and vitamin D inadequacy is hypothesized to play a role in the pathogenesis of chronic obstructive pulmonary disease. Further data are needed to clarify the relation between vitamin D status, genetic variation in vitamin D metabolic genes, and cross-sectional and longitudinal changes in lung function in healthy adults.

Methods

We estimated the association between serum 25-hydroxyvitamin D [25(OH)D] and cross-sectional forced expiratory volume in the first second (FEV₁) in Framingham Heart Study (FHS) Offspring and Third Generation participants and the association between serum 25(OH)D and longitudinal change in FEV₁ in Third Generation participants using linear mixed-effects models. Using a gene-based approach, we investigated the association between 241 SNPs in 6 select vitamin D metabolic genes in relation to longitudinal change in FEV₁ in Offspring participants and pursued replication of these findings in a meta-analyzed set of 4 independent cohorts.

Results

We found a positive cross-sectional association between 25(OH)D and FEV₁ in FHS Offspring and Third Generation participants (P = 0.004). There was little or no association between 25(OH)D and longitudinal change in FEV₁ in Third Generation participants (P = 0.97). In Offspring participants, the CYP2R1 gene, hypothesized to influence usual serum 25(OH)D status, was associated with longitudinal change in FEV₁ (gene-based P < 0.05). The most significantly associated SNP from CYP2R1 had a consistent direction of association with FEV₁ in the meta-analyzed set of replication cohorts, but the association did not reach statistical significance thresholds (P = 0.09).

Conclusions

Serum 25(OH)D status was associated with cross-sectional FEV₁, but not longitudinal change in FEV₁. The inconsistent associations may be driven by differences in the groups studied. CYP2R1 demonstrated a gene-based association with longitudinal change in FEV₁ and is a promising candidate gene for further studies.

Electronic supplementary material

The online version of this article (doi:10.1186/s12931-015-0238-y) contains supplementary material, which is available to authorized users.