A search for quantitative trait loci for ovulation rate in cattle

Seventy-seven polymorphic microsatellites were analysed in offspring of three elite sires that were part of the foundation of an experimental population selected for twinning rate at the US Meat Animal Research Center, Clay Center, Nebraska. All females were assessed for ovulation rate by rectal palpation of corpora lutea over 8–10 consecutive oestrous cycles from approximately 12 to 18 months of age, and associations between ovulation rate and sire allele were examined in each of the three sire groups. A preliminary analysis was performed using selectively genotyped daughters of each sire. Markers found significant or approaching significance were also genotyped in all daughters, sons and granddaughters of these sires. A test of marker associations limited to the granddaughter data provided an independent confirmation of marker effect and significance relative to the initial test with daughter data. Putative ovulation rate quantitative trait loci were detected on chromosomes 7 and 23. Marker UWCA20 on chromosome 7 was associated with an effect in excess of one phenotypic standard deviation and accounted for approximately 10% of phenotypic variation ovulation rate. Marker CYP21 (steroid 21-hydroxylase) on chromosome 23 was associated with an effect of slightly less than half a phenotypic standard deviation and accounted for approximately 4% of phenotypic variation.     

Segregation of a major gene influencing ovulation in progeny of Lacaune meat sheep

Inheritance of the ovulation rate (OR) in the Lacaune meat breed was studied through records from a small nucleus of 36 hyper-prolific ewes screened on farms on the basis of their natural litter size, and from progeny data of three selected Lacaune sires. These sires were chosen at the AI centre according to their breeding values estimated for the mean and the variability of their daughters'' litter size. Non-carrier Lacaune dairy ewes were inseminated to produce 121 F1 daughters and 27 F1 sons. Twelve sons (four from each sire) were used in turn to inseminate non-carrier Lacaune dairy ewes providing 260 BC progeny ewes. F1 and BC progeny were brought from private farms and gathered after weaning on an experimental farm where ovulation rates were recorded in the first and second breeding seasons. With an average of 6.5 records each, the mean OR of hyper-prolific ewes was very high (5.34), and 38.4% of records showed a rate of 6 or more. F1 data showed high repeatability of OR (r = 0.54) within ewe, with significant variability among ewes. High OR (≥ 4) were observed in each family. A segregation analysis provided a significant likelihood ratio and classified the three founders as heterozygous. BC ewes also displayed high repeatability of OR (r = 0.47) and the mean OR varied considerably between families (from 1.24 to 1.78). Seven of the 12 BC families presented high-ovulating ewes (at least one record ≥ 4) and segregation analysis yielded a highly significant likelihood ratio as compared to an empirical test distribution. The high variability of the mean ovulation rate shown by a small group of daughters of BC ewes inseminated by putative carrier F1 rams, and the very high ovulation rate observed for some of these ewe lambs, confirmed the segregation of a major gene with two co-dominant alleles borne by an autosome. The difference between homozygous non-carriers and heterozygous ewes was about one ovulation on the observed scale and 2.2 standard deviations on the underlying scale.     

Genome Assembly Anchored QTL Map of Bovine Chromosome 14

Bovine chromosome 14 (BTA14) has been widely explored for quantitative trait loci (QTL) and genes related to economically important traits in both dairy and beef cattle. We reviewed more than 40 investigations and anchored 126 QTL to the current genome assembly (Btau 4_0). Using this anchored QTL map, we observed that, in dairy cattle, the region spanning 0 – 10 Mb on BTA14 has the highest density QTL map with a total of 56 QTL, mainly for milk production traits. It is very likely that both somatic cell score (SCS) and clinical mastitis share some common QTL in two regions: 61.48 Mb - 73.84 Mb and 7.86 Mb – 39.55 Mb, respectively. As well, both ovulation rate and twinning rate might share a common QTL region from 34.16 Mb to 65.38 Mb. However, there are no common QTL locations in three pregnancy related phenotypes: non-return rate, pregnancy rate and daughter pregnancy rate. In beef cattle, the majority of QTL are located in a broad region of 15 Mb – 45 Mb on the chromosome. Functional genes, such as CRH, CYP11B1, DGAT1, FABP4 and TG, as potential candidates for some of these QTL, were also reviewed. Therefore, our review provides a standardized QTL map anchored within the current genome assembly, which would enhance the process of selecting positional and physiological candidate genes for many important traits in cattle.     

New insight into the SSC8 genetic determination of fatty acid composition in pigs

Background

Fat content and fatty acid composition in swine are becoming increasingly studied because of their effect on sensory and nutritional quality of meat. A QTL (quantitative trait locus) for fatty acid composition in backfat was previously detected on porcine chromosome 8 (SSC8) in an Iberian x Landrace F₂ intercross. More recently, a genome-wide association study detected the same genomic region for muscle fatty acid composition in an Iberian x Landrace backcross population. ELOVL6, a strong positional candidate gene for this QTL, contains a polymorphism in its promoter region (ELOVL6:c.-533C < T), which is associated with percentage of palmitic and palmitoleic acids in muscle and adipose tissues. Here, a combination of single-marker association and the haplotype-based approach was used to analyze backfat fatty acid composition in 470 animals of an Iberian x Landrace F₂ intercross genotyped with 144 SNPs (single nucleotide polymorphisms) distributed along SSC8.

Results

Two trait-associated SNP regions were identified at 93 Mb and 119 Mb on SSC8. The strongest statistical signals of both regions were observed for palmitoleic acid (C16:1(n-7)) content and C18:0/C16:0 and C18:1(n-7)/C16:1(n-7) elongation ratios. MAML3 and SETD7 are positional candidate genes in the 93 Mb region and two novel microsatellites in MAML3 and nine SNPs in SETD7 were identified. No significant association for the MAML3 microsatellite genotypes was detected. The SETD7:c.700G > T SNP, although statistically significant, was not the strongest signal in this region. In addition, the expression of MAML3 and SETD7 in liver and adipose tissue varied among animals, but no association was detected with the polymorphisms in these genes. In the 119 Mb region, the ELOVL6:c.-533C > T polymorphism showed a strong association with percentage of palmitic and palmitoleic fatty acids and elongation ratios in backfat.

Conclusions

Our results suggest that the polymorphisms studied in MAML3 and SETD7 are not the causal mutations for the QTL in the 93 Mb region. However, the results for ELOVL6 support the hypothesis that the ELOVL6:c.-533C > T polymorphism has a pleiotropic effect on backfat and intramuscular fatty acid composition and that it has a role in the determination of the QTL in the 119 Mb region.     

SMAD3 rs17228212 Gene Polymorphism Is Associated with Reduced Risk to Cerebrovascular Accidents and Subclinical Atherosclerosis in Anti-CCP Negative Spanish Rheumatoid Arthritis Patients

Rheumatoid arthritis (RA) is a complex polygenic inflammatory disease associated with accelerated atherosclerosis and increased risk of cardiovascular (CV) disease. Previous genome-wide association studies have described SMAD3 rs17228212 polymorphism as an important signal associated with CV events. The aim of the present study was to evaluate for the first time the relationship between this gene polymorphism and the susceptibility to CV manifestations and its potential association with the presence of subclinical atherosclerosis assessed by the evaluation of carotid intima-media thickness (cIMT) in patients with RA.

Methods

One thousand eight hundred and ninety-seven patients fulfilling classification criteria for RA were genotyped for SMAD3 rs17228212 gene polymorphism through TaqMan genotyping assay. Also, subclinical atherosclerosis determined by the assessment of cIMT was analyzed in a subgroup of these patients by carotid ultrasonography.

Results

No statistically significant differences were observed when allele frequencies of RA patients with or without CV events were compared. Nevertheless, when RA patients were stratified according to anti-cyclic citrullinated peptide (anti-CCP) status, we found that in RA patients who were negative for anti-CCP antibodies, the presence of C allele of SMAD3 rs17228212 polymorphism conferred a protective effect against the risk of cerebrovascular accident (CVA) after adjustment for demographic and classic CV risk factors (HR [95%CI]=0.36 [0.14–0.94], p=0.038) in a Cox regression model. Additionally, correlation between the presence of C allele of SMAD3 rs17228212 polymorphism and lower values of cIMT was found after adjustment for demographic and classic CV risk factors (p-value=0.0094) in the anti-CCP negative RA patients.

Conclusions

Our results revealed that SMAD3 rs17228212 gene variant is associated with lower risk of CVA and less severe subclinical atherosclerosis in RA patients negative for anti-CCP antibodies. These findings may have importance to establish predictive models of CV disease in RA patients according to anti-CCP status.     

An MspI polymorphism at the bovine growth hormone (bGH) gene is linked to a locus affecting milk protein percentage

SSCP analysis of the bovine growth hormone (bGH) gene in Israel Holstein dairy cattle uncovered five intragenic haplotypes, denoted A to E. Of these, Haplotype E differed from the others at six fragments; one of which corresponded to the polymorphic MspI site in intron III, at which haplotype E carried the disabled MspI (-) allele. Haplotype E was observed in a single sire only, carrying haplotype A as the second bGH allele. In 523 daughters of this sire genotyped for the MspI polymorphism, heterozygous (+/-) as compared to homozygous (+/+) daughters, showed a significant increasing effect on protein percentage and kg protein per year; and a decreasing effect (P < 0.10) on milk somatic cell counts (MSSC). None of the daughters were homozygous (-/-), indicating that the frequency of this allele in the general population was essentially zero. Calculated skewness (g1) values for the two daughter groups differed significantly with (+/-) daughters showing negative skewness (in the direction of lower protein percentage), and (+/+) daughters positive skewness (in the direction of higher protein percentage). The direction of skewness in each group is indicative of the presence of a QTL having an increasing effect on milk protein percentage in coupling linkage with the MspI (-) allele in this sire, but at some distance from it. Maximum likelihood estimates of the proportion of recombination (r) between the putative QTL and bGH, and the allele substitution effect at the QTL (d), were r = 0.33, a = 0.07% protein, with standard errors 0.058 and 0.009% protein, respectively.     

The BCL2 -938C>A Promoter Polymorphism Is Associated with Risk for and Time to Aseptic Loosening of Total Hip Arthroplasty

Aseptic loosening is a major cause of revision surgery of total hip arthroplasty (THA). Only few host factors affecting aseptic loosening have been identified until now, although they are urgently needed to identify and possibly treat those patients at higher risk for aseptic loosening. To determine whether the functional single nucleotide polymorphism (SNP) c.-938C>A (rs2279115), located in the promoter region of the BCL2 gene has an impact on aseptic loosening of THA we genotyped and analyzed 234 patients suffering from aseptic loosening and 231 patients after primary THA. The polymorphism is associated with risk for aseptic loosening with the CC genotype at highest risk for aseptic loosening, Odds Ratio CC vs. AA 1.93, 95%CI 1.15–3.25, p = 0.013. In contrast, low risk AA genotype carriers that still developed aseptic loosening showed a significantly shorter time to aseptic loosening than patients carrying the C allele (p = 0.004). These results indicate that the BCL2 -938C>A polymorphism influences the occurrence and course of aseptic loosening and suggests this polymorphism as an interesting candidate for prospective studies and analyses in THA registers.     

A Single Nucleotide Polymorphism within the Interferon Gamma Receptor 2 Gene Perfectly Coincides with Polledness in Holstein Cattle

Polledness is a high impact trait in modern milk and beef production to meet the demands of animal welfare and work safety. Previous studies have mapped the polled-locus to the proximal region of the bovine chromosome 1 (BTA1) and narrowed it down to approximately 1 Mb. Sequencing of the positional candidate genes within the 1 Mb polled region and whole genome sequencing of Holsteins revealed a single nucleotide polymorphism (SNP) AC000158: g.1390292G>A within intron 3 of the interferon gamma receptor 2 gene (IFNGR2) in perfect co-segregation with polledness in Holsteins. This complete association was validated in 443 animals of the same breed. This SNP allows reliable genotyping of horned, heterozygous and homozygous polled Holsteins, even in animals that could not be resolved using the previously published haplotype for Holstein.     

Characterization of SMAD3 Gene Variants for Possible Roles in Ventricular Septal Defects and Other Congenital Heart Diseases

Background

Nodal/TGF signaling pathway has an important effect at early stages of differentiation of human embryonic stem cells in directing them to develop into different embryonic lineages. SMAD3 is a key intracellular messenger regulating factor in the Nodal/TGF signaling pathway, playing important roles in embryonic and, particularly, cardiovascular system development. The aim of this work was to find evidence on whether SMAD3 variations might be associated with ventricular septal defects (VSD) or other congenital heart diseases (CHD).

Methods

We sequenced the SMAD3 gene for 372 Chinese Han CHD patients including 176 VSD patients and evaluated SNP rs2289263, which is located before the 5’UTR sequence of the gene. The statistical analyses were conducted using Chi-Square Tests as implemented in SPSS (version 13.0). The Hardy-Weinberg equilibrium test of the population was carried out using the online software OEGE.

Results

Three heterozygous variants in SMAD3 gene, rs2289263, rs35874463 and rs17228212, were identified. Statistical analyses showed that the rs2289263 variant located before the 5’UTR sequence of SMAD3 gene was associated with the risk of VSD (P value=0.013 <0.05).

Conclusions

The SNP rs2289263 in the SMAD3 gene is associated with VSD in Chinese Han populations.     

Survey of genome sequences in a wild sweet potato,Ipomoea trifida (H. B. K.) G. Don

Ipomoea trifida (H. B. K.) G. Don. is the most likely diploid ancestor of the hexaploid sweet potato, I. batatas (L.) Lam. To assist in analysis of the sweet potato genome, de novo whole-genome sequencing was performed with two lines of I. trifida, namely the selfed line Mx23Hm and the highly heterozygous line 0431-1, using the Illumina HiSeq platform. We classified the sequences thus obtained as either ‘core candidates’ (common to the two lines) or ‘line specific’. The total lengths of the assembled sequences of Mx23Hm (ITR_r1.0) was 513 Mb, while that of 0431-1 (ITRk_r1.0) was 712 Mb. Of the assembled sequences, 240 Mb (Mx23Hm) and 353 Mb (0431-1) were classified into core candidate sequences. A total of 62,407 (62.4 Mb) and 109,449 (87.2 Mb) putative genes were identified, respectively, in the genomes of Mx23Hm and 0431-1, of which 11,823 were derived from core sequences of Mx23Hm, while 28,831 were from the core candidate sequence of 0431-1. There were a total of 1,464,173 single-nucleotide polymorphisms and 16,682 copy number variations (CNVs) in the two assembled genomic sequences (under the condition of log₂ ratio of >1 and CNV size >1,000 bases). The results presented here are expected to contribute to the progress of genomic and genetic studies of I. trifida, as well as studies of the sweet potato and the genus Ipomoea in general.     

Precisely mapping a major gene conferring resistance to Hessian fly in bread wheat using genotyping-by-sequencing

Background

One of the reasons hard red winter wheat cultivar ‘Duster’ (PI 644016) is widely grown in the southern Great Plains is that it confers a consistently high level of resistance to biotype GP of Hessian fly (Hf). However, little is known about the genetic mechanism underlying Hf resistance in Duster. This study aimed to unravel complex structures of the Hf region on chromosome 1AS in wheat by using genotyping-by-sequencing (GBS) markers and single nucleotide polymorphism (SNP) markers.

Results

Doubled haploid (DH) lines generated from a cross between two winter wheat cultivars, ‘Duster’ and ‘Billings’ , were used to identify genes in Duster responsible for effective and consistent resistance to Hf. Segregation in reaction of the 282 DH lines to Hf biotype GP fit a one-gene model. The DH population was genotyped using 2,358 markers developed using the GBS approach. A major QTL, explaining 88% of the total phenotypic variation, was mapped to a chromosome region that spanned 178 cM and contained 205 GBS markers plus 1 SSR marker and 1 gene marker, with 0.86 cM per marker in genetic distance. The analyses of GBS marker sequences and further mapping of SSR and gene markers enabled location of the QTL-containing linkage group on the short arm of chromosome 1A. Comparative mapping of the common markers for the gene for QHf.osu-1A^d in Duster and the Hf-resistance gene for QHf.osu-1A⁷⁴ in cultivar ‘2174’ showed that the two Hf resistance genes are located on the same chromosome arm 1AS, only 11.2 cM apart in genetic distance. The gene at QHf.osu-1A^d in Duster has been delimited within a 2.7 cM region.

Conclusion

Two distinct resistance genes exist on the short arm of chromosome 1A as found in the two hard red winter cultivars, 2174 and Duster. Whereas the Hf resistance gene in 2174 is likely allelic to one or more of the previously mapped resistance genes (H9, H10, H11, H16, or H17) in wheat, the gene in Duster is novel and confers a more consistent phenotype than 2174 in response to biotype GP infestation in controlled-environment assays.

Electronic supplementary material

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

Chromosomal regions underlying noncoagulation of milk in Finnish Ayrshire cows

About 10% of Finnish Ayrshire cows produce noncoagulating milk, i.e., milk that does not form a curd in a standard 30-min testing time and is thus a poor raw material for cheese dairies. This phenomenon is associated with peak and midlactation, but some cows produce noncoagulating milk persistently. A genomewide scan under a selective DNA pooling method was carried out to locate genomic regions associated with the noncoagulation of milk. On the basis of the hypothesis of the same historical mutation, we pooled the data across sires. Before testing pools for homogeneity, allele intensities were corrected for PCR artifacts, i.e., shadow bands and differential amplification. Results indicating association were verified using daughter design and selective genotyping within families. Data consisted of 18 sire families with 477 genotyped daughters in total, i.e., 12% of each tail of the milk coagulation ability. Data were analyzed using interval mapping under maximum-likelihood and nonparametric methods. BMS1126 on chromosome 2 and BMS1355 on chromosome 18 were associated with noncoagulation of milk across families on an experimentwise 0.1% significance level. By scanning gene databases, we found two potential candidate genes: LOC538897, a nonspecific serine/threonine kinase on chromosome 2, and SIAT4B, a sialyltransferase catalyzing the last step of glycosylation of κ-casein on chromosome 18. Further studies to determine the role of the candidates in the noncoagulation of milk are clearly needed.     

Preferential Localization of the Limb-Girdle Muscular Dystrophy Type 2A Gene in the Proximal Part of a 1-cM 15q15.1-q15.3 Interval

A gene for a recessive form of limb-girdle muscular dystrophy (LGMD2A) has been localized to chromosome 15. A physical map of the 7-cM candidate 15q15.1-q21.1 region has been constructed by means of a 10–12-Mb continuum of overlapping YAC clones. New microsatellite markers developed from these YACs were genotyped on large, consanguineous LGMD2A pedigrees from different origins. The identification of recombination events in these families allowed the restriction of the LGMD2A region to an estimated 1-cM interval, equivalent to ~3–4 Mb. Linkage disequilibrium data on genetic isolates from the island of Réunion and from the Amish community suggest a preferential location of the LGMD2A gene in the proximal part of this region. Analysis of the interrelated pedigrees from Réunion revealed the existence of at least six different carrier haplotypes. This allelic heterogeneity is incompatible with the presumed existence of a founder effect and suggests that multiple LGMD2A mutations may segregate in this population.     

Identification of an ovulation rate QTL in cattle on BTA14 using selective DNA pooling and interval mapping

Increased twinning incidence in beef cattle has the potential to improve production efficiency. However, phenotypic selection for twinning rate is difficult because of the trait's low heritability and the long time interval necessary to collect phenotypic records. Therefore, this trait and the correlated trait of ovulation rate are ideal candidates for marker-assisted selection. The objective of this study was to complete a genome-wide search for ovulation rate quantitative trait loci (QTL) in two related sire families. The families (paternal halfsib sires 839802 and 839803) were from a population of cattle selected for ovulation rate at the USDA Meat Animal Research Center, Clay Center, Nebraska. Putative ovulation rate QTL have previously been identified in the 839802 family on chromosomes 7 and 19; however, marker coverage in the original scan was not complete. This study fills the gaps in marker coverage of the earlier study by adding approximately 60 informative microsatellites to each sire family. Each family was genotyped using selective DNA pooling. Sons and daughters were included in either the high or low pool based on their estimated breeding value deviations from the mid-parent average (EBVMD) for ovulation rate. Approximately 40% (839802) and 26% (839803) of available progeny comprised the high and low pools combined. Pooled typing revealed possible associations (nominal P < 0.05) between ovulation rate and marker genotype for 11 and 15 microsatellites in the 839802 and 839803 families, respectively. Subsequent interval mapping strengthened support for the presence of an ovulation rate QTL on BTA14 (chromosome-wise P < 0.02).     

Genome-wide association and linkage analyses localize a progressive retinal atrophy locus in Persian cats

Hereditary eye diseases of animals serve as excellent models of human ocular disorders and assist in the development of gene and drug therapies for inherited forms of blindness. Several primary hereditary eye conditions affecting various ocular tissues and having different rates of progression have been documented in domestic cats. Gene therapy for canine retinopathies has been successful, thus the cat could be a gene therapy candidate for other forms of retinal degenerations. The current study investigates a hereditary, autosomal recessive, retinal degeneration specific to Persian cats. A multi-generational pedigree segregating for this progressive retinal atrophy was genotyped using a 63 K SNP array and analyzed via genome-wide linkage and association methods. A multi-point parametric linkage analysis localized the blindness phenotype to a ~1.75 Mb region with significant LOD scores (Z ≈ 14, θ = 0.00) on cat chromosome E1. Genome-wide TDT, sib-TDT, and case–control analyses also consistently supported significant association within the same region on chromosome E1, which is homologous to human chromosome 17. Using haplotype analysis, a ~1.3 Mb region was identified as highly associated for progressive retinal atrophy in Persian cats. Several candidate genes within the region are reasonable candidates as a potential causative gene and should be considered for molecular analyses.     

AZFc Deletions and Spermatogenic Failure: A Population-Based Survey of 20,000 Y Chromosomes

Deletions involving the Y chromosome’s AZFc region are the most common known genetic cause of severe spermatogenic failure (SSF). Six recurrent interstitial deletions affecting the region have been reported, but their population genetics are largely unexplored. We assessed the deletions’ prevalence in 20,884 men in five populations and found four of the six deletions (presented here in descending order of prevalence): gr/gr, b2/b3, b1/b3, and b2/b4. One of every 27 men carried one of these four deletions. The 1.6 Mb gr/gr deletion, found in one of every 41 men, almost doubles the risk of SSF and accounts for ∼2% of SSF, although <2% of men with the deletion are affected. The 1.8 Mb b2/b3 deletion, found in one of every 90 men, does not appear to be a risk factor for SSF. The 1.6 Mb b1/b3 deletion, found in one of every 994 men, appears to increase the risk of SSF by a factor of 2.5, although <2% of men with the deletion are affected, and it accounts for only 0.15% of SSF. The 3.5 Mb b2/b4 deletion, found in one of every 2,320 men, increases the risk of SSF 145 times and accounts for ∼6% of SSF; the observed prevalence should approximate the rate at which the deletion arises anew in each generation. We conclude that a single rare variant of major effect (the b2/b4 deletion) and a single common variant of modest effect (the gr/gr deletion) are largely responsible for the AZFc region’s contribution to SSF in the population.     

Searching for Genomic Region of High-Fat Diet-Induced Type 2 Diabetes in Mouse Chromosome 2 by Analysis of Congenic Strains

SMXA-5 mice are a high-fat diet-induced type 2 diabetes animal model established from non-diabetic SM/J and A/J mice. By using F2 intercross mice between SMXA-5 and SM/J mice under feeding with a high-fat diet, we previously mapped a major diabetogenic QTL (T2dm2sa) on chromosome 2. We then produced the congenic strain (SM.A-T2dm2sa (R0), 20.8–163.0 Mb) and demonstrated that the A/J allele of T2dm2sa impaired glucose tolerance and increased body weight and body mass index in the congenic strain compared to SM/J mice. We also showed that the combination of T2dm2sa and other diabetogenic loci was needed to develop the high-fat diet-induced type 2 diabetes. In this study, to narrow the potential genomic region containing the gene(s) responsible for T2dm2sa, we constructed R1 and R2 congenic strains. Both R1 (69.6–163.0 Mb) and R2 (20.8–128.2 Mb) congenic mice exhibited increases in body weight and abdominal fat weight and impaired glucose tolerance compared to SM/J mice. The R1 and R2 congenic analyses strongly suggested that the responsible genes existed in the overlapping genomic interval (69.6–128.2 Mb) between R1 and R2. In addition, studies using the newly established R1A congenic strain showed that the narrowed genomic region (69.6–75.4 Mb) affected not only obesity but also glucose tolerance. To search for candidate genes within the R1A genomic region, we performed exome sequencing analysis between SM/J and A/J mice and extracted 4 genes (Itga6, Zak, Gpr155, and Mtx2) with non-synonymous coding SNPs. These four genes might be candidate genes for type 2 diabetes caused by gene-gene interactions. This study indicated that one of the genes responsible for high-fat diet-induced diabetes exists in the 5.8 Mb genomic interval on mouse chromosome 2.     

Comparative mapping of a region on chromosome 10 containing QTL for reproduction in swine

Several quantitative trait loci (QTL) for important reproductive traits (age of puberty, ovulation rate, nipple number and plasma FSH) have been identified on the long arm of porcine chromosome 10. Bi-directional chromosome painting has shown that this region is homologous to human chromosome 10p. Because few microsatellite or type I markers have been placed on SSC10, we wanted to increase the density of known ESTs mapped in this region of the porcine genome. Genes were chosen for their position on human chromosome 10, sequence availability from the TIGR pig gene indices, and their potential as a candidate gene. The PCR primers were designed to amplify across introns or 3'-UTR to maximize single nucleotide polymorphism (SNP) discovery. Parents of the mapping population (one sire and seven dams) were amplified and sequenced to find informative markers. The SNPs were genotyped using primer extension and mass spectrometry. These amplification products were also used to probe a BAC library (RPCI-44, Roswell Park Cancer Institute) for positive clones and screened for microsatellites. Six genes from human chromosome 10p (AKR1C2, PRKCQ, ITIH2, ATP5C1, PIP5K2A and GAD2) were mapped in the MARC swine mapping population. Gene order was conserved within these markers from centromere to telomere of porcine chromosome 10q, as compared with human chromosome 10p. Four of these genes (PIP5K2A, ITIH2, GAD2 and AKR1C2), which map under QTL, are potential candidate genes. Identification of porcine homologues near important QTL and development of a comparative map for this chromosome will allow further fine- mapping and positional cloning of candidate genes affecting reproductive traits.     

Deploying QTL-seq for rapid delineation of a potential candidate gene underlying major trait-associated QTL in chickpea

A rapid high-resolution genome-wide strategy for molecular mapping of major QTL(s)/gene(s) regulating important agronomic traits is vital for in-depth dissection of complex quantitative traits and genetic enhancement in chickpea. The present study for the first time employed a NGS-based whole-genome QTL-seq strategy to identify one major genomic region harbouring a robust 100-seed weight QTL using an intra-specific 221 chickpea mapping population (desi cv. ICC 7184 × desi cv. ICC 15061). The QTL-seq-derived major SW QTL (CaqSW1.1) was further validated by single-nucleotide polymorphism (SNP) and simple sequence repeat (SSR) marker-based traditional QTL mapping (47.6% R² at higher LOD >19). This reflects the reliability and efficacy of QTL-seq as a strategy for rapid genome-wide scanning and fine mapping of major trait regulatory QTLs in chickpea. The use of QTL-seq and classical QTL mapping in combination narrowed down the 1.37 Mb (comprising 177 genes) major SW QTL (CaqSW1.1) region into a 35 kb genomic interval on desi chickpea chromosome 1 containing six genes. One coding SNP (G/A)-carrying constitutive photomorphogenic9 (COP9) signalosome complex subunit 8 (CSN8) gene of these exhibited seed-specific expression, including pronounced differential up-/down-regulation in low and high seed weight mapping parents and homozygous individuals during seed development. The coding SNP mined in this potential seed weight-governing candidate CSN8 gene was found to be present exclusively in all cultivated species/genotypes, but not in any wild species/genotypes of primary, secondary and tertiary gene pools. This indicates the effect of strong artificial and/or natural selection pressure on target SW locus during chickpea domestication. The proposed QTL-seq-driven integrated genome-wide strategy has potential to delineate major candidate gene(s) harbouring a robust trait regulatory QTL rapidly with optimal use of resources. This will further assist us to extrapolate the molecular mechanism underlying complex quantitative traits at a genome-wide scale leading to fast-paced marker-assisted genetic improvement in diverse crop plants, including chickpea.