Genomic imputation and evaluation using high-density Holstein genotypes

Genomic evaluations for 161,341 Holsteins were computed by using 311,725 of 777,962 markers on the Illumina BovineHD Genotyping BeadChip (HD). Initial edits with 1,741 HD genotypes from 5 breeds revealed that 636,967 markers were usable but that half were redundant. Holstein genotypes were from 1,510 animals with HD markers, 82,358 animals with 45,187 (50K) markers, 1,797 animals with 8,031 (8K) markers, 20,177 animals with 6,836 (6K) markers, 52,270 animals with 2,683 (3K) markers, and 3,229 nongenotyped dams (0K) with >90% of haplotypes imputable because they had 4 or more genotyped progeny. The Holstein HD genotypes were from 1,142 US, Canadian, British, and Italian sires, 196 other sires, 138 cows in a US Department of Agriculture research herd (Beltsville, MD), and 34 other females. Percentages of correctly imputed genotypes were tested by applying the programs findhap and FImpute to a simulated chromosome for an earlier population that had only 1,112 animals with HD genotypes and none with 8K genotypes. For each chip, 1% of the genotypes were missing and 0.02% were incorrect initially. After imputation of missing markers with findhap, percentages of genotypes correct were 99.9% from HD, 99.0% from 50K, 94.6% from 6K, 90.5% from 3K, and 93.5% from 0K. With FImpute, 99.96% were correct from HD, 99.3% from 50K, 94.7% from 6K, 91.1% from 3K, and 95.1% from 0K genotypes. Accuracy for the 3K and 6K genotypes further improved by approximately 2 percentage points if imputed first to 50K and then to HD instead of imputing all genotypes directly to HD. Evaluations were tested by using imputed actual genotypes and August 2008 phenotypes to predict deregressed evaluations of US bulls proven after August 2008. For 28 traits tested, the estimated genomic reliability averaged 61.1% when using 311,725 markers vs. 60.7% when using 45,187 markers vs. 29.6% from the traditional parent average. Squared correlations with future data were slightly greater for 16 traits and slightly less for 12 with HD than with 50K evaluations. The observed 0.4 percentage point average increase in reliability was less favorable than the 0.9 expected from simulation but was similar to actual gains from other HD studies. The largest HD and 50K marker effects were often located at very similar positions. The single-breed evaluation tested here and previous single-breed or multibreed evaluations have not produced large gains. Increasing the number of HD genotypes used for imputation above 1,074 did not improve the reliability of Holstein genomic evaluations.

Safety and pharmacokinetics of an intramuscular humanized monoclonal antibody to respiratory syncytial virus in premature infants and infants with bronchopulmonary dysplasia. The MEDI-493 Study Group

BACKGROUND

Respiratory syncytial virus (RSV) is the leading cause of lower respiratory disease in infants and children. MEDI-493 (palivizumab, Synagis) is a humanized monoclonal IgG1 antibody to the fusion protein of RSV, and it is highly active in vitro against RSV A and B strains.

OBJECTIVE

To describe the safety, tolerance, immunogenicity and pharmacokinetics of monthly intramuscular injections of MEDI-493 among premature infants and children with bronchopulmonary dysplasia and to compare these data with information previously obtained with intravenous dosing.

DESIGN

A Phase I/II multicenter, open label, escalating dose clinical trial. PATIENT POPULATION AND DOSING REGIMEN: Children (n=65) born prematurely at < or =35 weeks of gestation who were < or =6 months of age (n=41) and children with bronchopulmonary dysplasia who were < or =24 months of age (n=24) were enrolled. From 1 to 5 monthly injections were given at doses of 5 mg/kg (n=11), 10 mg/kg (n=6) and 15 mg/kg (n=48). Serum was collected before administration of each dose, 30 days after the last dose, and 2, 7 and 14 days after the first and second doses for measurement of MEDI-493 concentrations by enzyme-linked immunosorbent assay.

RESULTS

The pharmacokinetics of MEDI-493 were similar to those of other human IgG1 antibodies. Mean serum MEDI-493 concentrations were 91.1 microg/ml (range, 52.3 to 174.0) 2 days after the initial dose of 15 mg/kg and 49.2 microg/ml (range, 13.5 to 132.0) at 30 days. Monthly dosing of 15 mg/kg maintained mean trough concentrations of approximately 70 microg/ml. These concentrations were similar to previously published trough concentrations after i.v. administration. MEDI-493 injections were well-tolerated. Only three children had adverse events judged to be possibly related to MEDI-493. Ten children had transient, low titer anti-MEDI-493 binding titers (1:10 to 1:40) which were not associated with a pattern of specific adverse events or alterations of MEDI-493 concentrations. Two patients in the 5-mg/kg dose group were hospitalized for RSV; no RSV hospitalizations were found in the higher dose groups.

CONCLUSIONS

MEDI-493 was safe and well-tolerated. Monthly intramuscular doses of 15 mg/kg maintained mean trough serum concentrations that were above 40 microg/ml (the value associated with 99% reduction of pulmonary RSV in the cotton rat model). These concentrations were similar to those previously reported with i.v. administration of MEDI-493.

Heritability of rectal temperature and genetic correlations with production and reproduction traits in dairy cattle

Genetic selection for body temperature during heat stress might be a useful approach to reduce the magnitude of heat stress effects on production and reproduction. Objectives of the study were to estimate the genetic parameters of rectal temperature (RT) in dairy cows in freestall barns under heat stress conditions and to determine the genetic and phenotypic correlations of rectal temperature with other traits. Afternoon RT were measured in a total of 1,695 lactating Holstein cows sired by 509 bulls during the summer in North Florida. Genetic parameters were estimated with Gibbs sampling, and best linear unbiased predictions of breeding values were predicted using an animal model. The heritability of RT was estimated to be 0.17 ± 0.13. Predicted transmitting abilities for rectal temperature changed 0.0068 ± 0.0020°C/yr from (birth year) 2002 to 2008. Approximate genetic correlations between RT and 305-d milk, fat, and protein yields, productive life, and net merit were significant and positive, whereas approximate genetic correlations between RT and somatic cell count score and daughter pregnancy rate were significant and negative. Rectal temperature during heat stress has moderate heritability, but genetic correlations with economically important traits mean that selection for RT could lead to lower productivity unless methods are used to identify genes affecting RT that do not adversely affect other traits of economic importance.

Explorations in genome-wide association studies and network analyses with dairy cattle fertility traits

The objective of this study was to identify single nucleotide polymorphisms and gene networks associated with 3 fertility traits in dairy cattle-daughter pregnancy rate, heifer conception rate, and cow conception rate-using different approaches. Deregressed predicted transmitting abilities were available for approximately 24,000 Holstein bulls and 36,000 Holstein cows sampled from the National Dairy Database with high-density genotypes. Of those, 1,732 bulls and 375 cows had been genotyped with the Illumina BovineHD Genotyping BeadChip (Illumina Inc., San Diego, CA). The remaining animals were genotyped with various chips of lower density that were imputed to high density. Univariate and trivariate genome-wide association studies (GWAS) with both medium- (60,671 markers) and high-density (312,614 markers) panels were performed for daughter pregnancy rate, heifer conception rate, and cow conception rate using GEMMA (version 0.94; http://www.xzlab.org/software.html). Analyses were conducted using bulls only, cows only, and a sample of both bulls and cows. The partial correlation and information theory algorithm was used to develop gene interaction networks. The most significant markers were further investigated to identify putatively associated genes. Little overlap in associated genes could be found between GWAS using different reference populations of bulls only, cows only, and combined bulls and cows. The partial correlation and information theory algorithm was able to identify several genes that were not identified by ordinary GWAS. The results obtained herein will aid in further dissecting the complex biology underlying fertility traits in dairy cattle, while also providing insight into the nuances of GWAS.

Genetic evaluation of lactation persistency for five breeds of dairy cattle

Cows with high lactation persistency tend to produce less milk than expected at the beginning of lactation and more than expected at the end. Best prediction of lactation persistency is calculated as a function of trait-specific standard lactation curves and linear regressions of test-day deviations on days in milk. Because regression coefficients are deviations from a tipping point selected to make yield and lactation persistency phenotypically uncorrelated it should be possible to use 305-d actual yield and lactation persistency to predict yield for lactations with later endpoints. The objectives of this study were to calculate (co)variance components and breeding values for best predictions of lactation persistency of milk (PM), fat (PF), protein (PP), and somatic cell score (PSCS) in breeds other than Holstein, and to demonstrate the calculation of prediction equations for 400-d actual milk yield. Data included lactations from Ayrshire, Brown Swiss, Guernsey (GU), Jersey (JE), and Milking Shorthorn (MS) cows calving since 1997. The number of sires evaluated ranged from 86 (MS) to 3,192 (JE), and mean sire estimated breeding value for PM ranged from 0.001 (Ayrshire) to 0.10 (Brown Swiss); mean estimated breeding value for PSCS ranged from -0.01 (MS) to -0.043 (JE). Heritabilities were generally highest for PM (0.09 to 0.15) and lowest for PSCS (0.03 to 0.06), with PF and PP having intermediate values (0.07 to 0.13). Repeatabilities varied considerably between breeds, ranging from 0.08 (PSCS in GU, JE, and MS) to 0.28 (PM in GU). Genetic correlations of PM, PF, and PP with PSCS were moderate and favorable (negative), indicating that increasing lactation persistency of yield traits is associated with decreases in lactation persistency of SCS, as expected. Genetic correlations among yield and lactation persistency were low to moderate and ranged from -0.55 (PP in GU) to 0.40 (PP in MS). Prediction equations for 400-d milk yield were calculated for each breed by regression of both 305-d yield and 305-d yield and lactation persistency on 400-d yield. Goodness-of-fit was very good for both models, but the addition of lactation persistency to the model significantly improved fit in all cases. Routine genetic evaluations for lactation persistency, as well as the development of prediction equations for several lactation end-points, may provide producers with tools to better manage their herds.

High-density genome-wide association study for residual feed intake in Holstein dairy cattle

Improving feed efficiency (FE) of dairy cattle may boost farm profitability and reduce the environmental footprint of the dairy industry. Residual feed intake (RFI), a candidate FE trait in dairy cattle, can be defined to be genetically uncorrelated with major energy sink traits (e.g., milk production, body weight) by including genomic predicted transmitting ability of such traits in genetic analyses for RFI. We examined the genetic basis of RFI through genome-wide association (GWA) analyses and post-GWA enrichment analyses and identified candidate genes and biological pathways associated with RFI in dairy cattle. Data were collected from 4,823 lactations of 3,947 Holstein cows in 9 research herds in the United States. Of these cows, 3,555 were genotyped and were imputed to a high-density list of 312,614 SNP. We used a single-step GWA method to combine information from genotyped and nongenotyped animals with phenotypes as well as their ancestors' information. The estimated genomic breeding values from a single-step genomic BLUP were back-solved to obtain the individual SNP effects for RFI. The proportion of genetic variance explained by each 5-SNP sliding window was also calculated for RFI. Our GWA analyses suggested that RFI is a highly polygenic trait regulated by many genes with small effects. The closest genes to the top SNP and sliding windows were associated with dry matter intake (DMI), RFI, energy homeostasis and energy balance regulation, digestion and metabolism of carbohydrates and proteins, immune regulation, leptin signaling, mitochondrial ATP activities, rumen development, skeletal muscle development, and spermatogenesis. The region of 40.7 to 41.5 Mb on BTA25 (UMD3.1 reference genome) was the top associated region for RFI. The closest genes to this region, CARD11 and EIF3B, were previously shown to be related to RFI of dairy cattle and FE of broilers, respectively. Another candidate region, 57.7 to 58.2 Mb on BTA18, which is associated with DMI and leptin signaling, was also associated with RFI in this study. Post-GWA enrichment analyses used a sum-based marker-set test based on 4 public annotation databases: Gene Ontology, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, Reactome pathways, and medical subject heading (MeSH) terms. Results of these analyses were consistent with those from the top GWA signals. Across the 4 databases, GWA signals for RFI were highly enriched in the biosynthesis and metabolism of amino acids and proteins, digestion and metabolism of carbohydrates, skeletal development, mitochondrial electron transport, immunity, rumen bacteria activities, and sperm motility. Our findings offer novel insight into the genetic basis of RFI and identify candidate regions and biological pathways associated with RFI in dairy cattle.

Single nucleotide polymorphisms associated with thermoregulation in lactating dairy cows exposed to heat stress

Dairy cows with increased rectal temperature experience lower milk yield and fertility. Rectal temperature during heat stress is heritable, so genetic selection for body temperature regulation could reduce effects of heat stress on production. One aim of the study was to validate the relationship between genotype and heat tolerance for single nucleotide polymorphisms (SNPs) previously associated with resistance to heat stress. A second aim was to identify new SNPs associated with heat stress resistance. Thermotolerance was assessed in lactating Holsteins during the summer by measuring rectal temperature (a direct measurement of body temperature regulation; n = 435), respiration rate (an indirect measurement of body temperature regulation, n = 450) and sweating rate (the major evaporative cooling mechanism in cattle, n = 455). The association between genotype and thermotolerance was evaluated for 19 SNPs previously associated with rectal temperature from a genomewide analysis study (GWAS), four SNPs previously associated with change in milk yield during heat stress from GWAS, 2 candidate gene SNPs previously associated with rectal temperature and respiration rate during heat stress (ATPA1A and HSP70A) and 66 SNPs in genes previously shown to be associated with reproduction, production or health traits in Holsteins. For SNPs previously associated with heat tolerance, regions of BTA4, BTA6 and BTA24 were associated with rectal temperature; regions of BTA6 and BTA24 were associated with respiration rate; and regions of BTA5, BTA26 and BTA29 were associated with sweating rate. New SNPs were identified for rectal temperature (n = 12), respiration rate (n = 8) and sweating rate (n = 3) from among those previously associated with production, reproduction or health traits. The SNP that explained the most variation were PGR and ASL for rectal temperature, ACAT2 and HSD17B7 for respiration rate, and ARL6IP1 and SERPINE2 for sweating rate. ARL6IP1 was associated with all three thermotolerance traits. In conclusion, specific genetic markers responsible for genetic variation in thermoregulation during heat stress in Holsteins were identified. These markers may prove useful in genetic selection for heat tolerance in Holstein cattle.

The use of high-frequency oscillatory ventilation (HFOV) and extracorporeal membrane oxygenation (ECMO) in the management of the term/near term infant with respiratory failure

The use of high-frequency oscillatory ventilation (HFOV) was evaluated as a rescue intervention in 122 consecutive infants meeting criteria for extracorporeal membrane oxygenation (ECMO). Fifty-three percent responded to HFOV and never required ECMO, 3/65 died. Infants who ultimately required ECMO had lower aortic and pulmonary peak flow velocities, and lower pulmonary acceleration, circumferential fiber shortening and shortening fraction than those who were successfully managed with HFOV. The use of HFOV with an appropriate strategy decreased the need for ECMO in this patient population.

Complications of extracorporeal membrane oxygenation in neonates

In cases of severe respiratory failure, cardiopulmonary bypass has been used as support until cardiac and pulmonary recovery occurs. We report the Wilford Hall USAF Medical Center experience with extracorporeal membrane oxygenation (ECMO) and its associated complications. From July 1985 to March 1989, 57 neonates were placed on membrane oxygenators. The overall survival was 79%. Technical complications encountered included catheter-related problems, mechanical complications of the pump apparatus, and hemorrhage. Hemorrhagic complications were the most frequent and devastating complications encountered. Intracranial hemorrhage accounted for six deaths associated with bypass. Extracorporeal membrane oxygenation is successful in significantly improving survival of neonates whose predicted mortality approaches 100% with conventional treatment. The rate of infant mortality using the membrane oxygenator is not affected by technical complications related to catheter position, mechanical problems with the circuit, or hemorrhage, excluding intracranial hemorrhage. The major cause of death of infants receiving extracorporeal membrane oxygenation is the underlying disease process leading to cardiopulmonary failure.

Reducing animal sequencing redundancy by preferentially selecting animals with low-frequency haplotypes

Many studies leverage targeted whole-genome sequencing (WGS) experiments to identify rare and causal variants within populations. As a natural consequence of their experimental design, many of these surveys tend to sequence redundant haplotype segments due to their high frequency in the base population, and the variants discovered within sequencing data are difficult to phase. We propose a new algorithm, called inverse weight selection (IWS), that preferentially selects individuals based on the cumulative presence of rare frequency haplotypes to maximize the efficiency of WGS surveys. To test the efficacy of this method, we used genotype data from 112,113 registered Holstein bulls derived from the US national dairy database. We demonstrate that IWS is at least 6.8% more efficient than previously published methods in selecting the least number of individuals required to sequence all haplotype segments ≥4% frequency in the US Holstein population. We also suggest that future surveys focus on sequencing homozygous haplotype segments as a first pass to achieve a 50% reduction in cost with an added benefit of phasing variant calls efficiently. Together, this new selection algorithm and experimental design suggestion significantly reduce the overall cost of variant discovery through WGS experiments, making surveys for causal variants influencing disease and production even more efficient.

Urinary leukotriene E(4) excretion during the first month of life and subsequent bronchopulmonary dysplasia in premature infants

BACKGROUND

Inflammation plays an important role in the pathogenesis of bronchopulmonary dysplasia (BPD), but the exact nature of this inflammatory process is incompletely understood. Older infants with established BPD have higher levels of urinary leukotriene E(4) (LTE(4)) compared to healthy infants of the same age. This suggests that cysteinyl leukotrienes may play a role in the abnormalities seen in BPD.

OBJECTIVES

To measure urinary LTE(4) levels during the first month of life in premature infants, and to determine whether there are significant differences in premature infants who develop BPD, as compared to those who do not develop BPD.

DESIGN

Prospective, blinded, controlled study.

SETTING

Neonatal ICUs of a tertiary-care university hospital.

METHODS

Thirty-seven premature infants (< 33 weeks of gestational age) were enrolled prospectively at birth. Urinary LTE(4) levels were measured blinded, using a standard radioimmunoassay technique at 2 days, 7 days, and 28 days of life. At 1 month of age, infants were classified as with or without BPD, based on need for supplemental oxygen, and characteristic chest radiographs. Clinical features and urinary LTE(4) were compared between the two groups.

RESULTS

Mean +/- SD gestational age was 29 +/- 2.6 weeks. None of the infants had a family history of asthma. Thirteen of 37 infants were classified as having BPD at 28 days after birth. Mean gestational age in infants who developed BPD was 27 +/- 2.4 weeks, compared to 30 +/- 2 weeks in infants who did not develop BPD (p < 0.05). In infants with BPD, mean urinary LTE(4) levels of urinary creatinine were 1,762 +/- 2,003 pg/mg, 1,236 +/- 992 pg/mg, and 5,541 +/- 5,146 pg/mg at days 2, 7, and 28, respectively, compared to 1,304 +/- 1,195 pg/mg, 1,158 +/- 1,133 pg/mg, and 2,800 +/- 2,080 pg/mg in infants without BPD. LTE(4) levels at 2 days, 7 days, and 28 days did not correlate with the subsequent development of BPD. LTE(4) levels at day 28 were significantly higher than LTE(4) levels at day 2 and day 7 in both groups, even after correcting for gestational age or birth weight (p < 0.05). There was significant inverse correlation between LTE(4) levels at day 2 with gestational age and birth weight (p < 0.05). All 13 infants with BPD received steroid pulses, compared to 3 of 26 infants without BPD. Gestational age and use of postnatal steroid pulses, diuretics, and theophylline (for apnea of prematurity) were significantly associated with each other and with the subsequent development of BPD.

CONCLUSION

Urinary LTE(4) levels measured on the second day of life in very-low-birth-weight infants inversely correlate with gestational age and birth weight. Urinary LTE(4) levels may reflect lung injury and/or inflammation in premature infants, not necessarily related to BPD as it is presently defined.

Short communication: Projecting milk yield using best prediction and the MilkBot lactation model

The accuracy and precision of 3 lactation models was estimated by summarizing means and variability in projection error for next-test milk and actual 305-d milk yield (M305) for 50-d intervals in a large Dairy Herd Improvement Association data set. Lactations were grouped by breed (Holstein, Jersey, and crossbred) and parity (first vs. later). A smaller, single-herd data set with both Dairy Herd Improvement Association data and daily milk weights was used to compare M305 calculated from test-day data with M305 computed by summing daily milk weights. The lactation models tested were best prediction (BP), the nonlinear MilkBot (MB) model, and a null model (NM) based on a stepwise function. The accuracy of the models was ranked (best to worst) MB, BP, and NM for later-parity cows and MB, NM, and BP for first-parity cows, with MB achieving accuracy in projecting daily milk of 0.5 kg or better in most groups. The models generally showed better accuracy after 50 d in milk. Best prediction and NM had low accuracy for crossbred cows and first-parity Holstein and Jersey cows. The MB model appears to be more precise than BP, and NM had low precision, especially for M305. Regression of model-generated M305 on summed M305 showed BP and MB to be equally efficient in ranking lactations, but MB was better at quantifying differences.

Short communication: best prediction of 305-day lactation yields with regional and seasonal effects

In the United States, lactation yields are calculated using best prediction (BP), a method in which test-day (TD) data are compared with breed- and parity-specific herd lactation curves that do not account for differences among regions of the country or seasons of calving. Complete data from 538,090 lactations of 348,123 Holstein cows with lactation lengths between 250 and 500 d, records made in a single herd, at least 5 reported TD, and twice-daily milking were extracted from the national dairy database and used to construct regional and seasonal lactation curves. Herds were assigned to 1 of 7 regions of the country, individual lactations were assigned to 3-mo seasons of calving, and lactation curves for milk, fat, and protein yields were estimated by parity group for regions, seasons, and seasons within regions. Multiplicative pre-adjustment factors (MF) also were computed. The resulting lactation curves and MF were tested on a validation data set of 891,806 lactations from 400,000 Holstein cows sampled at random from the national dairy database. Mature-equivalent milk, fat, and protein yields were calculated using the standard and adjusted curves and MF, and differences between 305-d mature-equivalent yields were tested for significance. Yields calculated using 50-d intervals from 50 to 250 d in milk (DIM) and using all TD to 500 DIM allowed comparisons of predictions for records in progress (RIP). Differences in mature-equivalent milk ranged from 0 to 51 kg and were slightly larger for first-parity than for later parity cows. Milk and components yields did not differ significantly in any case. Correlations of yields for 50-d intervals with those using all TD were similar across analyses. Yields for RIP were slightly more accurate when adjusted for regional and seasonal differences.

Major quantitative trait loci influencing milk production and conformation traits in Guernsey dairy cattle detected on Bos taurus autosome 19

The goal of this study was to identify potential quantitative trait loci (QTL) for 27 production, fitness, and conformation traits of Guernsey cattle through genome-wide association (GWA) analyses, with extra emphasis on BTA19, where major QTL were observed for several traits. Animals' de-regressed predicted transmitting abilities (PTA) from the December 2018 traditional US evaluation were used as phenotypes. All of the Guernsey cattle included in the QTL analyses were predictor animals in the reference population, ranging from 1,077 to 1,685 animals for different traits. Single-trait GWA analyses were carried out by a mixed-model approach for all 27 traits using imputed high-density genotypes. A major QTL was detected on BTA19, influencing several milk production traits, conformation traits, and livability of Guernsey cattle, and the most significant SNP lie in the region of 26.2 to 28.3 Mb. The myosin heavy chain 10 (MYH10) gene residing within this region was found to be highly associated with milk production and body conformation traits of dairy cattle. After the initial GWA analyses, which suggested that many significant SNP are in linkage with one another, conditional analyses were used for fine mapping. The top significant SNP on BTA19 were fixed as covariables in the model, one at a time, until no more significant SNP were detected on BTA19. After this fine-mapping approach was applied, only 1 significant SNP was detected on BTA19 for most traits, but multiple, independent significant SNP were found for protein yield, dairy form, and stature. In addition, the haplotype that hosts the major QTL on BTA19 was traced to a US Guernsey born in 1954. The haplotype is common in the breed, indicating a long-term influence of this QTL on the US Guernsey population.

Genetic and nongenetic profiling of milk pregnancy-associated glycoproteins in Holstein cattle

Pregnancy-associated glycoproteins (PAG) are secreted by the trophoblast and are detectable in maternal circulation around the time of attachment of the fetal placenta, as well as in blood and milk throughout gestation. The understanding of the genetic mechanisms controlling PAG levels can confer advantages for livestock breeding programs given the precocity and the ease of obtaining this phenotype from routine pregnancy diagnosis. The aim of this study was to characterize PAG levels by estimating genetic parameters and correlations with other dairy traits, and to identify genomic regions and candidate genes associated with PAG levels in milk. The PAG data consisted of pregnancy diagnoses using commercial assays from 2012 to 2017, and genotype data consisted of 54,123 SNP markers for 2,352 individuals (embryos and dams). The model included contemporary group (herd, year, and season) and embryo age as fixed effects, and random embryonic (direct) and maternal (indirect) genetic effects. Using genomic data, the estimated heritability for direct and maternal genetic effects (± standard deviations) were 0.23 ± 0.05 and 0.11 ± 0.05, respectively. The genetic correlation between these effects was almost zero (0.001 ± 0.06). A preliminary analysis revealed low correlations between milk PAG levels and other dairy traits. The genome-wide association analysis was performed using 2 approaches: single-marker and single-step using all markers. Four genomic regions with direct genetic effects were detected on Bos taurus autosome (BTA) 6, BTA7, BTA19, and BTA29 of the embryonic genome. The BTA29 locus was within the bovine PAG gene cluster, suggesting a cis-regulatory quantitative trait locus on the PAG expression. However, other associations, without an obvious link to PAG expression, could be related to the transportation of PAG and their concentration in milk. Only 1 region from the maternal genome, on BTA4, had a significant indirect effect, where WNT2 is a candidate gene related to placenta vascularization and gestation health. Collectively, our results suggest a moderate genetic control of milk PAG levels from both maternal and fetal genomes, but larger studies are needed to fully evaluate the usefulness of milk PAG in the genetic evaluation of fetal growth and cow fertility.

Marker selection and genomic prediction of economically important traits using imputed high-density genotypes for 5 breeds of dairy cattle

Marker sets used in US dairy genomic predictions were previously expanded by including high-density (HD) or sequence markers with the largest effects for Holstein breed only. Other non-Holstein breeds lacked enough HD genotyped animals to be used as a reference population at that time, and thus were not included in the genomic prediction. Recently, numbers of non-Holstein breeds genotyped using HD panels reached an acceptable level for imputation and marker selection, allowing HD genomic prediction and HD marker selection for Holstein plus 4 other breeds. Genotypes for 351,461 Holsteins, 347,570 Jerseys, 42,346 Brown Swiss, 9,364 Ayrshires (including Red dairy cattle), and 4,599 Guernseys were imputed to the HD marker list that included 643,059 SNP. The separate HD reference populations included Illumina BovineHD (San Diego, CA) genotypes for 4,012 Holsteins, 407 Jerseys, 181 Brown Swiss, 527 Ayrshires, and 147 Guernseys. The 643,059 variants included the HD SNP and all 79,254 (80K) genetic markers and QTL used in routine national genomic evaluations. Before imputation, approximately 91 to 97% of genotypes were unknown for each breed; after imputation, 1.1% of Holstein, 3.2% of Jersey, 6.7% of Brown Swiss, 4.8% of Ayrshire, and 4.2% of Guernsey alleles remained unknown due to lower density haplotypes that had no matching HD haplotype. The higher remaining missing rates in non-Holstein breeds are mainly due to fewer HD genotyped animals in the imputation reference populations. Allele effects for up to 39 traits were estimated separately within each breed using phenotypic reference populations that included up to 6,157 Jersey males and 110,130 Jersey females. Correlations of HD with 80K genomic predictions for young animals averaged 0.986, 0.989, 0.985, 0.992, and 0.978 for Jersey, Ayrshire, Brown Swiss, Guernsey, and Holstein breeds, respectively. Correlations were highest for yield traits (about 0.991) and lowest for foot angle and rear legs-side view (0.981and 0.982, respectively). Some HD effects were more than twice as large as the largest 80K SNP effect, and HD markers had larger effects than nearby 80K markers for many breed-trait combinations. Previous studies selected and included markers with large effects for Holstein traits; the newly selected HD markers should also improve non-Holstein and crossbred genomic predictions and were added to official US genomic predictions in April 2020.

Inheritance of a mutation causing neuropathy with splayed forelimbs in Jersey cattle

A new undesirable genetic factor, neuropathy with splayed forelimbs (JNS), has been identified recently in the Jersey breed. Calves affected with JNS are unable to stand on splayed forelimbs that exhibit significant extensor rigidity and excessive lateral abduction at birth. Affected calves generally are alert at birth but exhibit neurologic symptoms, including spasticity of head and neck and convulsive behavior. Other symptoms reported include dislocated shoulders, congenital craniofacial anomalies, and degenerative myelopathy. Inheritance of an undesirable genetic factor was determined from a study of 16 affected calves reported by Jersey breeders across the United States. All of their pedigrees traced back on both paternal and maternal sides to a common ancestor born in 1995. Genotypes revealed that JNS is attributable to a specific haplotype on Bos taurus autosome 6. Currently 8.2% of the genotyped US Jersey population are carriers of the haplotype. Sequencing of the region of shared homozygosity revealed missense variant rs1116058914 at base 60,158,901 of the ARS-UCD1.2 reference map as the most concordant with the genetic condition and the most likely cause. The single-base G to A substitution is in the coding region of the last exon of UCHL1, which is conserved across species. Mutations in humans and gene knockouts in mice cause similar recessive symptoms and muscular degeneration. Since December 2020, carrier status has been tracked using the identified haplotype and reported for all 459,784 genotyped Jersey animals. With random mating, about 2,200 affected calves per year with losses of about $250,000 would result from the 1.3 million US Jersey cows in the national population. Selection and mating programs can reduce numbers of JNS-affected births using either the haplotype status or a direct gene test in the future. Breeders should report calf abnormalities to their breed association to help discover new defects such as JNS.

Short communication: Phenotypic and genetic effects of the polled haplotype on yield, longevity, and fertility in US Brown Swiss, Holstein, and Jersey cattle

Phenotypes from the December 2018 US national genetic evaluations were used to compute effects of the polled haplotype in US Brown Swiss (BS), Holstein (HO), and Jersey (JE) cattle on milk, fat, and protein yields, somatic cell score, single-trait productive life, daughter pregnancy rate, heifer conception rate, and cow conception rate. Lactation records pre-adjusted for nongenetic factors and direct genomic values were used to estimate phenotypic and genetic effects of the polled haplotype, respectively. No phenotypic or direct genomic values effects were different from zero for any trait in any breed. Genomic PTA (gPTA) for the lifetime net merit (NM$) selection index of bulls born since January 1, 2012, that received a marketing code from the National Association of Animal Breeders (Madison, WI), and cows born on or after January 1, 2015, were compared to determine whether there was a systematic benefit to polled or horned genetics. Horned bulls had the highest average gPTA for NM$ in all 3 breeds, but that difference was significant only in HO and JE (HO: 615.4 ± 1.9, JE: 402.3 ± 3.4). Homozygous polled BS cows had significantly higher average gPTA for NM$ than their heterozygous polled or horned contemporaries (PP = 261.4 ± 43.5, Pp = 166.1 ± 13.7, pp = 174.1 ± 1.8), but the sample size was very small (n = 9). In HO and JE, horned cows had higher gPTA for NM$ (HO = 378.3 ± 0.2, JE = 283.3 ± 0.3). Selection for polled cattle should not have a detrimental effect on yield, fertility, or longevity, but these differences show that, in the short term, selection for polled over horned cattle will result in lower rates of genetic gain.

Illicit drug use and self-reported vision loss among military service members or veterans

INTRODUCTION

Little is known about differences in vision loss prevalence among service members or veterans (SMVs) and civilians; further, no study has compared vision loss risk factors in these two populations. As such, we seek to fill this gap in the literature.

METHODS

In this cross sectional study, we obtained data on 106 SMVs and 1572 civilians from the 2013-2018 National Health and Nutrition Examination Surveys. We compared the prevalence of or mean values of vision loss risk factors between SMVs and civilians using the Wald χ² statistic or Kruskal-Wallis test. Further, we examined the relative strength of 17 vision loss risk factors in predicting self-reported vision loss via Firth's logistic regression.

RESULTS

SMVs had a significantly higher prevalence of illicit drug use (20.75% vs 13.62%) and HIV (1.89% vs 0.41%), while civilians had a higher prevalence of poor dietary habits (7.61% vs 13.21%). SMVs also had higher mean values of systolic blood pressure (125.85 vs 122.53 mmHg), pack years of cigarette smoking (8.29 vs 4.25), and sedentary minutes per day (379.15 vs 337.07 min). More SMVs (8.49%) self-reported vision loss than civilians (4.48%). After adjustment for covariates, illicit drug use (adjusted β coefficient=0.72, p=0.02) was associated with self-reported vision loss.

CONCLUSIONS

This study indicates that self-reported vision loss among SMVs is more prevalent than among civilians, and vision loss in SMVs is associated with severe or prolonged illicit drug use.

Development of polymorphic markers in the immune gene complex loci of cattle

The addition of cattle health and immunity traits to genomic selection indices holds promise to increase individual animal longevity and productivity, and decrease economic losses from disease. However, highly variable genomic loci that contain multiple immune-related genes were poorly assembled in the first iterations of the cattle reference genome assembly and underrepresented during the development of most commercial genotyping platforms. As a consequence, there is a paucity of genetic markers within these loci that may track haplotypes related to disease susceptibility. By using hierarchical assembly of bacterial artificial chromosome inserts spanning 3 of these immune-related gene regions, we were able to assemble multiple full-length haplotypes of the major histocompatibility complex, the leukocyte receptor complex, and the natural killer cell complex. Using these new assemblies and the recently released ARS-UCD1.2 reference, we aligned whole-genome shotgun reads from 125 sequenced Holstein bulls to discover candidate variants for genetic marker development. We selected 124 SNPs, using heuristic and statistical models to develop a custom genotyping panel. In a proof-of-principle study, we used this custom panel to genotype 1,797 Holstein cows exposed to bovine tuberculosis (bTB) that were the subject of a previous GWAS study using the Illumina BovineHD array. Although we did not identify any significant association of bTB phenotypes with these new genetic markers, 2 markers exhibited substantial effects on bTB phenotypic prediction. The models and parameters trained in this study serve as a guide for future marker discovery surveys particularly in previously unassembled regions of the cattle genome.

256 GENOMIC EVALUATION OF FERTILITY TRAITS AND DISCOVERY OF HAPLOTYPES THAT AFFECT FERTILITY OF US DAIRY CATTLE

Genomic evaluations of dairy cattle became official in the United States in January 2009 for Holsteins and Jerseys, and later for Brown Swiss, Ayrshires, and Guernseys. Up to 33 yield, fitness, calving, and conformation traits are evaluated, and the fertility traits included daughter pregnancy rate and heifer and cow conception rates. Additional fertility traits, such as age at first calving and days from calving to first insemination, also are being studied. Male fertility (sire conception rate) is evaluated phenotypically rather than through genomics. Over 1 million animals have genotypes in the national database, which reflects collaboration with Canada and Europe. Most of the genotypes are from females and are from genotyping chips with <30 000 single nucleotide polymorphisms (SNP). To combine data across chips, genotypes are imputed to a set of >77 000 SNP. The imputation process involves dividing the chromosome into segments of approximately equal length and determining the paternal or maternal origin of the alleles. Because some segments were never homozygous, they were assumed to contain an abnormality that resulted in early embryonic death. If a decrease in sire conception rate could be associated with a bull that was a carrier of such a chromosomal segment, the haplotype was designated as affecting fertility. Once the region was identified, bioinformatic analysis was used to discover the causative variant for many of those haplotypes. Accuracy of genomic evaluations is determined by size of the reference population and heritability of the trait. The reference population for Holsteins includes >180 000 bulls and cows. Because fertility traits have low heritabilities, genomic information is particularly useful in improving evaluation accuracy. Accuracy of fertility evaluations is expected to increase further by discovering causative variants for various aspects of conception and gestation through investigation of sequence data.

Open lung biopsy in the critically ill newborn

Experience with 17 open lung biopsies in critically ill premature neonates was reviewed. Despite their small size, prematurity, and near maximal ventilator requirements, the infants suffered no significant complications. In three cases, an infectious agent was identified. In one case end-stage lung fibrosis associated with persistent (greater than 3 months), severe respiratory failure prompted termination of support. In the remainder of the cases, definitively ruling out infection allowed the confident trial of a course of steroids in an attempt to treat bronchopulmonary dysplasia. Unlike older patients, the definitive diagnosis of no infection in the premature neonate is just as informative as the diagnosis of an infection. Used judiciously, open lung biopsy can be performed in the premature infant with acceptable morbidity and mortality.