Systems genetics approach reveals candidate genes for parasite resistance from quantitative trait loci studies in agricultural species

A systems genetics approach combining pathway analysis of quantitative trait loci (QTL) and gene expression information has provided strong evidence for common pathways associated with genetic resistance to internal parasites. Gene data, collected from published QTL regions in sheep, cattle, mice, rats and humans, and microarray data from sheep, were converted to human Entrez Gene IDs and compared to the KEGG pathway database. Selection of pathways from QTL data was based on a selection index that ensured that the selected pathways were in all species and the majority of the projects overall and within species. Pathways with either up- and down-regulated genes, primarily up-regulated genes or primarily down-regulated genes, were selected from gene expression data. After comparing the data sets independently, the pathways from each data set were compared and the common set of pathways and genes was identified. Comparisons within data sets identified 21 pathways from QTL data and 66 pathways from gene expression data. Both selected sets were enriched with pathways involved in immune functions, disease and cell responses to signals. The analysis identified 14 pathways that were common between QTL and gene expression data, and four directly associated with IFNγ or MHCII, with 31 common genes, including three MHCII genes. In conclusion, a systems genetics approach combining data from multiple QTL and gene expression projects led to the discovery of common pathways associated with genetic resistance to internal parasites. This systems genetics approach may prove significant for the discovery of candidate genes for many other multifactorial, economically important traits.     

A functional genomics approach to evaluate candidate genes located in a QTL interval for milk production traits on BTA6

The potential genetic and economic advantage of marker-assisted selection for enhanced production in dairy cattle has provided an impetus to conduct numerous genome scans in order to identify associations between DNA markers and future productive potential. One area of focus has been a quantitative trait locus on bovine chromosome 6 (BTA6) found to be associated with milk yield, milk protein and fat percentage, which has been subsequently fine-mapped to six positional candidate genes. Subsequent investigations have yet to resolve which of the potential positional candidate genes is responsible for the observed associations with productive performance. In this study, we analysed candidate gene expression and the effects of gene knockdown on expression of β- and κ-casein mRNA in a small interfering RNA transfected bovine in vitro mammosphere model. From our expression studies in vivo , we observed that four of the six candidates ( ABCG2 , SPP1 , PKD2 and LAP3 ) exhibited differential expression in bovine mammary tissue over the lactation cycle, but in vitro functional studies indicate that inhibition of only one gene, SPP1 , had a significant impact on milk protein gene expression. These data suggest that the gene product of SPP1 (also known as osteopontin) has a significant role in the modulation of milk protein gene expression. While these findings do not exclude other positional candidates from influencing lactation, they support the hypothesis that the gene product of SPP1 is a significant lactational regulatory molecule.     

Genome-wide association mapping in Norwegian Red cattle identifies quantitative trait loci for fertility and milk production on BTA12

Reproductive performance is a critical trait in dairy cattle. Poor reproductive performance leads to prolonged calving intervals, higher culling rates and extra expenses related to multiple inseminations, veterinary treatments and replacements. Genetic gain for improved reproduction through traditional selection is often slow because of low heritability and negative correlations with production traits. Detection of DNA markers associated with improved reproductive performance through genome-wide association studies could lead to genetic gain that is more balanced between fertility and production. Norwegian Red cattle are well suited for such studies, as very large numbers of detailed reproduction records are available. We conducted a genome-wide association study for non-return rate, fertility treatments and retained placenta using almost 1 million records on these traits and 17 343 genome-wide single-nucleotide polymorphisms. Genotyping costs were minimized by genotyping the sires of the cows recorded and by using daughter averages as phenotypes. The genotyped sires were assigned to either a discovery or a validation population. Associations were only considered to be validated if they were significant in both groups. Strong associations were found and validated on chromosomes 1, 5, 8, 9, 11 and 12. Several of these were highly supported by findings in other studies. The most important result was an association for non-return rate in heifers in a region of BTA12 where several associations for milk production traits have previously been found. Subsequent fine-mapping verified the presence of a quantitative trait loci (QTL) having opposing effects on non-return rate and milk production at 18 Mb. The other reproduction QTL did not have pleiotropic effects on milk production, and these are therefore of considerable interest for use in marker-assisted selection.     

Construction of a cucumber genetic linkage map with SRAP markers and location of the genes for lateral branch traits

Using SRAP (sequence-related amplified polymorphism) markers a genetic linkage map of cucumber was constructed with a population consisting of 138 F₂ individuals derived from a cross of the two cucumber lines, S06 and S52. In the survey of parental polymorphisms with 182 primer combinations, 64 polymorphism-revealing primer pairs were screened out, which generated totally 108 polymorphic bands with an average of 1.7 bands per primer pair and at most 6 bands from one primer pair. The constructed molecular linkage map included 92 loci, distributed in seven linkage groups and spanning 1164.2 cM in length with an average genetic distance of 12.6 cM between two neighboring loci. Based on this linkage map, the quantitative trait loci (QTL) for the lateral branch number (lbn) and the lateral branch average length (lbl) in cucumber were identified by QTLMapper1.6. A major QTL lbn1 located between ME11SA4B and ME5EM5 in LG2 could explain 10.63% of the total variation with its positively effecting allele from S06. A major QTL lbl1 located between DC1OD3 and DC1EM14 in LG2 could account for 10.38% of the total variation with its positively effecting allele from S₀₆.     

Whole-genome scan to detect QTL for milk production, conformation, fertility and functional traits in two US Holstein families

A genome scan was conducted in two US Holstein half-sib families to identify quantitative trait loci (QTL) affecting milk production and conformation traits using the granddaughter design. The sires of the two studied families were related as sire and son and had 96 and 212 sons respectively. A total of 221 microsatellite loci were scored in both families. Statistical analysis was performed using two different analytical methods; half-sib least squares regression and Bayesian Monte Carlo Markov Chain. Traits analysed included five traditional milk production traits, somatic cell count, daughter pregnancy rate, male fertility and 20 conformation traits. A total of 47 tests achieved at least genome-wise significance. However, results from the two methods of analysis were only concordant for QTL location and level of significance in eight instances.     

奶牛产奶性状与乳房炎相关基因CpG含量及分布特征的比较分析

DNA甲基化是表观遗传学的重要组成部分,基因启动子区及第一外显子区的CpG甲基化通常抑制该基因的表达,而去甲基化则促进基因表达。已有的研究发现荷斯坦牛的乳房炎指标SCC(Somatic cell count)与产奶量呈较强负相关。文章分析并比较了这两类性状的相关基因的启动子区、第一外显子、下游2 000 bp序列中CpG含量及分布特征。结果表明,乳房炎相关基因的启动子、第一外显子中CpG含量显著低于产奶性状相关基因,而两类性状基因下游2 000 bp序列中CpG含量无显著性差异。另外,文中提出了两个量化基因序列中CpG特征的指标,一个是CpG平均距离,用来衡量序列中的CpG分布;另一个是条件概率p(G|C),用以量化序列中二核苷酸CpG随碱基C出现的可能性,并对两类基因的启动子和第一外显子区域的这两个指标做了统计检验。研究结果对产奶性状与乳房炎相关基因的DNA甲基化调控研究奠定了基础。     

Detection of loci affecting milk production and health traits in an elite US Holstein population using microsatellite markers

Quantitative trait loci (QTL) affecting health and milk production traits were studied in seven large half-sib US Holstein families by using the granddaughter design. Genotyping for 16 markers was completed and marker allele differences within and pooled-across families were analysed. Potential QTL were identified for somatic cell score (SCS), fat yield, fat percentage, protein yield and protein percentage. Three markers (BM203, BM4505 and BM2078) were associated with significant effects for different traits and, after further analysis, may be useful in marker-assisted selection in specific families. Comparisons between these data and previously identified QTL support the location of a QTL for milk yield and protein yield on chromosome 21.     

Association analysis of adiponectin and somatostatin polymorphisms on BTA1 with growth and carcass traits in Angus cattle

This study tested positional candidate genes adiponectin (ADIPOQ) and somatostatin (SST) for effects on carcass traits in a commercially relevant cattle population. Both genes are located within a region of BTA1 previously reported to harbour quantitative trait loci (QTL) that affect marbling, quality grade, yield grade, ribeye area and weaning weight in Bos taurus x Bos indicus crosses. Except for the first intron of ADIPOQ, both genes, including over 2 kb upstream of the promoters, were sequenced in five registered Angus sires to identify polymorphisms. A variable copy duplication and three single nucleotide polymorphisms (SNPs) in ADIPOQ and one SNP in SST were genotyped and tested for association with 19 traits in a 14-generation pedigree of 1697 registered Angus artificial insemination sires representing all the major USA lineages of the breed. Linear models that parameterized predicted genetic merits in terms of allele substitution effects were fit by weighted least squares, and goodness-of-fit tests were employed to differentiate causal mutations or polymorphisms in strong linkage disequilibrium (LD) with causal mutations from markers in weak LD with QTL. We confirmed the presence of QTL affecting marbling, ribeye muscle area and fat thickness in the vicinity of SST and ADIPOQ on BTA1 in Angus; excluded SST as underlying the ribeye muscle area QTL; and excluded ADIPOQ as underlying the marbling score QTL. However, association analysis provides very limited information about QTL location and has little intrinsic value when performed in the absence of linkage or LD analysis using flanking marker data to localize the QTL effect relative to positional candidate genes.     

Associations of functional candidate genes derived from gene-expression profiles of prenatal porcine muscle tissue with meat quality and muscle deposition

Ten genes (ANK1, bR10D1, CA3, EPOR, HMGA2, MYPN, NME1, PDGFRA, ERC1, TTN), whose candidacy for meat-quality and carcass traits arises from their differential expression in prenatal muscle development, were examined for association in 1700 performance-tested fattening pigs of commercial purebred and crossbred herds of Duroc, Pietrain, Pietrain x (Landrace x Large White), Duroc x (Landrace x Large White) as well as in an experimental F(2) population based on a reciprocal cross of Duroc and Pietrain. Comparative sequencing revealed polymorphic sites segregating across commercial breeds. Genetic mapping results corresponded to pre-existing assignments to porcine chromosomes or current human-porcine comparative maps. Nine of these genes showed association with meat-quality and carcass traits at a nominal P-value of < or = 0.05; PDGFRA revealed no association reaching the P < or = 0.05 threshold. In particular, HMGA2, CA3, EPOR, NME1 and TTN were associated with meat colour, pH and conductivity of loin 24 h postmortem; CA3 and MYPN exhibited association with ham weight and lean content (FOM) respectively at P-values of < 0.003 that correspond to false discovery rates of < 0.05. However, none of the genes showed significant associations for a particular trait across all populations. The study revealed statistical-genetic evidence for association of the functional candidate genes with traits related to meat quality and muscle deposition. The polymorphisms detected are not likely causal, but markers were identified that are in linkage disequilibrium with causal genetic variation within particular populations.     

Confirmation and refinement of a QTL on BTA5 affecting milk production traits in the Fleckvieh dual purpose cattle breed

We analysed a QTL affecting milk yield (MY), milk protein yield (PY) and milk fat yield (FY) in the dual purpose cattle breed Fleckvieh on BTA5. Twenty-six microsatellite markers covering 135 cM were selected to analyse nine half-sib families containing 605 sons in a granddaughter design. We thereby assigned two new markers to the public linkage map using the CRI-MAP program. Phenotypic records were daughter yield deviations (DYD) originating from the routinely performed genetic evaluations of breeding animals. To determine the position of the QTL, three different approaches were applied: interval mapping (IM), linkage analysis by variance component analysis (LAVC), and combined linkage disequilibrium (LD) and linkage (LDL) analysis. All three methods mapped the QTL in the same marker interval ( BM2830-ETH152 ) with the greatest test-statistic value at 118, 119.33 and 119.33 cM respectively. The positive QTL allele simultaneously increases DYD in the first lactation by 272 kg milk, 7.1 kg milk protein and 7.0 kg milk fat. Although the mapping accuracy and the significance of a QTL effect increased from IM over LAVC to LDL, the confidence interval was large (13, 20 and 24 cM for FY, MY and PY respectively) for the positional cloning of the causal gene. The estimated averages of pair wise marker LD with a distance <5 cM were low (0.107) and reflect the large effective population size of the Fleckvieh subpopulation analysed. This low level of LD suggests a need for increase in marker density in following fine mapping steps.     

The effect of neurospecific knockdown of candidate genes for locomotor behavior and sound production in Drosophila melanogaster

Molecular mechanisms underlying the functioning of central pattern generators (CPGs) are poorly understood. Investigations using genetic approaches in the model organism Drosophila may help to identify unknown molecular players participating in the formation or control of motor patterns. Here we report Drosophila genes as candidates for involvement in the neural mechanisms responsible for motor functions, such as locomotion and courtship song. Twenty-two Drosophila lines, used for gene identification, were isolated from a previously created collection of 1064 lines, each carrying a P element insertion in one of the autosomes. The lines displayed extreme deviations in locomotor and/or courtship song parameters compared with the whole collection. The behavioral consequences of CNS-specific RNAi-mediated knockdowns for 10 identified genes were estimated. The most prominent changes in the courtship song interpulse interval (IPI) were seen in flies with Sps2 or CG15630 knockdown. Glia-specific knockdown of these genes produced no effect on the IPI. Estrogen-induced knockdown of CG15630 in adults reduced the IPI. The product of the CNS-specific gene, CG15630 (a predicted cell surface receptor), is likely to be directly involved in the functioning of the CPG generating the pulse song pattern. Future studies should ascertain its functional role in the neurons that constitute the song CPG. Other genes (Sps2, CG34460), whose CNS-specific knockdown resulted in IPI reduction, are also worthy of detailed examination.     

SNP detection and genetic mapping of porcine genes encoding enzymes in hepatic metabolic pathways and evaluation of linkage with carcass traits

We have previously identified and mapped porcine expressed sequence tags (ESTs) derived from genes that are preferentially expressed in liver. The aim of the present study was to identify single nucleotide polymorphisms (SNPs) in porcine genes encoding enzymes in hepatic metabolic pathways and use the SNPs for mapping. Furthermore, these genes, which are involved in utilization and partitioning of nutrients, were examined for their effects on carcass and meat quality traits by linkage analyses. In total, 100 ESTs were screened for SNPs by single strand conformation polymorphism analyses across a diverse panel of animals with a 36% success rate. Twelve of 36 polymorphic loci segregated in a three-generation Duroc x Berlin Miniature Pig (F2) resource population, the DUMI resource population, and were genetically mapped. Interval mapping of the corresponding chromosomes was performed to verify mapping of the genes within quantitative trait loci (QTL) regions detected in this resource population. QTL with genome-wide significance were detected in the vicinity of GNMT, ESTL147 and HGD. These loci therefore are positional candidate genes.     

A pilot replication study of two PER3 single nucleotide polymorphisms as potential genetic markers for morning and evening earliness-lateness

Polymorphisms in genes of circadian system family seem to be of most importance for understanding of mechanisms underlying self-assessed individual variation in morning–evening preference. A review of earlier reported positive findings indicated that, at least, four polymorphisms in period circadian clock 3 (PER3) showed significant association with, at least, one of sub-constructs of a morningness–eveningness scale. However, similar to other candidate gene studies, these studies suffer from increased likelihood of false positive findings. We tried to replicate some of the most recently published positive results on associations of sub-traits of morningness–eveningness with two PER3 non-synonymous single nucleotide polymorphisms. DNA from buccal swabs was collected from healthy residents of three Russian cities, Moscow (N = 149) and Novosibirsk and Stavropol (N = 248). The tested hypotheses were formulated in accord with the earlier reported positive findings: the rare alleles might be linked to a higher score on (i) morning earliness–lateness scale (rs2640909, Moscow data-set) and (ii) both morning and evening earliness–lateness scales (rs228729, Novosibirsk and Stavropol datasets). The results provided support for the former hypothesis. These and earlier reported results highlighted several critical issues that remained to be addressed in future independent replications of positive findings on potential genetic markers for morning and evening earliness–lateness.     

Mitochondrial DNA of the sea anemone,Metridium senile (Cnidaria): Prokaryote-like genes for tRNAf-Met and small-subunit ribosomal RNA,and standard genetic code specificities for AGR and ATA codons

The nucleotide sequence of a segment of the mitochondrial DNA (mtDNA) molecule of the sea anemone Metridium senile (phylum Cnidaria, class Anthozoa, order Actiniaria) has been determined, within which have been identified the genes for respiratory chain NADH dehydrogenase subunit 2 (ND2), the small-subunit rRNA (s-rRNA), cytochrome c oxidase subunit II(COII), ND4, ND6, cytochrome b (Cyt b), tRNA^f-Met, and the large-subunit rRNA (1-rRNA). The eight genes are arranged in the order given and are all transcribed from the same strand of the molecule. The overall order of the M. senile mt-genes differs from that of other metazoan mtDNAs. In M. senile mt-protein genes, AGA and AGG codons appear to have the standard genetic code specification of arginine, rather than serine as found for other invertebrate mt-genetic codes. Also, ATA has the standard genetic code specification of isoleucine. TGA occurs in three M. senile mt-protein genes and may specify tryptophan as in other metazoan, protozoan, and some fungal mt-genetic codes. The M. senile mt-rRNA^f-Met gene has primary and secondary structure features closely resembling those of the Escherichia coli initiator tRNA, including standard dihydrouridine and TC loop sequences and a mismatch pair at the top of the aminoacyl stem. Determinations of the 5 and 3 end nucleotides of the M. senile mt-srRNAs indicated that these molecules have a homogenous size of 1,081 ntp, larger than any other known metazoan mt-s-rRNAs. Consistent with its larger size, the M. senile mt-s-rRNA can be folded into a secondary structure that more closely resembles that of the E. coli 16S rRNA than can any other metazoan mt-s-rRNA. These findings concerning M. senile mtDNA indicate that most of the unusual features regarding metazoan mt-genetic codes, rRNAs, and probably tRNAs developed after divergence of the Cnidarian line from the ancestral line common to other metazoa.Correspondence to: D.R. Wolstenholme     

Polymorphism analysis in genes of the somatotropic axis in Nellore cattle selected for growth

The somatotropic axis consists of genes that are involved in muscular development. These genes are potential regions of study to identify possible QTL for economically important traits in beef cattle. The aim of this study was to verify the existence of GH1, POU1F1, and GHR polymorphisms in Nellore cattle to verify the influence of selection in these mutations and to analyse the association between molecular markers and body weight at different ages, yearling hip height, carcass fat thickness and loin eye area. Six hundred forty-five animals from the Centro APTA Bovinos de Corte, were genotyped by PCR–RFLP techniques. The association analyses were performed with general mixed models taking into consideration the effect of one marker, and other model taking into consideration interactions between two molecular markers. Only the molecular markers rs81109601 on GH1 and rs109136815 on GHR were polymorphic; however, they were not found to be under selection. The association of the GHR rs109136815 marker and loin eye area was observed (p < 0.05), as well as the effect of interaction between the markers and the female body weight at 550 days of age (p < 0.04). The interaction effect should be considered in situations where the interactivity between two genes is known.     

Alpha-defensin DEFA1A3 gene copy number variation in Asians and its genetic association study in Chinese systemic lupus erythematosus patients

Neutrophil extracellular traps (NETs) were closely associated with activation of type I interferon (IFN) pathway in systemic lupus erythematosus (SLE). We aimed to study the genetic basis of NETs-DEFA1A3 copy number variations (CNV) in SLE and HapMap CHB + JPT populations by quantitative real-time PCR and whole genome sequences data. DEFA1A3 CNs did not differ significantly between SLE patients and controls. DEFA1A3 CNs ranged from 3 to 11 in CHB and 4 to 16 in JPT. The median of DEFA1A3 CNV of CHB (6 copies) was significantly lower than that of JPT (9 copies). Associations of genotype of tag SNP rs2738113 with DEFA1A3 CNs and mRNA expression of IFNα were observed in CHB and JPT populations. Our data provided a genetic reference of DEFA1A3 CNV for further studies and suggested that the genetic pathogenesis of NETs, as well as DEFA1A3 in SLE should be further evaluated, specially in different populations.     

Expression,SNV identification,linkage disequilibrium,and combined genotype association analysis of the muscle-specific gene CSRP3 in Chinese cattle

The cysteine and glycine-rich protein 3 (CSRP3) plays an important role in the myofiber differentiation. Here, we identified five SNVs in all exon and intron regions of the CSRP3 gene using DNA sequencing, PCR-RFLP and forced-PCR-RFLP methods in 554 cattle. Four of the five SNVs were significantly associated with growth performance and carcass traits of the cattle. In addition, we evaluated haplotype frequency and linkage disequilibrium coefficient of five sequence variants. The result of haplotype analysis demonstrated 28 haplotypes present in Qinchuan and two haplotypes in Chinese Holstein. Only haplotypes 1 and 8 were being shared by two populations, haplotype 14 had the highest haplotype frequency in Qinchuan (17.4%) and haplotype 8 had the highest haplotype frequency in Chinese Holstein (94.4%). Statistical analyses of combined genotypes indicated that some combined genotypes were significantly or highly significantly associated with growth and carcass traits in the Qinchuan cattle population. qPCR analyses also showed that bovine CSRP3 gene was exclusively expressed in longissimus dorsi muscle and heart tissues. The data support the high potential of the CSRP3 as a marker gene for the improvement of growth performance and carcass traits in selection programs.