Case report: The effective response to pembrolizumab in combination with bevacizumab in the treatment of a recurrent glioblastoma with multiple extracranial metastases

Multiple extracranial metastases of recurrent glioblastoma are rare and often indicate a very poor prognosis. The main conventional treatments are chemotherapy, radiotherapy, chemoradiotherapy or antiangiogenic therapy. Median overall survival is 2.3 to 6 months after the detection of extracranial metastases, and to date, there is no effective treatment for these patients. Herein, we report a recurrent glioblastoma patient with lung metastasis treated with a combination therapy containing bevacizumab and pembrolizumab due to overexpression of PD-L1 and the absence of driver mutations. The progression-free survival was 11 months from lung metastases to bone metastases. This combination treatment was further used as maintenance therapy for another 11 months after bone metastasis and secondary dorsal metastasis because there was no suitable treatment alternative. The overall survival was 27 months after lung metastases, which is much longer than previously reported cases. To our knowledge, this was the first effective use of bevacizumab plus pembrolizumab in a glioblastoma patient with extracranial metastases. Furthermore, this was the first time that bevacizumab plus pembrolizumab was used as a maintenance treatment in glioblastoma, with 11 months of response. Importantly, we showed that such combination therapy may be a novel and effective therapy for glioblastoma patients with extracranial metastases.

Rapid subcutaneous progression after immunotherapy in pretreated patients with metastatic carcinoma: two case reports

There is heterogeneity in cancer patients' responses to immune checkpoint inhibitors (ICIs), including hyperprogression, which is very rapid tumor progression following immunotherapy, and pseudoprogression, which is an initial increase followed by a decrease in tumor burden or in the number of tumor lesions. This heterogeneity complicates clinical decisions because either premature withdrawal of the treatment or prolonged ineffective treatment harms patients. We presented two patients treated with ICIs with heterogeneous responses. One patient had Merkel cell carcinoma in the right thigh, and the other had nasopharyngeal squamous carcinoma. The first patient was treated with sintilimab and the second with sintilimab combined with abraxane. In the first patient, subcutaneous lesions grew substantially after the first cycle of treatment with sintilimab. In the second patient, subcutaneous lesions grew gradually after the second cycle of treatment with sintilimab combined with abraxane. In both cases, biopsy examination confirmed that newly emerged lesions were metastases of the primary tumor. These two cases remind clinicians that when subcutaneous nodules appear after treatment with ICIs, pathological biopsy is needed to determine the nature-pseudoprogression or rapid progression-of the disease course.

Successful Treatment of Patients with Refractory High-Grade Serous Ovarian Cancer with GOPC-ROS1 Fusion Using Crizotinib: A Case Report

BACKGROUND

Recently, multiple poly (ADP-ribose) polymerase (PARP) inhibitors have demonstrated excellent efficacy among patients with ovarian cancer with or without BRCA mutations. However, alternative therapeutic options are urgently required for patients who cannot benefit from conventional chemotherapy or PARP inhibitors.

CASE PRESENTATION

A patient with high-grade serous ovarian carcinoma presented to our clinic after developing resistance to chemotherapy. Paired tumor-normal next-generation sequencing (NGS) was performed using peripheral blood to identify potential actionable mutations. NGS revealed the patient harboring a GOPC-ROS1 fusion, which was subsequently verified using a reverse transcription polymerase chain reaction assay. No germline or somatic mutation in BRCA1/2 or mismatch repair genes was detected. Therefore, the patient received crizotinib treatment. A rapid, favorable clinical response (partial response at 1 month) was observed, with further pathological response monitored and evaluated in follow-up interrogation.

CONCLUSION

This study suggested that crizotinib was an off-the-shelf, practical, and ostensibly effective treatment option for patients with ovarian cancer with ROS1 rearrangement. NGS-based genetic testing may guide to plan therapeutic paradigms, and render precision medicine promising in ovarian cancer treatment.

IMPLICATIONS FOR PRACTICE

Despite the previous report of ROS1 fusion in patients with ovarian cancer, it remains unknown whether patients can benefit from targeted therapeutic drugs. This study reports a GOPC-ROS1 fusion identified by next-generation sequencing in a patient with chemotherapy-resistant ovarian cancer. The patient was administered crizotinib and showed rapid, remarkable response. This study suggests that comprehensive sequencing should be offered for patients with ovarian cancer without effective therapeutic strategies, and crizotinib can be used to treat ROS1-rearranged ovarian carcinomas.

Immunomodulatory features of radiotherapy in lung carcinoma

Background

As a standard treatment modality, radiotherapy (RT) has been widely employed for cancer treatment. In addition to directly killing tumor cells, RT is known for its immunomodulatory effects. Nevertheless, the effects of ionizing radiation on immune reactions, as well as their underpinning mechanisms, are complex and remain unclear.

Methods

The immunomodulatory effects of ionizing irradiation were dynamically monitored by concomitant assessment of peripheral blood lymphocytes, and Th1 (CD3+CD4+IFN-r+), Th2 (CD3+CD4+IL-4+), Tc1 (CD3+CD8+IFN-r+) and Tc2 (CD3+CD8+IL-4+) cells, along with CD25, CD28, CTLA-4, PD-1, Foxp3, TGF-β, and IL-10 gene expression levels in 30 lung cancer patients who underwent RT.

Results

Local cancer RT activated cellular immune reactions, which was reflected by an obvious reduction of B cells and increased CD8 and NK cell amounts, and consequent increase of suppressor T cells (Ts). Further investigation showed that Th2 and Tc2 responses were significantly increased while Th1 and Tc1 were decreased.

Conclusions

The immunomodulatory effects mediated by RT are characterized by a shift from humoral to cellular immunity, significant augmentation of CD8 and Ts subpopulation, and Th2 and Tc2 responses, indicating an immuno-activating response, which might be beneficial for initial antitumor immune reactions, but may not affect the later ones. Immunomodulatory therapy should be performed upon RT to restore the immune balance for maintaining a Th1/Tc1 dominant immunity to achieve long-term anticancer immunity.

[Retrospective analysis for 104 cases of early-stage Hodgkin's Lymphoma treated with different modality therapies]

This paper explored the curative effect of combined modality therapy and extended field radiotherapy for early-stage Hodgkin's Lymphoma. 104 cases of early-stage Hodgkin's Lymphoma from Jan 1987 to Dec 2010 in PLA Hospital 307 were retrospectively analyzed, including 76 cases in combined modality therapy group and 28 cases in extended field radiotherapy group, and the long-term efficacy and toxicity of two therapy modalities were evaluated. The results showed that the median survival time of 104 cases was 85.42 months, the complete remission rates of combined modality therapy and extended field radiotherapy groups were 72.4 and 71.4 respectively (P = 0.924); the overall response rates of combined modality therapy and extended field radiotherapy groups were 97.4 and 96.4 respectively (P = 0.779); the 5-year overall survival (OS) rates in the 2 groups were 89.5 and 89.1 respectively, and the 8-year OS rates of the 2 groups were 81.3 and 70.6. No statistical difference was found in above-mentioned 2 groups. Moreover, the 5-year progression free survival (PFS) rates of these 2 groups were 84.2 and 69.0 (P = 0.04), and 8-year PFS rates of these 2 groups were 80.0 and 55.5 (P = 0.04) respectively, the 5-year relapse rates of these 2 groups were 28.1 and 45.6 (P = 0.023) respectively. It is concluded that the combined modality therapy can raise the PFS rate and reduce the relapse rate as compared with extended field radiotherapy for early-stage Hodgkin's Lymphoma, but there is no difference in the overall survival rate between the 2 groups.

Identification of quantitative trait loci for shank length and growth at different development stages in chicken

Shank length affects chicken leg health and longer shanks are a source of leg problems in heavy-bodied chickens. Identification of quantitative trait loci (QTL) affecting shank length traits may be of value to genetic improvement of these traits in chickens. A genome scan was conducted on 238 F(2) chickens from a reciprocal cross between the Silky Fowl and the White Plymouth Rock breeds using 125 microsatellite markers to detect static and developmental QTL affecting weekly shank length and growth (from 1 to 12 weeks) in chickens. Static QTL affected shank length from birth to time t, while developmental QTL affected shank growth from time t-1 to time t. Seven static QTL on six chromosomes (GGA2, GGA3, GGA4, GGA7, GGA9 and GGA23) were detected at ages of 2, 3, 4, 5, 6, 7, 9 and 12 weeks, and six developmental QTL on five chromosomes (GGA1, GGA2, GGA4, GGA5 and GGA23) were detected for five shank growth periods, weeks 2-3, 4-5, 5-6, 10-11 and 11-12. A static QTL and a developmental QTL (SQSL1 and DQSL2) were identified at GGA2 (between ADL0190 and ADL0152). SQSL1 explained 2.87-5.30% of the phenotypic variation in shank length from 3 to 7 weeks. DQSL2 explained 2.70% of the phenotypic variance of shank growth between 2 and 3 weeks. Two static and two developmental QTL were involved chromosome 4 and chromosome 23. Two chromosomes (GGA7 and GGA9) had static QTL but no developmental QTL and another two chromosomes (GGA1 and GGA5) had developmental QTL but no static QTL. The results of this study show that shank length and shank growth at different developmental stages involve different QTL.

Role of selection and inbreeding on the incidence of cutaneous malignant melanoma in Sinclair swine

This paper reports the quantitative analysis of the historical database of a herd of Sinclair swine affected by cutaneous malignant melanoma. The herd was under partial and non-systematic selection for melanoma susceptibility (animals having at least one tumour during the first 6 weeks of life). Weighted selection differentials for the number of tumours at birth and the number of tumours at 6 weeks were generally positive and between -0.43 and 4.76 tumours for the number of tumours at 6 weeks. Estimates of the heritability for number of tumours at birth and at 6 weeks using 1934 animals were 0.27 (+/-0.03) and 0.25 (+/-0.03), respectively. The estimate of the genetic correlation between these two traits was 0.95 (+/-0.03). Genetic trends were positive for the number of tumours at birth and at 6 weeks. In spite of positive selection differentials and a moderate heritability, there was a negative phenotypic trend in the number of tumours. Natural selection might be acting in a direction opposite to artificial selection in the Sinclair herd. The slopes of the regression of the number of tumours at birth, at 6 weeks, and melanoma susceptibility on individual inbreeding coefficients were non-significant, indicating no evidence of dominance. The number of live-born pigs was lower in litters from parents susceptible to the disease (p < 0.01).

A genome scan for quantitative trait loci associated with body weight at different developmental stages in chickens

A genome scan to detect quantitative trait loci (QTL) affecting body weight in chickens was conducted on 238 F(2) chickens from a reciprocal cross of Silky Fowl and White Plymouth Rock using 125 microsatellite markers covering 23 autosomes and the Z chromosome. Two types of QTL were considered: static QTL (SQ) and developmental QTL (DQ). Static QTL affected the body weight from hatch to time t, and DQ affected the body weight from time t-1 to time t. Six SQ and nine DQ were detected. Of these QTL, four reached a genome-wide significance of 5% or better, with SQ1 and DQ1 being the most significant QTL. Static QTL1 was on chromosome 1 between GCT0006 and MCW0106 and explained 4.05-9.80% of the phenotypic variation in body weights from 3 to 12 weeks of age. At 9, 10 and 11 weeks, the genome-wide significance thresholds of SQ1 were <1%. Developmental QTL1 was located on chromosome 1 between MCW0168 and GCT0006, and explained 2.75% of the phenotypic variation for body weight from week 7 to 8 with a genome-wide significance level <1%. The results suggest that body weight from hatch to time t and developmental growth from time t-1 to time t may involve two different sets of genes or gene actions.

Allelic variation and genetic linkage of avian microsatellites in a new turkey population for genetic mapping

Efforts to build a comprehensive genetic linkage map for the turkey (Meleagris gallopavo) have focused on development of genetic markers and experimental resource families. In this study, PCR amplification was attempted for 772 microsatellite markers that had been previously developed for three avian species (chicken, quail and turkey). Allelic polymorphism at 410 markers (53.1% of total examined) was determined by genotyping ten individuals (six F1 parents and four grandparents) in a new resource population specifically developed for genetic linkage mapping. Of these 410 markers, 109 (26.6%) were polymorphic in the tested individuals, with an average of 2.3 alleles per marker. Higher levels of polymorphism were found for the turkey-specific markers (61.1%) than for the chicken (22.7%) or quail-specific markers (33.3%). To test the fidelity of the matings, demonstrate the power of these families for linkage analysis, and determine genetic linkage relationships, 86 polymorphic markers were genotyped for up to 224 birds including founder grandparents, parents and F2 progeny. Linkage relationships for many of the chicken markers elucidated in the turkey were comparable to those observed in the chicken. These data demonstrate that the new UMN/NTBF resource population will provide a solid foundation for constructing a comparative genetic map of the turkey.

Detection and parameter estimation for quantitative trait loci using regression models and multiple markers

A strategy of multi-step minimal conditional regression analysis has been developed to determine the existence of statistical testing and parameter estimation for a quantitative trait locus (QTL) that are unaffected by linked QTLs. The estimation of marker-QTL recombination frequency needs to consider only three cases: 1) the chromosome has only one QTL, 2) one side of the target QTL has one or more QTLs, and 3) either side of the target QTL has one or more QTLs. Analytical formula was derived to estimate marker-QTL recombination frequency for each of the three cases. The formula involves two flanking markers for case 1), two flanking markers plus a conditional marker for case 2), and two flanking markers plus two conditional markers for case 3). Each QTL variance and effect, and the total QTL variance were also estimated using analytical formulae. Simulation data show that the formulae for estimating marker-QTL recombination frequency could be a useful statistical tool for fine QTL mapping. With 1 000 observations, a QTL could be mapped to a narrow chromosome region of 1.5 cM if no linked QTL is present, and to a 2.8 cM chromosome region if either side of the target QTL has at least one linked QTL.

Detection of Quantitative Trait Loci Affecting Milk Production, Health, and Reproductive Traits in Holstein Cattle

We report putative quantitative trait loci affecting female fertility and milk production traits using the merged data from two research groups that conducted independent genome scans in Dairy Bull DNA Repository grandsire families to identify quantitative trait loci (QTL) affecting economically important traits. Six families used by both groups had been genotyped for 367 microsatellite markers covering 2713.5 cM of the cattle genome (90%), with an average spacing of 7.4 cM. Phenotypic traits included PTA for pregnancy rate and daughter deviations for milk, protein and fat yields, protein and fat percentages, somatic cell score, and productive life. Analysis of the merged dataset identified putative quantitative trait loci that were not detected in the separate studies, and the pregnancy rate PTA estimates that recently became available allowed detection of pregnancy rate QTL for the first time. Sixty-one putative significant marker effects were identified within families, and 13 were identified across families. Highly significant effects were found on chromosome 3 affecting fat percentage and protein yield, on chromosome 6 affecting protein and fat percentages, on chromosome 14 affecting fat percentage, on chromosome 18 affecting pregnancy rate, and on chromosome 20 affecting protein percentage. Within-family analysis detected putative QTL associated with pregnancy rate on six chromosomes, with the effect on chromosome 18 being the most significant statistically. These findings may help identify the most useful markers available for QTL detection and, eventually, for marker-assisted selection for improvement of these economically important traits.

A software tool for the graphical visualization of large and complex populations

Pedigree drawing is an essential tool in genetic and genealogical studies. A genetic analysis of a disease often starts with drawing a pedigree to show the overall population structure, the relationship among individuals, and gene flows from generation to generation. Such a graphical presentation of a pedigree is valuable for understanding the nature of the disease such as inheritance mode and familial trends and of the disease, for tracing the source of a detrimental gene, and for identifying the founders of the population. A genealogy study typically requires a pedigree drawing to show relationships among individuals. However, as the size and complexity of a pedigree increase, drawing a clear pedigree becomes a challenge. Pedigraph has been created to solve this problem. We developed a software tool named Pedigraph capable of creating artistic graphical pedigree drawings of large and complex populations with flexible options for pedigree analysis. Preliminary tests show that this software tool has great potential to be a useful tool for research in breeding, genetics and genomics in plants, animals, and zoo species, as well as a useful tool for studying history and human populations.

Statistical analysis and experimental design for mapping genes of complex traits in domestic animals

Gene mapping for complex traits has been an active and challenging research area in humans, agricultural and laboratory species. In domestic animals and poultry, the goal of gene mapping is to find genes associated with production, reproduction and health traits. In humans, gene mapping for complex diseases has been a significant challenge. Although gene mapping results are accumulating rapidly, confirmed results are scarce. In domestic animals, gene mapping results are also accumulating rapidly and several large scale mapping projects in chickens, swine, and dairy cattle are currently in progress. Although new and more reliable results can be expected in the near future, the precise understanding of genes underlying complex traits is still some distance away. Recently, the approach of mapping genes of complex traits has been applied to gene expression data to map transcription regulatory elements, adding a new dimension to gene mapping. Statistical analysis of gene mapping data and experimental design are critical components of a gene mapping research. The purpose of this review article is to discuss current status as well as future directions and challenges in statistical analysis and experimental design for mapping genes of complex traits in domestic animals and poultry.

Linkage analysis using direct and indirect counting and relative efficiencies for codominant and dominant loci

A method based on direct and indirect counting is developed for rapid and accurate linkage analysis for codominant and dominant loci. Methods for estimating gender-specific recombination frequencies are available for cases where at least one of the two loci is multiallelic and for biallelic loci with mixed parental linkage phases where at least one locus is codominant. Most of the estimates of gender-average and gender-specific recombination frequencies required iterative solutions. The new method makes use of the full data set, yields exact estimates of the recombination frequencies when the observed and expected genotypic frequencies are equal, and are computationally efficient. Relative efficiency of various data types is affected by the inheritance mode and by parental linkage phases of biallelic loci, but unaffected by the locus polymorphism when using the full data set for linkage analysis. The ability to determine parental linkage phases is affected by the locus polymorphism as well as inheritance mode. Intercross (or F-2 design) is more efficient for mapping codominant loci, whereas backcross is more efficient if dominance is involved. Mixed parental linkage phases of biallelic loci are less efficient than coupling or repulsion linkage phases. Ignoring noninformative offspring results in biased estimates of recombination frequency for biallelic loci only and reduced LOD scores for all cases.

Generation and exploration of a dense genetic map in a region of a QTL affecting corpora lutea in a Meishan x Yorkshire cross

Previously genomic scans revealed quantitative trait loci (QTL) on porcine Chromosome 8 (SSC8) as significantly affecting the number of corpora lutea (CL) in swine. In one study, statistical evidence for the putative QTL was found in the chromosomal region defined by the microsatellites (MS) SW205, SW444, SW206, and SW29. A Yeast Artificial Chromosome library was screened by using the corresponding primers for clones containing these MS by PCR. From five positive YAC clones, 10 additional MS were isolated and mapped to SSC8 with the INRA-University of Minnesota porcine Radiation Hybrid (IMpRH) panel. The genetic map position of the QTL has been refined by addition of these 10 markers. The QTL evaluation included pedigrees of F2-intercross Meishan x Yorkshire design, with phenotypic data of 108 F2 female offspring and genotypic data for 29 MS markers on SSC8. The analysis was performed by using the least squares regression method. The calculated QTL effect for CL obtained by the multilocus least squares method showed a maximum test statistic (F value = 13.98) at position 99 cM between three MS derived from YACs containing SW205 and SW1843 spanning an interval of 7.1 cM. The point-wise (nominal) P-value was 5.21 x 10-6 corresponding to a genome-wide P-value of 0.009. The additive QTL effect explained 17.4% of the phenotypic variance.

A new mutation in the connexin 32 gene was found in Charcot- Marie-Tooth disease in Chinese patients

OBJECTIVE

To investigate the characteristics of gene mutations of connexin 32 exon 2 in Charcot-Marie-Tooth disease in Chinese patients.

METHODS

Screening for connexin 32 gene mutation was conducted in 6 unrelated CMT1 patients without duplication and 10 unrelated CMT2 patients. Mobility shift of exon 2 was analyzed by SSCP and further confirmed by sequencing. The PCR products were cut by appropriate restricted enzyme in 50 normal controls.

RESULTS

One missense mutation at nucleotides 62(G-->A) was found in a CMT1 patient. 50 normal controls were analyzed by the enzyme HaeIII and no abnormality was found. This proved that the mutation was the cause of disease.

CONCLUSION

This mutation has not been reported previously. A proportion of CMTX patients may exist in the group of CMT1 patients in China.

Standardization and conversion of marker polymorphism measures

Large scale gene mapping efforts in domestic animals have generated and mapped a large number of genetic markers that are useful for mapping quantitative trait and disease loci and for DNA diagnostic purposes such as parentage testing. Marker polymorphism is an important criterion for selecting genetic markers in planning experiment for mapping quantitative trait loci or for DNA diagnostic purposes. Current formulations of marker polymorphism measures are functions of marker allele frequencies. In this study, two measures of marker polymorphism that are available from gene mapping studies and do not require allele frequencies were proposed and analyzed: the observed polymorphic information content (PIC) and the observed family information content (FIC). The observed FIC was more stable than the observed PIC because the observed FIC is unaffected by the variation in the frequency of heterozygous parents. However, both FIC and PIC are dependent on the gene mapping design. The effective number of alleles is recommended as a tool to standardize marker polymorphism measures so that polymorphism of different markers can be compared on an equal basis, and to obtain a new polymorphism measure (such an exclusion probability) from an existing measure (such as FIC). The usage of the effective number of alleles to standardize FIC, PIC and exclusion probabilities is illustrated using genetic markers in a published linkage map.

A genome scan for QTL influencing milk production and health traits in dairy cattle

A genome scan was conducted in the North American Holstein-Friesian population for quantitative trait loci (QTL) affecting production and health traits using the granddaughter design. Resource families consisted of 1,068 sons of eight elite sires. Genome coverage was estimated to be 2,551 cM (85%) for 174 genotyped markers. Each marker was tested for effects on milk yield, fat yield, protein yield, fat percentage, protein percentage, somatic cell score, and productive herd life using analysis of variance. Joint analysis of all families identified marker effects on 11 chromosomes that exceeded the genomewide, suggestive, or nominal significance threshold for QTL effects. Large marker effects on fat percentage were found on chromosomes 3 and 14, and multimarker regression analysis was used to refine the position of these QTL. Half-sibling families from Israeli Holstein dairy herds were used in a daughter design to confirm the presence of the QTL for fat percentage on chromosome 14. The QTL identified in this study may be useful for marker-assisted selection and for selection of a refined set of candidate genes affecting these traits.

Designs of reference families for the construction of genetic linkage maps

The reference family panel is the foundation of a gene mapping program because it affects the cost and quality of the genetic linkage maps, and should be designed to yield reliable linkage detection and locus ordering at minimal gene mapping cost. A map cost function was defined as the number of genotypes required per marker per unit of genome coverage and was used to obtain optimal designs with respect to linkage detection. An ordering reliability function was defined as the likelihood ratio of the most likely order to the second most likely order of genetic markers and was used to find optimal designs with respect to locus ordering. Optimum levels of recombination frequency were found to be in the neighborhood of 0.11-0.15 for linkage detection and were in the region of 0.05-0.20 for locus ordering. Therefore, recombination frequencies optimal for linkage detection are also optimal for locus ordering. Based on the optimal detection levels, sample size (number of offspring) and map cost requirements were derived for six representative designs, assuming gender-specific linkage maps and two alleles with equal frequency for each marker. The sample size required for linkage detection ranged from 168 to 432 offspring for full-sib designs and ranged from 350 to 600 offspring for half-sib designs depending on the family size and the target LOD score, with corresponding minimal map costs of 10-20 genotypes per marker per centiMorgan map coverage. Locus ordering generally requires more genotypes than linkage detection. For full-sib designs, meioses from both genders should be used for locus ordering even when the maps are gender-specific. For half-sib designs, additional families may be needed for locus ordering. Sample size for ordering closely linked loci as required by positional cloning were provided. Effects of family size, grandparents, and marker polymorphism on design efficiency were analyzed.

Detection of putative loci affecting milk production and composition, health, and type traits in a United States Holstein population

Quantitative trait loci affecting milk yield and composition, health, and type traits were studied for seven large grandsire families of US Holstein using the granddaughter design. The families were genotyped at 20 microsatellite markers on 15 chromosomes, and the effects of the marker alleles were analyzed for 28 traits (21 type traits, 5 milk yield and composition traits, somatic cell score, and productive herd life). Markers BM415 on chromosome 6 and BM6425 on chromosome 14 were associated with effects on protein percentage in a single grandsire family. The latter marker had a lower probability of being associated with changes in milk yield and fat percentage in the same family. Increases in productive herd life were associated with an allele at marker BM719 on chromosome 16 in one grandsire family.

Reduced bovine leukaemia virus proviral load in genetically resistant cattle

The bovine leukaemia virus (BLV) is an exogenous retrovirus that is closely related to the human T cell leukaemia viruses. Genetic resistance and susceptibility to persistent lymphocytosis (PL), an advanced subclinical stage of infection characterized by a polyclonal expansion of the infected B cell population, have been mapped to structural motifs in bovine MHC DRB3 (class II) alleles. To determine whether alleles of DRB3 influence the number of BLV-infected B cells in peripheral blood, seven pairs of Holstein cows naturally infected with BLV were matched on the basis of DRB3 genotype (resistance or susceptibility to PL), age, and year of seroconversion. Flow cytometry was used to separate B cell populations that then were tested for the presence of provirus by a single-cell PCR methodology. Animals with the PL-resistance associated DRB3.2*11 allele had significantly fewer BLV-infected B cells than did age- and seroconversion-matched cows with DRB3 alleles associated with susceptibility to PL. Our results demonstrate that DRB3 or a closely linked gene may play a direct role in controlling the number of BLV-infected peripheral B cells in vivo. Association of MHC class II alleles with resistance to disease progression in cattle naturally infected with BLV provides a unique immunogenetic model for the study of host resistance to human and other animal retroviral infections.

Detection of putative loci affecting conformational type traits in an elite population of United States Holsteins using microsatellite markers

Quantitative trait loci affecting conformational type traits were studied in seven large grandsire families of US Holsteins using the granddaughter design and 16 microsatellite markers on 10 chromosomes. The most significant marker effect was marker BM203 (chromosome 27) for dairy form in a single grandsire family. A multivariate analysis for dairy form and milk yield was also conducted, and the result was highly significant, indicating that a segregating quantitative trait locus or loci affecting dairy form and milk yield could exist near BM203 on chromosome 27. Marker BM1258 (chromosome 23) had a significant effect on udder depth. A multivariate analysis on udder depth and somatic cell score was conducted for markers 513 and BM1258, and both markers showed significant effects on these two traits, indicating that one or several quantitative trait loci affecting udder depth and mastitis might exist on chromosome 23. Marker BM4204 (chromosome 9) had a significant effect on foot angle and on the composite index of traits pertaining to feet and legs, indicating that one or several quantitative trait loci affecting traits pertaining to feet and legs might exist on chromosome 9. Selection on these markers could increase genetic progress within these families.

Exclusion probabilities of 22 bovine microsatellite markers in fluorescent multiplexes for semiautomated parentage testing

Six multiplexes developed for semiautomated fluorescence genotyping were evaluated for parentage testing. These multiplexes contained primer pairs for the amplification of 22 microsatellites on 17 bovine autosomes. Exclusion probabilities were determined from genotypes of 1022 Holstein cattle and 311 beef cattle belonging to five breeds. Two cases were considered: case 1, genotypes are known for an alleged parent and an offspring but genotypes of a confirmed parent are unknown; and case 2, genotypes are known for an alleged parent, a confirmed parent and an offspring. If the alleged parent is not the true parent, then the 22 markers will exclude the alleged parent with a probability of > 0.9986 for case 1 and with a probability of > 0.99999 for case 2. On the basis of these exclusion probabilities, the probability that an alleged parent will be falsely included as a parent is in the range of 1/716 to 1/2845 for case 1 and 1/1.2 million to 1/159753 for case 2. In addition to these results, a rapid and efficient non-organic method for extraction of DNA from semen is described.

A male linkage map of the cattle (Bos taurus) genome

A male linkage map of the cattle (Bos taurus) genome was constructed using nine large half-sib families. The map consists of 269 loci, of which 249 are microsatellites and 20 are structural genes. Among the 249 microsatellites, 140 are markers selected from other maps and 98 are new assignments. Chromosome assignment were established for 35 new markers by somatic cell hybrid analysis, of which 26 were confirmed by linkage analysis. Genome coverage is 1975 cM contained within terminal markers on all 29 autosomes. The average distance between adjacent loci is 9.7 cM, with 72.1% of the map intervals < or = 15 cM and 4.9% of the intervals > or = 25 cM. The inclusion of mapped markers permitted integration and comparisons with other maps, facilitating the identification of discrepancies in chromosome assignment, gene order, and map distance. The inclusion of Type I and blood group markers in the map was useful for comparative mapping, revealing possible blood group orthologies between humans and cattle. The map generated will serve as a useful tool for comparative mapping, mapping of quantitative trait loci and marker assisted selection.

The prevalence of proviral bovine leukemia virus in peripheral blood mononuclear cells at two subclinical stages of infection

The bovine leukemia virus (BLV) is an oncogenic retrovirus that is associated with the development of persistent lymphocytosis (PL) and lymphoma in cattle. While B lymphocytes have been shown to be the primary cellular target of BLV, recent studies suggest that some T lymphocytes and monocytes may be infected by the virus. Because virally altered functions of monocytes and/or T cells could contribute to the development of lymphoproliferative disease, we sought to clarify the distribution of the BLV provirus in subpopulations of peripheral blood mononuclear cells in seropositive cows with and without PL. CD2+ T cells, monocytes, and CD5+ and CD5- B cells were sorted by flow cytometry and tested for the presence of BLV by single-cell PCR. We did not obtain convincing evidence that peripheral blood monocytes or T lymphocytes contain the BLV provirus in seropositive cows with or without PL. In seropositive cows without PL (n=14), BLV-infected CD5+ and CD5- B cells accounted for 9.2% +/- 19% and 0.1% +/- 1.8% of circulating B lymphocytes, respectively. In cows with PL (n=5), BLV-infected CD5+ and CD5- B cells accounted for 66% +/- 4.8% and 13.9% +/- 6.6% of circulating B lymphocytes, respectively. The increase in lymphocyte numbers in cows with PL was entirely attributable to the 45-fold and 99-fold expansions of infected CD5+ and CD5- B-cell populations, respectively. Our results demonstrate that B cells are the only mononuclear cells in peripheral blood that are significantly infected with BLV. On the basis of the absolute numbers of infected cells in seropositive, hematologically normal animals, there appear to be differences in susceptibility to viral spread in vivo that may be under the genetic control of the host.

Prediction of informativeness for microsatellite markers among progeny of sires used for detection of economic trait loci in dairy cattle

Individual loci affecting economic traits can be located using genetic linkage. Application of either daughter or granddaughter design requires determination of allele origin in the progeny. If only the sires and their progeny are genotyped, the paternal allele origin of progeny having the same genotype as the sire cannot be determined. The expected frequency of informative sons can be predicted for each sire and genetic marker from the allele frequencies in the population. The accuracy of a predictor of the frequency of informative progeny was tested on 103 grandsire x microsatellite combinations. Number of sons per grandsire varied from 24 to 129. Allele frequencies in the population were estimated by genotyping seven sires. The regression of the frequency of informative sons on the predicted frequency was 1.04 with a zero intercept model. Thus, considering the large number of genetic markers available for analysis, predicted informative frequency is a useful criterion for selection of genetic markers.

A panel of polymorphic bovine, ovine and caprine microsatellite markers

A panel of 81 new polymorphic bovine microsatellite markers is described, together with further information on a previously reported group of 16 markers. The mean polymorphism information content of the 97 markers determined in 20 cattle was 0.66. Seventy-three of these markers have been assigned to chromosomes by either linkage analysis or use of hybrid cell panels. Thirty-nine of the markers were polymorphic in sheep, and 32 were polymorphic in goat. This study identified a set of 18 robust markers that were polymorphic in all three species and that covered 14 bovine chromosomes. This provides a single group of markers, which would be suited to genetic distance analysis and parentage control in cattle, sheep and goat.

Genetic mapping of F13A to BTA23 by sperm typing: difference in recombination rate between bulls in the DYA-PRL interval

Our objective was to extend the linkage and comparative maps of BTA23 by determining whether the structural gene for the A subunit of blood coagulation factor XIII (F13A) is linked to BoLA-DYA, a centromeric marker, or the distally located gene encoding prolactin (PRL). Bovine F13A was mapped relative to DYA and PRL in an experiment that examined segregation of alleles in 176 sperm. Genotyping was performed by PCR-RFLP for all loci, following amplification of the haploid genome by primer extension preamplification. F13A was found to be linked to PRL (theta = 0.314 +/- 0.038). The most likely order is DYA-PRL-F13A (odds > 10(4):1). This result demonstrates conservation of synteny between BTA23 and most of HSA6p. Surprisingly, theta DYA-PRL was 0.310 +/- 0.039, 83.4% greater (P < 0.02) than we found for another bull (Van Eijk et al., Mamm. Genome 4: 113, 1993). The difference in recombination rate in the DYA-PRL interval provides further evidence for an unusual recombination hot spot between the bovine Mhc class IIa and class IIb subregions and suggests that bull-specific maps may be necessary for marker-assisted selection.

Gene-centromere mapping of bovine DYA, DRB3, and PRL using secondary oocytes and first polar bodies: evidence for four-strand double crossovers between DYA and DRB3

A genetic map consisting of three loci anchored at the centromere of bovine chromosome 23 was constructed by genotyping secondary oocytes (SO) and first polar bodies (PB1) using a polymerase chain reaction (PCR)-based approach. Primary oocytes arrested in prophase of meiosis I were stimulated in vitro to resume division and extrude PB1. Sixty SO and their matched PB1 were collected from 14 cows by micro-manipulation, subjected to amplification of the whole genome by primer extension preamplification, and genotyped independently for the linked genes PRL, DRB3, and DYA by PCR-RFLP analysis. Single-locus analysis of gene-centromere recombination rates were theta cen-DYA = 0.042, theta cen-DRB3 = 0.113, and theta cen-PRL = 0.166. The most likely order is cen-DYA-DRB3-PRL. Analysis of typing data from 3 cows revealed three meiotic divisions consistent with "linkage phase exchange" between DYA and DRB3 or PRL. One of the three linkage phase exchanges was confirmed by complementary genotypes in a matched secondary oocyte-first polar body pair. Such linkage phase exchanges could result from four-strand double crossovers between homologous chromosomes. Because all four gametes produced by four-strand double crossovers will be recombinant, more frequent occurrence of such events in females may explain the sexual dimorphism in genetic maps. Alternatively, four-strand crossovers could represent a type of recombination hotspot between DYA and DRB3, suggesting a mechanism for the high recombination frequency (15%) between these two class II genes of the bovine major histocompatibility complex.

Linkage information content and efficiency of full-sib and half-sib designs for gene mapping

The accuracy of a genetic map depends on the amount of linkage information contained in the data set used for construction of the map. The amount of linkage information is related to the designs employed for linkage analysis. The purpose of this study was to provide general formulations for various genotyping schemes and family structures in order to evaluate the amount of linkage information in a data set. Linkage information content (LIC) was defined as the frequency of fully informative gametes, which are gametes from doubly heterozygous parents with known linkage phases. Depending on the design, LIC is based on two generations if the parental phases are determined statistically, or three generations if the parental phases are determined genetically. Different schemes were considered in deriving LIC: (1) genotyping of one parent or two parents, and (2) genotyping of two or three generation families. The LIC for a full-sib design was found to be generally greater than for a half-sib design but requires typing a large number of individuals when at least one locus has only two alleles. The efficiency of the full-sib design is reduced significantly if a sex-specific linkage map is sought.

Multilocus analysis for gene-centromere mapping using first polar bodies and secondary oocytes

Polar body and oocyte typing is a new technique for gene-centromere mapping and for generating female linkage maps. A maximum likelihood approach is presented for ordering multiple markers relative to the centromere and for estimating recombination frequencies between markers and between the centromere and marker loci. Three marker-centromere orders are possible for each pair of markers: two orders when the centromere flanks the two markers and one order when the centromere is flanked by the two markers. For each possible order, the likelihood was expressed as a function of recombination frequencies for two adjacent intervals. LOD score for recombination frequency between markers or between the centromere and a marker locus was derived based on the likelihood for each gene-centromere order. The methods developed herein provide a general solution to the problem of multilocus gene-centromere mapping that involves all theoretical crossover possibilities, including four-strand double crossovers.

A genetic map of nine loci on bovine chromosome 2

A male-specific genetic linkage map of nine loci on bovine Chromosome (Chr) 2 (BTA2) was constructed from 306 offspring belonging to six paternal half-sib families. Loci studied were the structural genes for liver/bone/kidney alkaline phosphatase (ALPL). Gardner-Rasheed feline sarcoma (v-fgr) oncogene homolog (FGR), alpha-L-fucosidase 1 (FUCA1), and fibronectin 1 (FN1), and the microsatellite loci ARO28, DU17S2, DU17S3, DU17S4, and DU17S5. Genotyping was performed by restriction fragment length polymorphism (RFLP) for structural genes and polymerase chain reaction (PCR) for the microsatellites. Two genetically independent linkage groups were identified. The order of genes in the first linkage group, L31, is (ARO28-FN1)-FGR-FUCA1-ALPL, covering a map distance of 34.1 cM between terminal markers. The second linkage group, L32, consists of DU17S2-DU17S5-DU17S4-DU17S3 and is 41.3 cM in length. Genetic linkage between FN1 and FGR confirms previous physical assignment of these genes to the same synteny group. Currently, the genetic linkage of FN1 and FGR is unique to cattle and thus localizes a site of chromosomal evolution to a 22-cM interval between the two loci.

The bovine B and C blood group systems are not likely to be the orthologues of human RH: an interesting twist in the comparative map

We tested the hypothesis that either the bovine B or C blood group system is the orthologue of human RH. A comparative linkage mapping strategy was applied, using blood typing and restriction fragment length polymorphism (RFLP) analysis of four loci linked to RH on HSA1; PGD, FGR, ALPL and FUCA1. Four sires with a total of 255 half-sib offspring were used for the linkage analysis. Strong support for linkage between ALPL, FUCA1 and FGR was obtained for all sire families (lod scores > 11 for all pairwise comparisons). This new linkage group was assigned to bovine synteny group U17 based on previous somatic cell mapping of the FGR locus. The most favoured order is ALPL-FUCA1-FGR (2.18:1), with ALPL and FGR 5.4 CM and 6.3 CM, respectively, from FUCA1. The B and C blood group systems and PGD were genetically independent of each other and all other markers, indicating that neither B nor C is likely to be the bovine orthologue of human RH. However, given available comparative mapping data, there is some chance that the bovine orthologue of RH is on bovine synteny group U6. Although gene order appears to be conserved with humans, the differences in recombination rates between these three loci in cattle, humans and mice strongly suggest that it is not possible to use human map distances to predict map distances in cattle, making it imperative that bovine gene mappers continue to emphasize adding type I markers to the bovine linkage map.

Milk and fat yields decline in bovine leukemia virus-infected Holstein cattle with persistent lymphocytosis

Effects of bovine leukemia virus (BLV) infection on milk and fat yields were studied by using data collected from Holstein cows over a 6-year period. Milk and fat yields in BLV-infected cows with persistent lymphocytosis (PL) declined significantly relative to their BLV-infected non-PL herdmates. Declines were most pronounced in cows older than 6 years. The estimated loss to the dairy industry due to PL is more than $42 million annually. A major histocompatibility complex class I (BoLA-A) allele that has been previously associated with resistance to PL was associated with longevity and realization of milk production potentials, indicating that genetic resistance to PL will have an economic benefit in herds where BLV is endemic.

Detection of bovine leukemia virus proviral DNA in individual cells

We have developed a method of analyzing individual cells to detect proviral DNA of the bovine leukemia virus (BLV) using flow cytometry and PCR. Individual cells of the BL3* cell line, which contain multiple integrated copies of the BLV provirus, and the uninfected cell line BL3(0), were sorted into wells of a 96-well plate. Following cell lysis, portions of the BLV envelope (ENV) and cellular prolactin (PRL) genes were amplified simultaneously using PCR. Viral and cellular products of first-round PCR were amplified separately in a second round of PCR using "heminested" primers. Separation of the PCR products by polyacrylamide gel electrophoresis yielded distinct fragments of the predicted sizes. The operational sensitivity of this method for the detection of virus was > 90% when testing single infected cells. In addition, we were able to reliably amplify DNA from a single BL3* cell among as many as 10(5) BL3(0) cells and established that the sensitivity for detecting a single infected cell among 20, 100, or 1000 uninfected cells was at least 90%. Estimates of low percentages of infected cells were obtained by applying probability theory to results of experiments conducted on wells containing more than one cell. Using these methods, B lymphocytes obtained from the peripheral blood of BLV-infected cattle were tested for proviral DNA. BLV ENV was identified in 76.9 +/- 4.9% of single B cells tested from a seropositive animal with persistent lymphocytosis (PL), but in only 0.033 +/- 0.009% of B cells from another seropositive cow without PL.(ABSTRACT TRUNCATED AT 250 WORDS)

Estimation of genetic parameters for somatic cell score in Holsteins

Genetic parameters of somatic cell scores for Holstein cows were estimated using an animal model and REML for two data sets. Set 1, with 13,017 records from 5278 cows, was used to obtain variance components, heritability, and repeatability for two lactation measures: the simple average and the weighted average of test day data. Set 2, with 14,418 records from 4806 cows, was used to obtain genetic correlations for the simple average between lactations 1 and 2, between lactations 1 and 3, and between lactations 2 and 3. Simple and weighted average of test day somatic cell scores had the same heritabilities (.12) and repeatabilities (.35). Phenotypic variances were about 1.2, and herd-sire interaction variances were small (.002). Genetic correlation for somatic cell score was about .55 between lactations 1 and 2 and between lactations 1 and 3 and .65 between lactations 2 and 3. Phenotypic correlation was .20 between lactations 1 and 2, .16 between lactations 1 and 3, and .31 between lactations 2 and 3.

Multitrait animal model with genetic groups

Genetic groups of unknown parents are extended to the multitrait animal model. Computationally feasible mixed model equations are obtained. A strategy to include genetic groups for missing data is proposed. Canonical and triangular transformations can be applied to the multitrait animal model with groups if the transformations can be applied to the same model without groups. Formulations for REML estimation with groups are derived, and the results are almost as feasible as REML estimation without groups. A numerical example is given to illustrate computations of REML formulations.

Prediction error variance and restricted maximum likelihood estimation for animal model with relationship grouping

Genetic evaluation using animal model with relationship grouping has been shown to be feasible. However, algorithms were unavailable for prediction error variance and REML estimation of variance components. This paper shows that prediction error variance of an estimable function of the total merit of additive genetic and group effects is a simple function of a generalized inverse of the coefficient matrix for a transformed mixed model equation or of the inverse of the coefficient matrix when it is restricted to full rank. The REML algorithms, using the transformed equation, having slightly more complicated expressions than usual but could be more feasible computationally. Formulae for prediction error variance apply in general. The REML algorithms are extended to an animal model with an arbitrary number of random factors and can be extended to estimate covariance components.