Whole genome sequencing reveals a complex introgression history and the basis of adaptation to subarctic climate in wild sheep.

To predict species responses to anthropogenic disturbances and climate change, it is reasonable to use species with high sensitivity to such factors. Snow sheep (Ovis nivicola) could represent a good candidate for this; as the only large herbivore species adapted to the cold and alpine habitats of northeastern Siberia, it plays a crucial role in its ecosystem. Despite having an extensive geographical distribution among all ovine species, it is one of the least studied. In this study, we sequenced and analysed six genomes of snow sheep in combination with all other wild sheep species to infer key aspects of their evolutionary history and unveil the genetic basis of their adaptation to subarctic environments. Despite their large census population size, snow sheep genomes showed remarkably low heterozygosity, which could reflect the effect of isolation and historical bottlenecks that we inferred using the pairwise sequential Markovian coalescent and runs of homozygosity. F4 -statistics indicated instances of introgression involving snow sheep with argali (Ovis ammon) and Dall (Ovis dalli) sheep, suggesting that these species might have been more widespread during the Pleistocene. Furthermore, the introgressed segments, which were identified using mainly minimum relative node depth, covered genes associated with immunity, adipogenesis and morphology-related traits, representing potential targets of adaptive introgression. Genes related to mitochondrial functions and thermogenesis associated with adipose tissue were identified to be under selection. Overall, our data suggest introgression as a mechanism facilitating adaptation in wild sheep species and provide insights into the genetic mechanisms underlying cold adaptation in snow sheep.

Genomic diversity and population structure of the indigenous Greek and Cypriot cattle populations.

BACKGROUND: The indigenous cattle populations from Greece and Cyprus have decreased to small numbers and are currently at risk of extinction due to socio-economic reasons, geographic isolation and crossbreeding with commercial breeds. This study represents the first comprehensive genome-wide analysis of 10 indigenous cattle populations from continental Greece and the Greek islands, and one from Cyprus, and compares them with 104 international breeds using more than 46,000 single nucleotide polymorphisms (SNPs). RESULTS: We estimated several parameters of genetic diversity (e.g. heterozygosity and allelic diversity) that indicated a severe loss of genetic diversity for the island populations compared to the mainland populations, which is mainly due to the declining size of their population in recent years and subsequent inbreeding. This high inbreeding status also resulted in higher genetic differentiation within the Greek and Cyprus cattle group compared to the remaining geographical breed groups. Supervised and unsupervised cluster analyses revealed that the phylogenetic patterns in the indigenous Greek breeds were consistent with their geographical origin and historical information regarding crosses with breeds of Anatolian or Balkan origin. Cyprus cattle showed a relatively high indicine ancestry. Greek island populations are placed close to the root of the tree as defined by Gir and the outgroup Yak, whereas the mainland breeds share a common historical origin with Busa. Unsupervised clustering and D-statistics analyses provided strong support for Bos indicus introgression in almost all the investigated local cattle breeds along the route from Anatolia up to the southern foothills of the Alps, as well as in most cattle breeds along the Apennine peninsula to the southern foothills of the Alps. CONCLUSIONS: All investigated Cyprus and Greek breeds present complex mosaic genomes as a result of historical and recent admixture events between neighbor and well-separated breeds. While the contribution of some mainland breeds to the genetic diversity pool seems important, some island and fragmented mainland breeds suffer from a severe decline of population size and loss of alleles due to genetic drift. Conservation programs that are a compromise between what is feasible and what is desirable should focus not only on the still highly diverse mainland breeds but also promote and explore the conservation possibilities for island breeds.

Conservation of a domestic metapopulation structured into related and partly admixed strains.

Preservation of genetic diversity is one of the most pressing challenges in the planetary boundaries concept. Within this context, we focused on genetic diversity in a native, unselected and highly admixed domesticated metapopulation. A set of 1,828 individuals from 60 different cattle breeds was analysed using a medium density SNP chip. Among these breeds, 14 Busa strains formed a metapopulation represented by 350 individuals, while the remaining 46 breeds represented the global cattle population. Genetic analyses showed that the scarcely selected and less differentiated Busa metapopulation contributed a substantial proportion (52.6%) of the neutral allelic diversity to this global taurine population. Consequently, there is an urgent need for synchronized maintenance of this highly fragmented domestic metapopulation, which is distributed over several countries without sophisticated infrastructure and highly endangered by continuous replacement crossing as part of the global genetic homogenization process. This study collected and evaluated samples, data and genomewide information and developed genome-assisted cross-border conservation concepts. To detect and maintain genetic integrity of the metapopulation strains, we designed and applied a composite test that combines six metrics based on additive genetic relationships, a nearest neighbour graph and the distribution of semiprivate alleles. Each metric provides distinct information components about past admixture events and offers an objective and powerful tool for the detection of admixed outliers. The here developed conservation methods and presented experiences could easily be adapted to comparable conservation programmes of domesticated or other metapopulations bred and kept in captivity or under some other sort of human control.

Remapping of the belted phenotype in cattle on BTA3 identifies a multiplication event as the candidate causal mutation.

BACKGROUND: It has been known for almost a century that the belted phenotype in cattle follows a pattern of dominant inheritance. In 2009, the approximate position of the belt locus in Brown Swiss cattle was mapped to a 922-kb interval on bovine chromosome 3 and, subsequently, assigned to a 336-kb haplotype block based on an animal set that included, Brown Swiss, Dutch Belted (Lakenvelder) and Belted Galloway individuals. A possible candidate gene in this region i.e. HES6 was investigated but the causal mutation remains unknown. Thus, to elucidate the causal mutation of this prominent coat color phenotype, we decided to remap the belted phenotype in an independent animal set of several European bovine breeds, i.e. Gurtenvieh (belted Brown Swiss), Dutch Belted and Belted Galloway and to systematically scan the candidate region. We also checked the presence of the detected causal mutation in the genome of belted individuals from a Siberian cattle breed. RESULTS: A combined linkage disequilibrium and linkage analysis based on 110 belted and non-belted animals identified a candidate interval of 2.5 Mb. Manual inspection of the haplotypes in this region identified four candidate haplotypes that consisted of five to eight consecutive SNPs. One of these haplotypes overlapped with the initial 922-kb interval, whereas two were positioned proximal and one was positioned distal to this region. Next-generation sequencing of one heterozygous and two homozygous belted animals identified only one private belted candidate allele, i.e. a multiplication event that is located between 118,608,000 and 118,614,000 bp. Targeted locus amplification and quantitative real-time PCR confirmed an increase in copy number of this region in the genomes of both European (Belted Galloway, Dutch Belted and Gurtenvieh) and Siberian (Yakutian cattle) breeds. Finally, using nanopore sequencing, the exact breakpoints were determined at 118,608,362 and 118,614,132 bp. The closest gene to the candidate causal mutation (16 kb distal) is TWIST2. CONCLUSIONS: Based on our findings and those of a previously published study that identified the same multiplication event, a quadruplication on bovine chromosome 3 between positions 118,608,362 and 118,614,132 bp is the most likely candidate causal mutation for the belted phenotype in cattle.

Population structure and genetic diversity of 25 Russian sheep breeds based on whole-genome genotyping.

BACKGROUND: Russia has a diverse variety of native and locally developed sheep breeds with coarse, fine, and semi-fine wool, which inhabit different climate zones and landscapes that range from hot deserts to harsh northern areas. To date, no genome-wide information has been used to investigate the history and genetic characteristics of the extant local Russian sheep populations. To infer the population structure and genome-wide diversity of Russian sheep, 25 local breeds were genotyped with the OvineSNP50 BeadChip. Furthermore, to evaluate admixture contributions from foreign breeds in Russian sheep, a set of 58 worldwide breeds from publicly available genotypes was added to our data. RESULTS: We recorded similar observed heterozygosity (0.354-0.395) and allelic richness (1.890-1.955) levels across the analyzed breeds and they are comparable with those observed in the worldwide breeds. Recent effective population sizes estimated from linkage disequilibrium five generations ago ranged from 65 to 543. Multi-dimensional scaling, admixture, and neighbor-net analyses consistently identified a two-step subdivision of the Russian local sheep breeds. A first split clustered the Russian sheep populations according to their wool type (fine wool, semi-fine wool and coarse wool). The Dagestan Mountain and Baikal fine-fleeced breeds differ from the other Merino-derived local breeds. The semi-fine wool cluster combined a breed of Romanian origin, Tsigai, with its derivative Altai Mountain, the two Romney-introgressed breeds Kuibyshev and North Caucasian, and the Lincoln-introgressed Russian longhaired breed. The coarse-wool group comprised the Nordic short-tailed Romanov, the long-fat-tailed outlier Kuchugur and two clusters of fat-tailed sheep: the Caucasian Mountain breeds and the Buubei, Karakul, Edilbai, Kalmyk and Tuva breeds. The Russian fat-tailed breeds shared co-ancestry with sheep from China and Southwestern Asia (Iran). CONCLUSIONS: In this study, we derived the genetic characteristics of the major Russian local sheep breeds, which are moderately diverse and have a strong population structure. Pooling our data with a worldwide genotyping set gave deeper insight into the history and origin of the Russian sheep populations.

Detection of two non-synonymous SNPs in SLC45A2 on BTA20 as candidate causal mutations for oculocutaneous albinism in Braunvieh cattle.

BACKGROUND: Cases of albinism have been reported in several species including cattle. So far, research has identified many genes that are involved in this eye-catching phenotype. Thus, when two paternal Braunvieh half-sibs with oculocutaneous albinism were detected on a private farm, we were interested in knowing whether their phenotype was caused by an already known gene/mutation. RESULTS: Analysis of genotyping data (50K) of the two albino individuals, their mothers and five other relatives identified a 47.61-Mb candidate haplotype on Bos taurus chromosome BTA20. Subsequent comparisons of the sequence of this haplotype with sequence data from four Braunvieh sires and the Aurochs genome identified two possible candidate causal mutations at positions 39,829,806 bp (G/A; R45Q) and 39,864,148 bp (C/T; T444I) that were absent in 1682 animals from various bovine breeds included in the 1000 bull genomes project. Both polymorphisms represent coding variants in the SLC45A2 gene, for which the human equivalent harbors numerous variants associated with oculocutaneous albinism type 4. We demonstrate an association of R45Q and T444I with the albino phenotype by targeted genotyping. CONCLUSIONS: Although the candidate gene SLC45A2 is known to be involved in albinism in different species, to date in cattle only mutations in the TYR and MITF genes were reported to be associated with albinism or albinism-like phenotypes. Thus, our study extends the list of genes that are associated with bovine albinism. However, further research and more samples from related animals are needed to elucidate if only one of these two single nucleotide polymorphisms or the combination of both is the actual causal variant.

Confirmation of a non-synonymous SNP in PNPLA8 as a candidate causal mutation for Weaver syndrome in Brown Swiss cattle.

BACKGROUND: Bovine progressive degenerative myeloencephalopathy (Weaver syndrome) is a neurodegenerative disorder in Brown Swiss cattle that is characterized by progressive hind leg weakness and ataxia, while sensorium and spinal reflexes remain unaffected. Although the causal mutation has not been identified yet, an indirect genetic test based on six microsatellite markers and consequent exclusion of Weaver carriers from breeding have led to the complete absence of new cases for over two decades. Evaluation of disease status by imputation of 41 diagnostic single nucleotide polymorphisms (SNPs) and a common haplotype published in 2013 identified several suspected carriers in the current breeding population, which suggests a higher frequency of the Weaver allele than anticipated. In order to prevent the reemergence of the disease, this study aimed at mapping the gene that underlies Weaver syndrome and thus at providing the basis for direct genetic testing and monitoring of today's Braunvieh/Brown Swiss herds. RESULTS: Combined linkage/linkage disequilibrium mapping on Bos taurus chromosome (BTA) 4 based on Illumina Bovine SNP50 genotypes of 43 Weaver-affected, 31 Weaver carrier and 86 Weaver-free animals resulted in a maximum likelihood ratio test statistic value at position 49,812,384 bp. The confidence interval (0.853 Mb) determined by the 2-LOD drop-off method was contained within a 1.72-Mb segment of extended homozygosity. Exploitation of whole-genome sequence data from two official Weaver carriers and 1145 other bulls that were sequenced in Run4 of the 1000 bull genomes project showed that only a non-synonymous SNP (rs800397662) within the PNPLA8 gene at position 49,878,773 bp was concordant with the Weaver carrier status. Targeted SNP genotyping confirmed this SNP as a candidate causal mutation for Weaver syndrome. Genotyping for the candidate causal mutation in a random sample of 2334 current Braunvieh animals suggested a frequency of the Weaver allele of 0.26 %. CONCLUSIONS: Through combined use of exhaustive sequencing data and SNP genotyping results, we were able to provide evidence that supports the non-synonymous mutation at position 49,878,773 bp as the most likely causal mutation for Weaver syndrome. Further studies are needed to uncover the exact mechanisms that underlie this syndrome.

Phylogenetic Analysis of Russian Native Sheep Breeds Based on mtDNA Sequences.

Eurasia is represented by all climatic zones and various environments. A unique breed variety of farm animals has been developed in Russia, whose territory covers a large area of the continent. A total of 69 local breeds and types of dairy, wool, and meat sheep (Ovis aries) are maintained here. However, the genetic diversity and maternal origin of these local breeds have not been comprehensively investigated. In this study, we describe the diversity and phylogeny of Russian sheep breeds inhabiting different geographical regions based on the analysis of complete sequences of mitochondrial genomes (mtDNA). Complete mtDNA sequences of the studied sheep were obtained using next-generation sequencing technology (NGS). All investigated geographical groups of sheep were characterized by high haplotype (Hd = 0.9992) and nucleotide diversity (π = 0.00378). Analysis of the AMOVA results showed that genetic diversity was majorly determined by within-population differences (77.87%). We identified 128 haplotypes in all studied sheep. Haplotypes belonged to the following haplogroups: B (64.8%), A (28.9%), C (5.5%), and D (0.8%). Haplogroup B was predominant in the western part of Russia. A high level of mtDNA polymorphism in the studied groups of local sheep indicates the presence of a significant reserve of unique genotypes in Russia, which is to be explored.

[Persistent, excessive crying in 5-month-old infants and the pre-, peri- and postnatal adversities of their mothers in a high-risk sample]. / Das exzessive Schreien bei 5 Monate alten Säuglingen und prä-, peri- und postnatale Belastungen ihrer Mütter in einer Hoch-Risiko-Stichprobe.

OBJECTIVE: To investigate the relationship between persistent, excessive crying in 5-month-old infants and the pre- and perinatal adversities as well as postpartal mood of their mothers. METHOD: A sample of 300 mother-child dyads was examined at infants' age of 18.5 weeks. All mothers exhibited psychosocial risks such as poverty, lack of social support, being underage, drug abuse or mental disorders. Excessive crying was assessed by the Wessel's «rule of threes¼. Pre-, peri- and postnatal problems were measured by self-report questionnaires. RESULTS: Multivariate data analysis revealed an increased risk for social adversities during pregnancy (OR = 17.66) and unwanted pregnancy (OR = 13.77). For the postnatal period persistent crying was associated with a higher rate of maternal postpartum depressive symptoms, maternal stress, dysfunctional mother-child interactions, perception of the infant as being «difficult¼ as well as bonding problems. CONCLUSIONS: The results point to the influence of prenatal stress in mothers and a primarily unwanted pregnancy on infants' persistent crying. An increased need for support is suggested in cases of considerably reduced well-being of the mother or bonding problems. The combined effect of prenatal variables and relationship variables influences the occurrence and perpetuation of early regulation problems.

Large-scale mitogenome sequencing reveals consecutive expansions of domestic taurine cattle and supports sporadic aurochs introgression.

The contribution of domestic cattle in human societies is enormous, making cattle, along with other essential benefits, the economically most important domestic animal in the world today. To expand existing knowledge on cattle domestication and mitogenome diversity, we performed a comprehensive complete mitogenome analysis of the species (802 sequences, 114 breeds). A large sample was collected in South-east Europe, an important agricultural gateway to Europe during Neolithization and a region rich in cattle biodiversity. We found 1725 polymorphic sites (810 singletons, 853 parsimony-informative sites and 57 indels), 701 unique haplotypes, a haplotype diversity of 0.9995 and a nucleotide diversity of 0.0015. In addition to the dominant T3 and several rare haplogroups (Q, T5, T4, T2 and T1), we have identified maternal line in Austrian Murbodner cattle that possess surviving aurochs' mitochondria haplotype P1 that diverged prior to the Neolithization process. This is convincing evidence for rare female-mediated adaptive introgression of wild aurochs into domesticated cattle in Europe. We revalidated the existing haplogroup classification and provided Bayesian phylogenetic inference with a more precise estimated divergence time than previously available. Occasionally, classification based on partial mitogenomes was not reliable; for example, some individuals with haplogroups P and T5 were not recognized based on D-loop information. Bayesian skyline plot estimates (median) show that the earliest population growth began before domestication in cattle with haplogroup T2, followed by Q (~10.0-9.5 kyBP), whereas cattle with T3 (~7.5 kyBP) and T1 (~3.0-2.5 kyBP) expanded later. Overall, our results support the existence of interactions between aurochs and cattle during domestication and dispersal of cattle in the past, contribute to the conservation of maternal cattle diversity and enable functional analyses of the surviving aurochs P1 mitogenome.

Mitochondrial DNA Analysis Clarifies Taxonomic Status of the Northernmost Snow Sheep (Ovis nivicola) Population.

Currently, the intraspecific taxonomy of snow sheep (Ovis nivicola) is controversial and needs to be specified using DNA molecular genetic markers. In our previous work using whole-genome single nucleotide polymorphism (SNP) analysis, we found that the population inhabiting Kharaulakh Ridge was genetically different from the other populations of Yakut subspecies to which it was usually referred. Here, our study was aimed at the clarification of taxonomic status of Kharaulakh snow sheep using mitochondrial cytochrome b gene. A total of 87 specimens from five different geographic locations of Yakut snow sheep as well as 20 specimens of other recognized subspecies were included in this study. We identified 19 haplotypes, two of which belonged to the population from Kharaulakh Ridge. Median-joining network and Bayesian tree analyses revealed that Kharaulakh population clustered separately from all the other Yakut snow sheep. The divergence time between Kharaulakh population and Yakut snow sheep was estimated as 0.48 ± 0.19 MYA. Thus, the study of the mtDNA cytb sequences confirmed the results of genome-wide SNP analysis. Taking into account the high degree of divergence of Kharaulakh snow sheep from other groups, identified by both nuclear and mitochondrial DNA markers, we propose to classify the Kharaulakh population as a separate subspecies.

Complete mitochondrial genomes of Karchaev goat (Capra hircus).

Karachaev goat (Capra hircus) is a local breed from North-Caucasus region, Russia. Here we present complete mitochondrial genome of Karachaev goat from the republic of Karachaevo-Cherkessia, Russia. The length of the studied sequence was 16,624 bp in size. It was shown that the studied specimen belonged to haplogroup A.

A de novo frameshift mutation in ZEB2 causes polledness, abnormal skull shape, small body stature and subfertility in Fleckvieh cattle.

Polledness in cattle is an autosomal dominant trait. Previous studies have revealed allelic heterogeneity at the polled locus and four different variants were identified, all in intergenic regions. In this study, we report a case of polled bull (FV-Polled1) born to horned parents, indicating a de novo origin of this polled condition. Using 50K genotyping and whole genome sequencing data, we identified on chromosome 2 an 11-bp deletion (AC_000159.1:g.52364063_52364073del; Del11) in the second exon of ZEB2 gene as the causal mutation for this de novo polled condition. We predicted that the deletion would shorten the protein product of ZEB2 by almost 91%. Moreover, we showed that all animals carrying Del11 mutation displayed symptoms similar to Mowat-Wilson syndrome (MWS) in humans, which is also associated with genetic variations in ZEB2. The symptoms in cattle include delayed maturity, small body stature and abnormal shape of skull. This is the first report of a de novo dominant mutation affecting only ZEB2 and associated with a genetic absence of horns. Therefore our results demonstrate undoubtedly that ZEB2 plays an important role in the process of horn ontogenesis as well as in the regulation of overall development and growth of animals.

The First Draft Genome Assembly of Snow Sheep (Ovis nivicola).

The snow sheep, Ovis nivicola, which is endemic to the mountain ranges of northeastern Siberia, are well adapted to the harsh cold climatic conditions of their habitat. In this study, using long reads of Nanopore sequencing technology, whole-genome sequencing, assembly, and gene annotation of a snow sheep were carried out. Additionally, RNA-seq reads from several tissues were also generated to supplement the gene prediction in snow sheep genome. The assembled genome was ∼2.62 Gb in length and was represented by 7,157 scaffolds with N50 of about 2 Mb. The repetitive sequences comprised of 41% of the total genome. BUSCO analysis revealed that the snow sheep assembly contained full-length or partial fragments of 97% of mammalian universal single-copy orthologs (n = 4,104), illustrating the completeness of the assembly. In addition, a total of 20,045 protein-coding sequences were identified using comprehensive gene prediction pipeline. Of which 19,240 (∼96%) sequences were annotated using protein databases. Moreover, homology-based searches and de novo identification detected 1,484 tRNAs; 243 rRNAs; 1,931 snRNAs; and 782 miRNAs in the snow sheep genome. To conclude, we generated the first de novo genome of the snow sheep using long reads; these data are expected to contribute significantly to our understanding related to evolution and adaptation within the Ovis genus.

Population Structure and Genetic Diversity of Sheep Breeds in the Kyrgyzstan.

Sheep are a main livestock species of Kyrgyzstan, a Central Asian country with predominating mountain terrain. The current gene pool of local sheep resources has been forming under diverse climate conditions from the era of the trading caravans of the Great Silk Road, through the Soviet period of large-scale livestock improvements, which was followed by the deep crisis at the end of the 20th century, up to now. However, not much is known about the genetic background and variability of the local sheep populations. Therefore, our aims were to provide a characterization of the population structure and genetic relations within the Kyrgyz sheep breeds and to study their genetic connections with the global sheep breeds using SNP analysis. Samples of the Alai (n = 31), Gissar (n = 30), Kyrgyz coarse wool (n = 13), Aykol (n = 31), and Tien-Shan (n = 24) breeds were genotyped with the OvineSNP50 BeadChip or the Ovine Infinium HD BeadChip (Illumina Inc., USA). The measure of inbreeding based on runs of homozygosity showed a minimum value in the Aykol breed (FROH = 0.034), while the maximum was found in the Alai breed (FROH = 0.071). Short ROH segments (ROH ≤ 4 Mb) were predominant in all breeds. Long ROH segments (ROH > 16 Mb) were absent in the Gissar breed. The Gissar and Aykol breeds had the highest values of the effective population sizes estimated for five generations ago (Ne5 = 660 and 563), whereas the Alai and Kyrgyz coarse wool displayed lower values (Ne5 = 176 and 128, respectively). The synthetic origin of the Aykol breed was clearly evidenced by all analyses applied. Based on the network and admixture analyses of the Kyrgyz and global sheep breeds, the Tien-Shan and the Russian semi-fine wool breeds demonstrated a common ancestry that most likely is due to a contribution of the Lincoln breed. The Gissar, Aykol, and Kyrgyz coarse wool breeds showed a genetic background predominating in sheep populations from Iran and China whereas the Alai demonstrated the different ancestry type. The revealed admixture patterns probably resulted from the exchange and trade during the era of the Great Silk Road, which partly overlapped with historical and archeological findings.

Links between maternal postpartum depressive symptoms, maternal distress, infant gender and sensitivity in a high-risk population.

BACKGROUND: Maternal postpartum depression has an impact on mother-infant interaction. Mothers with depression display less positive affect and sensitivity in interaction with their infants compared to non-depressed mothers. Depressed women also show more signs of distress and difficulties adjusting to their role as mothers than non-depressed women. In addition, depressive mothers are reported to be affectively more negative with their sons than with daughters. METHODS: A non-clinical sample of 106 mother-infant dyads at psychosocial risk (poverty, alcohol or drug abuse, lack of social support, teenage mothers and maternal psychic disorder) was investigated with EPDS (maternal postpartum depressive symptoms), the CARE-Index (maternal sensitivity in a dyadic context) and PSI-SF (maternal distress). The baseline data were collected when the babies had reached 19 weeks of age. RESULTS: A hierarchical regression analysis yielded a highly significant relation between the PSI-SF subscale "parental distress" and the EPDS total score, accounting for 55% of the variance in the EPDS. The other variables did not significantly predict the severity of depressive symptoms. A two-way ANOVA with "infant gender" and "maternal postpartum depressive symptoms" showed no interaction effect on maternal sensitivity. CONCLUSIONS: Depressive symptoms and maternal sensitivity were not linked. It is likely that we could not find any relation between both variables due to different measuring methods (self-reporting and observation). Maternal distress was strongly related to maternal depressive symptoms, probably due to the generally increased burden in the sample, and contributed to 55% of the variance of postpartum depressive symptoms.