Accurate measurement of gene copy number for human alpha-defensin DEFA1A3

Background

Multi-allelic copy number variants include examples of extensive variation between individuals in the copy number of important genes, most notably genes involved in immune function. The definition of this variation, and analysis of its impact on function, has been hampered by the technical difficulty of large-scale but accurate typing of genomic copy number. The copy-variable alpha-defensin locus DEFA1A3 on human chromosome 8 commonly varies between 4 and 10 copies per diploid genome, and presents considerable challenges for accurate high-throughput typing.

Results

In this study, we developed two paralogue ratio tests and three allelic ratio measurements that, in combination, provide an accurate and scalable method for measurement of DEFA1A3 gene number. We combined information from different measurements in a maximum-likelihood framework which suggests that most samples can be assigned to an integer copy number with high confidence, and applied it to typing 589 unrelated European DNA samples. Typing the members of three-generation pedigrees provided further reassurance that correct integer copy numbers had been assigned. Our results have allowed us to discover that the SNP rs4300027 is strongly associated with DEFA1A3 gene copy number in European samples.

Conclusions

We have developed an accurate and robust method for measurement of DEFA1A3 copy number. Interrogation of rs4300027 and associated SNPs in Genome-Wide Association Study SNP data provides no evidence that alpha-defensin copy number is a strong risk factor for phenotypes such as Crohn’s disease, type I diabetes, HIV progression and multiple sclerosis.

Electronic supplementary material

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

SNP rs10248565 in HDAC9 as a novel genomic aberration biomarker of lung adenocarcinoma in non-smoking women

Background

Numerous efforts have been made to elucidate the etiology and improve the treatment of lung cancer, but the overall five-year survival rate is still only 15%. Although cigarette smoking is the primary risk factor for lung cancer, only 7% of female lung cancer patients in Taiwan have a history of smoking. Since cancer results from progressive accumulation of genetic aberrations, genomic rearrangements may be early events in carcinogenesis.

Results

In order to identify biomarkers of early-stage adenocarcinoma, the genome-wide DNA aberrations of 60 pairs of lung adenocarcinoma and adjacent normal lung tissue in non-smoking women were examined using Affymetrix Genome-Wide Human SNP 6.0 arrays. Common copy number variation (CNV) regions were identified by ≥30% of patients with copy number beyond 2 ± 0.5 of copy numbers for each single nucleotide polymorphism (SNP) and at least 100 continuous SNP variant loci. SNPs associated with lung adenocarcinoma were identified by McNemar’s test. Loss of heterozygosity (LOH) SNPs were identified in ≥18% of patients with LOH in the locus. Aberration of SNP rs10248565 at HDAC9 in chromosome 7p21.1 was identified from concurrent analyses of CNVs, SNPs, and LOH.

Conclusion

The results elucidate the genetic etiology of lung adenocarcinoma by demonstrating that SNP rs10248565 may be a potential biomarker of cancer susceptibility.     

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.     

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.     

Multiple origins of the determinate growth habit in domesticated common bean (Phaseolus vulgaris)

Background and Aims

The actual number of domestications of a crop is one of the key questions in domestication studies. Answers to this question have generally been based on relationships between wild progenitors and domesticated descendants determined with anonymous molecular markers. In this study, this question was investigated by determining the number of instances a domestication phenotype had been selected in a crop species. One of the traits that appeared during domestication of common bean (Phaseolus vulgaris) is determinacy, in which stems end with a terminal inflorescence. It has been shown earlier that a homologue of the arabidopsis TFL1 gene – PvTFL1y – controls determinacy in a naturally occurring variation of common bean.

Methods

Sequence variation was analysed for PvTFL1y in a sample of 46 wild and domesticated accessions that included determinate and indeterminate accessions.

Key Results

Indeterminate types – wild and domesticated – showed only synonymous nucleotide substitutions. Determinate types – observed only among domesticated accessions – showed, in addition to synonymous substitutions, non-synonymous substitutions, indels, a putative intron-splicing failure, a retrotransposon insertion and a deletion of the entire locus. The retrotransposon insertion was observed in 70 % of determinate cultivars, in the Americas and elsewhere. Other determinate mutants had a more restricted distribution in the Americas only, either in the Andean or in the Mesoamerican gene pool of common bean.

Conclusions

Although each of the determinacy haplotypes probably does not represent distinct domestication events, they are consistent with the multiple (seven) domestication pattern in the genus Phaseolus. The predominance of determinacy in the Andean gene pool may reflect domestication of common bean prior to maize introduction in the Andes.     

Methanol extract from Vietnamese Caesalpinia sappan induces apoptosis in HeLa cells

Background

This study evaluated the cytotoxic activity of extracts from Caesalpinia sappan heartwood against multiple cancer cell lines using an MTT cell viability assay. The cell death though induction of apoptosis was as indicated by DNA fragmentation and caspase-3 enzyme activation.

Results

A methanol extract from C. sappan (MECS) showed cytotoxic activity against several of the cancer cell lines. The most potent activity exhibited by the MECS was against HeLa cells with an IC₅₀ value of 26.5 ± 3.2 μg/mL. Treatment of HeLa cells with various MECS concentrations resulted in growth inhibition and induction of apoptosis, as indicated by DNA fragmentation and caspase-3 enzyme activation.

Conclusion

This study is the first report of the anticancer properties of the heartwood of C. sappan native to Vietnam. Our findings demonstrate that C. sappan heartwood may have beneficial applications in the field of anticancer drug discovery.     

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.     

Essential oils from two Eucalyptus from Tunisia and their insecticidal action on Orgyia trigotephras (Lepidotera,Lymantriidae)

Background

Essential oils extracted from aromatic and medicinal plants have many biological properties and are therefore an alternative to the use of synthetic products. The chemical composition of essential oils from two medicinal plants (Eucalyptus globulus and E. lehmannii) was determined and, their insecticidal effects on the third and fourth larval stages of Orgyia trigotephras were assessed.

Results

Larvae were collected from Jebel Abderrahmane (North-East of Tunisia), conserved in groups of 50/box (21 × 10 × 10 cm) at a temperature of 25°C. Larvae were tested for larvicidal activities of essential oils. Each oil was diluted in ethanol (96%) to prepare 3 test solutions (S1 = 0.05%, S2 = 0.10% and S3 = 0.50%). Essential oils were used for contact, ingestion and Olfactory actions and compared to reference products (Bacillus thuringiensis and Decis). Olfactory action of essential oils shows that larvae mortality is higher than contact action, lower than ingestion action. MTM and FTM of S3 of E. lehmannii were respectively 1 h 32 min and 1 h 39 min are higher than those of E. globulus (MTM = 51 min and FTM = 1 h 22 min 34 sec). Contact action of E. lehmannii oil shows low insecticidal activity compared to E. globulus. MTM are respectively (1 min 52 sec and 1 min 7 sec), FTM are (2 min 38 sec, 1 min 39 sec), are the shortest recorded for S3, on the third stage of larvae. The fourth stage of larvae, MTM are (2 min 20 sec and 2 min 9 sec), FTM are (3 min 25 sec, 3 min 19 sec). Ingestion action of essential oils is longer than the contact action, since the time of death exceeds 60 minutes for all species.

Conclusion

Results shows that essential oils have a toxic action on nerves leading to a disruption of vital system of insects. High toxic properties make these plant-derived compounds suitable for incorporation in integrated pest management programs.     

Identification,characterization, and utilization of single copy genes in 29 angiosperm genomes

Background

Single copy genes are common across angiosperm genomes. With the sufficiently high quality sequenced genomes, the identification of large-scale single copy genes among multiple species is possible. Although some characteristics have been reported, our study provides novel insights into single copy genes.

Results

We identified single copy genes across 29 angiosperm genomes. A significant negative correlation was found between the number of duplicate blocks and the number of single copy genes. We found that a considerable number of single copy genes are located in organelles, showing a preference for binding and catalytic activity. The analysis of effective number of codons (Nc) illustrates that single copy genes have a stronger codon bias than non-single copy genes in eudicots. The relative high expression level of single copy genes was partially confirmed by the RNA-seq data, rather than the Codon Adaptation Index (CAI). Unlike in most other species, a strongly negatively correlation occurs between Nc and GC3 among single copy genes in grass genomes. When compared to all non-single copy genes, single copy genes indicate more conservation (as indicated by Ka and Ks values). But our alternative splicing (AS) results reveal that selective constraints are weaker in single copy genes than in low copy family genes (1–10 in-paralogs) and stronger than high copy family genes (>10 in-paralogs). Using concatenated shared single copy genes, we obtained a well-resolved phylogenetic tree. With the addition of intron sequences, the branch support is improved, but striking incongruences are also evident. Therefore, it is noteworthy that inclusion of intron sequences seems more appropriate for the phylogenetic reconstruction at lower taxonomic levels.

Conclusions

Our analysis provides insight into the evolutionary characteristics of single copy genes across 29 angiosperm genomes. The results suggest that there are key differences in evolutionary constraints between single copy genes and non-single copy genes. And to some extent, these evolutionary constraints show some species-specific differences, especially between eudicots and monocots. Our preliminary evidence also suggests that the concatenated shared single copy genes are well suited for use in resolving phylogenetic relationships.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-504) 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.     

A comprehensive genetic map of sugarcane that provides enhanced map coverage and integrates high-throughput Diversity Array Technology (DArT) markers

Background

Sugarcane genetic mapping has lagged behind other crops due to its complex autopolyploid genome structure. Modern sugarcane cultivars have from 110-120 chromosomes and are in general interspecific hybrids between two species with different basic chromosome numbers: Saccharum officinarum (2n = 80) with a basic chromosome number of 10 and S. spontaneum (2n = 40-128) with a basic chromosome number of 8. The first maps that were constructed utilised the single dose (SD) markers generated using RFLP, more recent maps generated using AFLP and SSRs provided at most 60% genome coverage. Diversity Array Technology (DArT) markers are high throughput allowing greater numbers of markers to be generated.

Results

Progeny from a cross between a sugarcane variety Q165 and a S. officinarum accession IJ76-514 were used to generate 2467 SD markers. A genetic map of Q165 was generated containing 2267 markers, These markers formed 160 linkage groups (LGs) of which 147 could be placed using allelic information into the eight basic homology groups (HGs) of sugarcane. The HGs contained from 13 to 23 LGs and from 204 to 475 markers with a total map length of 9774.4 cM and an average density of one marker every 4.3 cM. Each homology group contained on average 280 markers of which 43% were DArT markers 31% AFLP, 16% SSRs and 6% SNP markers. The multi-allelic SSR and SNP markers were used to place the LGs into HGs.

Conclusions

The DArT array has allowed us to generate and map a larger number of markers than ever before and consequently to map a larger portion of the sugarcane genome. This larger number of markers has enabled 92% of the LGs to be placed into the 8 HGs that represent the basic chromosome number of the ancestral species, S. spontaneum. There were two HGs (HG2 and 8) that contained larger numbers of LGs verifying the alignment of two sets of S. officinarum chromosomes with one set of S. spontaneum chromosomes and explaining the difference in basic chromosome number between the two ancestral species. There was also evidence of more complex structural differences between the two ancestral species.

Electronic supplementary material

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

Genetic architecture of resistance to Septoria tritici blotch in European wheat

Background

Septoria tritici blotch is an important leaf disease of European winter wheat. In our survey, we analyzed Septoria tritici blotch resistance in field trials with a large population of 1,055 elite hybrids and their 87 parental lines. Entries were fingerprinted with the 9 k SNP array. The accuracy of prediction of Septoria tritici blotch resistance achieved with different genome-wide mapping approaches was evaluated based on robust cross validation scenarios.

Results

Septoria tritici blotch disease severities were normally distributed, with genotypic variation being significantly (P < 0.01) larger than zero. The cross validation study revealed an absence of large effect QTL for additive and dominance effects. Application of genomic selection approaches particularly designed to tackle complex agronomic traits allowed to double the accuracy of prediction of Septoria tritici blotch resistance compared to calculation methods suited to detect QTL with large effects.

Conclusions

Our study revealed that Septoria tritici blotch resistance in European winter wheat is controlled by multiple loci with small effect size. This suggests that the currently achieved level of resistance in this collection is likely to be durable, as involvement of a high number of genes in a resistance trait reduces the risk of the resistance to be overcome by specific pathogen isolates or races.

Electronic supplementary material

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