Enlarging a training set for genomic selction by imputation of un-genotyped animals in populations of varying genetic architecture

Background: The most common application of imputation is to infer genotypes of a high-density panel of markers on animals that are genotyped for a low-density panel. However, the increase in accuracy of genomic predictions resulting from an increase in the number of markers tends to reach a plateau beyond a certain density. Another application of imputation is to increase the size of the training set with un-genotyped animals. This strategy can be particularly successful when a set of closely related individuals are genotyped. ----- Methods: Imputation on completely un-genotyped dams was performed using known genotypes from the sire of each dam, one offspring and the offspring’s sire. Two methods were applied based on either allele or haplotype frequencies to infer genotypes at ambiguous loci. Results of these methods and of two available software packages were compared. Quality of imputation under different population structures was assessed. The impact of using imputed dams to enlarge training sets on the accuracy of genomic predictions was evaluated for different populations, heritabilities and sizes of training sets. ----- Results: Imputation accuracy ranged from 0.52 to 0.93 depending on the population structure and the method used. The method that used allele frequencies performed better than the method based on haplotype frequencies. Accuracy of imputation was higher for populations with higher levels of linkage disequilibrium and with larger proportions of markers with more extreme allele frequencies. Inclusion of imputed dams in the training set increased the accuracy of genomic predictions. Gains in accuracy ranged from close to zero to 37.14%, depending on the simulated scenario. Generally, the larger the accuracy already obtained with the genotyped training set, the lower the increase in accuracy achieved by adding imputed dams. ----- Conclusions: Whenever a reference population resembling the family configuration considered here is available, imputation can be used to achieve an extra increase in accuracy of genomic predictions by enlarging the training set with completely un-genotyped dams. This strategy was shown to be particularly useful for populations with lower levels of linkage disequilibrium, for genomic selection on traits with low heritability, and for species or breeds for which the size of the reference population is limited.

Physiological analysis of central and peripheral insect circadian pacemaker neurons

Alle bisher untersuchten Lebewesen besitzen (circadiane) innere Uhren, die eine endogene Perioden-länge von ungefähr 24 Stunden generieren. Eine innere Uhr kann über Zeitgeber mit der Umwelt synchronisiert werden und ermöglicht dem Organismus, rhythmische Umweltveränderungen vorweg zu nehmen. Neben einem zentralen Schrittmacher, der Physiologie und Verhalten des Organismus steuert, gibt es in unterschiedlichen Organen auch periphere Uhren, die die zeitlichen Abläufe in der spezifischen Funktion dieser Organe steuern. In dieser Arbeit sollten zentrale und periphere Schrittmacherneurone von Insekten physiologisch untersucht und verglichen werden. Die Neurone der akzessorischen Medulla (AME) von Rhyparobia maderae dienten als Modellsystem für zentrale Schrittmacher, während olfaktorische Rezeptorneurone (ORNs) von Manduca sexta als Modellsystem für periphere Schrittmacher dienten. Die zentralen Schrittmacherneurone wurden in extrazellulären Ableitungen an der isolierten AME (Netzwerkebene) und in Patch-Clamp Experimenten an primären AME Zellkulturen (Einzelzellebene) untersucht. Auf Netzwerkebene zeigten sich zwei charakteristische Aktivitätsmuster: regelmäßige Aktivität und Wechsel zwischen hoher und niedriger Aktivität (Oszillationen). Es wurde gezeigt, dass Glutamat ein Neurotransmitter der weitverbreiteten inhibitorischen Synapsen der AME ist, und dass in geringem Maße auch exzitatorische Synapsen vorkommen. Das Neuropeptid pigment-dispersing factor (PDF), das von nur wenigen AME Neuronen exprimiert wird und ein wichtiger Kopplungsfaktor im circadianen System ist, führte zu Hemmungen, Aktivierungen oder Oszillationen. Die Effekte waren transient oder langanhaltend und wurden wahrscheinlich durch den sekundären Botenstoff cAMP vermittelt. Ein Zielmolekül von cAMP war vermutlich exchange protein directly activated by cAMP (EPAC). Auf Einzelzellebene wurde gezeigt, dass die meisten AME Neurone depolarisiert waren und deshalb nicht feuerten. Die Analyse von Strom-Spannungs-Kennlinien und pharmakologische Experimente ergaben, dass unterschiedliche Ionenkanäle vorhanden waren (Ca2+, Cl-, K+, Na+ Kanäle sowie nicht-spezifische Kationenkanäle). Starke, bei hohen Spannungen aktivierende Ca2+ Ströme (ICa) könnten eine wichtige Rolle bei Ca2+-abhängiger Neurotransmitter-Ausschüttung, Oszillationen, und Aktionspotentialen spielen. PDF hemmte unterschiedliche Ströme (ICa, IK und INa) und aktivierte nicht-spezifische Kationenströme (Ih). Es wurde angenommen, dass simultane PDF-abhängige Hyper- und Depolarisationen rhythmische Membranpotential-Oszillationen verursachen. Dieser Mechanismus könnte eine Rolle bei PDF-abhängigen Synchronisationen spielen. Die Analyse peripherer Schrittmacherneurone konzentrierte sich auf die Charakterisierung des olfaktorischen Corezeptors von M. sexta (MsexORCO). In anderen Insekten ist ORCO für die Membran-Insertion von olfaktorischen Rezeptoren (ORs) erforderlich. ORCO bildet Komplexe mit den ORs, die in heterologen Expressionssystemen als Ionenkanäle fungieren und Duft-Antworten vermitteln. Es wurde die Hypothese aufgestellt, dass MsexORCO in pheromonsensitiven ORNs in vivo nicht als Teil eines ionotropen Rezeptors sondern als Schrittmacherkanal fungiert, der unterschwellige Membranpotential-Oszillationen generiert. MsexORCO wurde mit vermeintlichen Pheromonrezeptoren in human embryonic kidney (HEK 293) Zellen coexprimiert. Immuncytochemie und Ca2+ Imaging Experimente zeigten sehr schwache Expressionsraten. Trotzdem war es möglich zu zeigen, dass MsexORCO wahrscheinlich ein spontan-aktiver, Ca2+-permeabler Ionenkanal ist, der durch den ORCO-Agonisten VUAA1 und cyclische Nucleotide aktiviert wird. Außerdem wiesen die Experimente darauf hin, dass MsexOR-1 offensichtlich der Bombykal-Rezeptor ist. Eine weitere Charakterisierung von MsexORCO in primären M. sexta ORN Zellkulturen konnte nicht vollendet werden, weil die ORNs nicht signifikant auf ORCO-Agonisten oder -Antagonisten reagierten.

Heat stress response for physiological traits in dairy and dual purpose cattle populations on phenotypic and genetic scales

The main objective of this thesis was to determine the potential impact of heat stress (HS) on physiological traits of lactating cows and semen quality of bulls kept in a temperate climate. The thesis is comprised of three studies. An innovative statistical modeling aspect common to all three studies was the application of random regression methodology (RRM) to study the phenotypic and genetic trajectory of traits in dependency of a continuous temperature humidity index (THI). In the first study, semen quality and quantity traits of 562 Holstein sires kept on an AI station in northwestern Germany were analyzed in the course of THI calculated from data obtained from the nearest weather station. Heat stress was identified based on a decline in semen quality and quantity parameters. The identified general HS threshold (THI = 60) and the thermoneutal zone (THI in the range from 50 to 60) for semen production were lower than detected in studies conducted in tropical and subtropical climates. Even though adult bulls were characterized by higher semen productivity compared to younger bulls, they responded with a stronger semen production loss during harsh environments. Heritabilities (low to moderate range) and additive genetic variances of semen characteristics varied with different levels of THI. Also, based on genetic correlations genotype, by environment interactions were detected. Taken together, these findings suggest the application of specific selection strategies for specific climate conditions. In the second study, the effect of the continuous environmental descriptor THI as measured inside the barns on rectal temperatures (RT), skin temperatures (ST), vaginal temperatures (VT), respiration rates (RR), and pulse rate (PR) of lactating Holstein Friesian (HF) and dual-purpose German black pied cattle (DSN) was analyzed. Increasing HS from THI 65 (threshold) to THI 86 (maximal THI) resulted in an increase of RT by 0.6 °C (DSN) and 1 °C (HF), ST by 3.5 °C (HF) and 8 °C (DSN), VT by 0.3 °C (DSN), and RR by 47 breaths / minute (DSN), and decreased PR by 7 beats / minute (DSN). The undesired effects of rising THI on physiological traits were most pronounced for cows with high levels of milk yield and milk constituents, cows in early days in milk and later parities, and during summer seasons in the year 2014. In the third study of this dissertation, the genetic components of the cow’s physiological responses to HS were investigated. Heat stress was deduced from indoor THI measurements, and physiological traits were recorded on native DSN cows and their genetically upgraded crosses with Holstein Friesian sires in two experimental herds from pasture-based production systems reflecting a harsh environment of the northern part of Germany. Although heritabilities were in a low range (from 0.018 to 0.072), alterations of heritabilities, repeatabilities, and genetic components in the course of THI justify the implementation of genetic evaluations including heat stress components. However, low repeatabilities indicate the necessity of using repeated records for measuring physiological traits in German cattle. Moderate EBV correlations between different trait combinations indicate the potential of selection for one trait to simultaneously improve the other physiological attributes. In conclusion, bulls of AI centers and lactating cows suffer from HS during more extreme weather conditions also in the temperate climate of Northern Germany. Monitoring physiological traits during warm and humid conditions could provide precious information for detection of appropriate times for implementation of cooling systems and changes in feeding and management strategies. Subsequently, the inclusion of these physiological traits with THI specific breeding values into overall breeding goals could contribute to improving cattle adaptability by selecting the optimal animal for extreme hot and humid conditions. Furthermore, the recording of meteorological data in close distance to the cow and visualizing the surface body temperature by infrared thermography techniques might be helpful for recognizing heat tolerance and adaptability in cattle.