DIA-Based Proteomic Analysis Reveals MYOZ2 as a Key Protein Affecting Muscle Growth and Development in Hybrid Sheep

Hybridization of livestock can be used to improve varieties, and different hybrid combinations produce unique breeding effects. In this study, male Southdown and Suffolk sheep were selected to hybridize with female Hu sheep to explore the effects of male parentage on muscle growth and the development of offspring. Using data-independent acquisition technology, we identified 119, 187, and 26 differentially abundant proteins (DAPs) between Hu × Hu (HH) versus Southdown × Hu (NH), HH versus Suffolk × Hu (SH), and NH versus SH crosses. Two DAPs, MYOZ2 and MYOM3, were common to the three hybrid groups and were mainly enriched in muscle growth and development-related pathways. At the myoblast proliferation stage, MYOZ2 expression decreased cell viability and inhibited proliferation. At the myoblast differentiation stage, MYOZ2 expression promoted myoblast fusion and enhanced the level of cell fusion. These findings provide new insights into the key proteins and metabolic pathways involved in the effect of male parentage on muscle growth and the development of hybrid offspring in sheep.

Genome-wide detection of structural variation in some sheep breeds using whole-genome long-read sequencing data

Genomic structural variants (SVs) constitute a significant proportion of genetic variation in the genome. The rapid development of long-reads sequencing has facilitated the detection of long-fragment SVs. There is no published study to detect SVs using long-read data from sheep. We applied a long-read mapping approach to detect SVs and characterized a total of 30,771 insertions, deletions, inversions and translocations. We identified 716, 916, 842 and 303 specific SVs in Southdown sheep, Alpine merino sheep, Qilian White Tibetan sheep and Oula sheep, respectively. We annotated these SVs and found that these SV-related genes were primarily enriched in the well-established pathways involved in the regulation of the immune system, growth and development and environmental adaptability. We detected and annotated SVs based on NGS resequencing data to validate the accuracy based on third-generation detection. Moreover, five candidate SVs were verified using the PCR method in 50 sheep. Our study is the first to use a long-reads sequencing approach to construct a novel structural variation map in sheep. We have completed a preliminary exploration of the potential effects of SVs on sheep.

Single-Cell Transcriptome Sequence Profiling on the Morphogenesis of Secondary Hair Follicles in Ordos Fine-Wool Sheep

The Ordos fine-wool sheep is a high-quality breed in China that produces superior natural textiles and raw materials such as wool and lamb meat. However, compared to the Australian Merino sheep, there is still a gap in terms of the wool fiber fineness and wool yield. The hair follicle is the main organ that controls the type of wool fiber, and the morphological changes in the secondary hair follicle are crucial in determining wool quality. However, the process and molecular mechanisms of hair follicle morphogenesis in Ordos fine-wool sheep are not yet clear. Therefore, analyzing the molecular mechanisms underlying the process of follicle formation is of great significance for improving the fiber diameter and wool production of Ordos fine-wool sheep. The differential expressed genes, APOD, POSTN, KRT5, and KRT15, which related to primary hair follicles and secondary hair follicles, were extracted from the dermal papillae. Based on pseudo-time analysis, the differentiation trajectories of dermal lineage cells and epidermal lineage cells in the Ordos fine-wool sheep were successfully constructed, providing a theoretical basis for breeding research in Ordos fine-wool sheep.

Genomic Selection for Live Weight in the 14th Month in Alpine Merino Sheep Combining GWAS Information

Alpine Merino Sheep is a novel breed reared from Australian Merino Sheep as the father and Gansu Alpine Fine-Wool Sheep as the mother, living all year in cold and arid alpine areas with exceptional wool quality and meat performance. Body weight is an important economic trait of the Alpine Merino Sheep, but there is limited research on identifying the genes associated with live weight in the 14th month for improving the accuracy of the genomic prediction of this trait. Therefore, this study's sample comprised 1310 Alpine Merino Sheep ewes, and the Fine Wool Sheep 50K Panel was used for genome-wide association study (GWAS) analysis to identify candidate genes. Moreover, the trial population (1310 ewes) in this study was randomly divided into two groups. One group was used as the population for GWAS analysis and screened for the most significant top 5%, top 10%, top 15%, and top 20% SNPs to obtain prior marker information. The other group was used to estimate the genetic parameters based on the weight assigned by heritability combined with different prior marker information. The aim of this study was to compare the accuracy of genomic breeding value estimation when combined with prior marker information from GWAS analysis with the optimal linear unbiased prediction method for genome selection (GBLUP) for the breeding value of target traits. Finally, the accuracy was evaluated using the five-fold cross-validation method. This research provides theoretical and technical support to improve the accuracy of sheep genome selection and better guide breeding. The results demonstrated that eight candidate genes were associated with GWAS analysis, and the gene function query and literature search results suggested that FAM184B, NCAPG, MACF1, ANKRD44, DCAF16, FUK, LCORL, and SYN3 were candidate genes affecting live weight in the 14th month (WT), which regulated the growth of muscle and bone in sheep. In genome selection analysis, the heritability of GBLUP to calculate the WT was 0.335-0.374, the accuracy after five-fold cross-verification was 0.154-0.190, and after assigning different weights to the top 5%, top 10%, top 15%, and top 20% of the GWAS results in accordance with previous information to construct the G matrix, the accuracy of the WT in the GBLUP model was improved by 2.59-7.79%.

Transcriptomics and metabolomics reveal improved performance of Hu sheep on hybridization with Southdown sheep

Consumers are increasingly demanding high-quality mutton. Cross breeding can improve meat quality and is widely used in sheep breeding. However, little is known about the molecular mechanism of cross breeding sheep meat quality. In this study, male Southdown and female Hu sheep were hybridized. The slaughter performance and longissimus dorsi quality of the 6-month-old hybrid offspring were measured, and the longissimus dorsi of the hybrid offspring was analyzed by transcriptomics and metabolomics to explore the effect of cross breeding on meat quality. The results showed that the production performance of Southdown × Hu F1 sheep was significantly improved, the carcass fat content was significantly decreased, and the eating quality of Southdown × Hu F1 sheep were better. Compared with the HS group (Hu × Hu), the NH group (Southdown × Hu) had 538 differentially expressed genes and 166 differentially expressed metabolites (P < 0.05), which were significantly enriched in amino acid metabolism and other related pathways. Up-regulated genes METTL21C, PPARGC1A and down-regulated gene WFIKKN2 are related to muscle growth and development. Among them, the METTL21C gene, which is related to muscle development, was highly correlated with carnosine, a metabolite related to meat quality (correlation > 0.6 and P < 0.05). Our results provide further understanding of the molecular mechanism of cross breeding for sheep muscle growth and meat quality optimization.

Spatial and temporal distribution of endotoxins, antibiotic resistance genes and mobile genetic elements in the air of a dairy farm in Germany

Antimicrobial resistance (AMR) is a serious issue that is continuously growing and spreading, leading to a dwindling number of effective treatments for infections that were easily treatable with antibiotics in the past. Animal farms are a major hotspot for AMR, where antimicrobials are often overused, misused, and abused, in addition to overcrowding of animals. In this study, we investigated the risk of AMR transmission from a farm to nearby residential areas by examining the overall occurrence of endotoxins, antibiotic resistance genes (ARGs), and mobile genetic elements (MGEs) in the air of a cattle farm. We assessed various factors, including the season and year, day and nighttime, and different locations within the farm building and its vicinity. The most abundant ARGs detected were tetW, aadA1, and sul2, genes that encode for resistances towards antibiotics commonly used in veterinary medicine. While there was a clear concentration gradient for endotoxin from the middle of the farm building to the outside areas, the abundance of ARGs and MGEs was relatively uniform among all locations within the farm and its vicinity. This suggests that endotoxins preferentially accumulated in the coarse particle fraction, which deposited quickly, as opposed to the ARGs and MGEs, which might concentrate in the fine particle fraction and remain longer in the aerosol phase. The occurrence of the same genes found in the air samples and in the manure indicated that ARGs and MGEs in the air mostly originated from the cows, continuously being released from the manure to the air. Although our atmospheric dispersion model indicated a relatively low risk for nearby residential areas, farm workers might be at greater risk of getting infected with resistant bacteria and experiencing overall respiratory tract issues due to continuous exposure to elevated concentrations of endotoxins, ARGs and MGEs in the air of the farm.

Pulse Oximetry and Perfusion Index Screening for Congenital Heart Defects: A Systematic Review and Meta-analysis

Congenital heart defects (CHDs) are the most common neonatal malformations and are a leading cause of infant death in developed countries. Finding safe and effective diagnostic methods to screen for CHDs is important. The aim of this study was to evaluate the effectiveness of pulse oximetry (PO) and perfusion index (PI) in screening CHD. We conducted a systematic review of studies in PubMed, Embase, and the Cochrane Library published on or before October 1, 2021. Studies based on PICOS were included in this systematic review. The flow chart is made by PRISMA software. The quality of included studies was assessed by RevMan5 software (QUADAS-2: Quality Assessment of Diagnostic Accuracy Studies-2). The sensitivity, specificity, and other measurements of accuracy were pooled using Stata/SE 12.0 software. Five studies containing 46,965 neonates were included in this study. A randomized-effects model was used for the meta-analysis because of significant heterogeneity. The combined sensitivity and specificity were 0.82 (95% confidence interval [CI], 0.53-0.95) and 0.97 (95% CI, 0.57-1.00), respectively. The area under the curve was 0.92 (95% CI, 0.89-0.94). The combination PO and PI was significant in CHD screening. Once diagnosed by the combined method, it means that the neonate is most likely to have a CHD. KEY POINTS: · Pulse oximetry and PI screening.. · Congenital heart defects.. · A systematic review and meta-analysis..

Global DNA Methylation, miRNA, and mRNA Profiles in Sheep Skeletal Muscle Promoted by Hybridization

With the development of high-throughput sequencing technology, several nongenetic variations, including noncoding RNAs such as miRNAs, and DNA methylation, have been found to play an important role in animal muscle development and fat metabolism. In this study, Southdown and Suffolk were selected as male parents for hybridization with Hu sheep (Southdown × Hu (NH), Suffolk × Hu (SH), and Hu × Hu (HH)). RNA sequencing, bisulfite sequencing, and small-RNA sequencing were used to study the methylation patterns and differences in miRNA and mRNA expression in the F1 sheep longissimus dorsi muscle tissue. We identified 765 differentially expressed genes (DEGs), 10,161 differentially methylated regions (DMRs), and 164 differentially expressed miRNAs, which were significantly enriched in AMPK signaling, fatty acid degradation, metabolism, and other related pathways (P < 0.05). In addition, we constructed a DNA methylation-mRNA and miRNA-mRNA coexpression network. A total of 42 common genes were identified from DMRs and DEGs. Importantly, we predicted that 33 differentially expressed miRNAs directly or indirectly targeted the SLC27A6. The data obtained in this study provide useful information and evidence to support further understanding of the miRNA and DNA methylation of key genes regulating muscle growth and fat metabolism in hybrid sheep populations.

Dose Coverage Variation Caused by Setup Uncertainties in VMAT-TBI Treatment

PURPOSE/OBJECTIVE(S)

Volumetric Modulated Arc Therapy Total Body Irradiation (VMAT TBI) offers several advantages over conventional TBI techniques, including reduced cost of the treatment room (a normal-sized vault versus a large shielded vault), increased patient comfort during treatment, lung sparing without the need for physical blocks, and a more homogeneous dose distribution. In VMAT-TBI treatment, plans at multiple isocenters are utilized to cover the entire body, and patients are positioned with the aid of image guidance. However, aligning the patient's entire body during setup can be challenging due to its pliability, and the setup time is heavily dependent on the tolerance allowed by the image guidance. Therefore, studying the variation in dose coverage caused by setup uncertainties in VMAT-TBI treatments can help optimize the clinical workflow and determine the optimal tolerance for patient positioning.

MATERIALS/METHODS

New plans were generated to simulate the uncertainties that occur during treatment setup for each patient. These plans were created by shifting the original VMAT TBI plans at the head, chest, abdomen (or pelvis) isocenters by 5mm and 1cm in the left-right (LR), inferior-superior (IS), and anterior-posterior (AP) directions, respectively. Dose DicomRT files were exported, and the dose change due to the shifts was analyzed. The statistical quantification of the percentage of the body that experienced a dose change of over 2%, 5%, and 10% of the prescription due to the shifts from the original plans was calculated for all patients. Histograms were generated, showing the percentage of body getting dose change of 1-2%, 2-3%, 3-5%, 5-7%, 7-10%, 10-15%.

RESULTS

The table below displays the percent volume receiving a dose change of 2%, 5%, and 10% of the prescription for a 5mm shift. Among the shift directions, the dose change is most sensitive in the IS direction, with similar impact observed the in LR and AP directions. Among different sites, the chest experiences the largest dose change, followed by the pelvis. For a 5mm shift in the IS direction, the average percent volume receiving a dose change of 2%, 5%, and 10% in the chest is 9.25%, 2.64%, and 0.27%, respectively. For a 1cm shift, the numbers are 12.23%, 6.75%, and 1.29%. In the pelvis (abdomen), these values are 9.03%, 1.67%, and 0.17% for a 5mm shift and 13.28%, 6.1%, and 0.85% for a 1cm shift. For head plans, the values are 2.72%, 0.9%, and 0.14% for a 5mm shift and 3.77%, 1.66%, and 0.53% for a 1cm shift.

CONCLUSION

Accurate alignment in the chest region is crucial in VMAT TBI treatment. Efforts should be made to minimize shifts over 1cm in the IS direction.

Dosimetric Impact of Seed Segmentation in GammaTile Surgically Targeted Radiation Therapy for Gliomas and Brain Metastases Cases

PURPOSE/OBJECTIVE(S)

To assess the variability of dose reporting variability due to uncertainty in segmentation of Cs-131 seeds in GammaTile therapy for gliomas and brain metastases.

MATERIALS/METHODS

Ten patients with either glioma or brain metastases had 4-11 GammaTiles placed along resection bed during craniotomy. A dose of 60 Gy is prescribed to 5 mm depth. Each GammaTile has four Cs-131 seeds imbedded in a biodegradable collagen sponge. GammaTile Post-Op workflow in MIM Symphony software is used for post-implant dose evaluation and reporting. This workflow requires a post-surgery CT to identify seeds, and a post-surgery MR for residual disease and OAR contours. Seeds are segmented using a threshold tool. Threshold levels may change depending on the CT used, thus users need to manually change the HU threshold value in each data set. Since GammaTiles are lined along the resection bed, PTVs are generated automatically by adding 8 mm expansion on the seed contours and later combined with residual disease contours. We simulate the seed contour uncertainty by applying -0.5 mm, -1.0 mm, +0.5 mm & +1.0 mm concentric margins to the current seed contours to create 4 new seed contours per patient. New PTVs are generated by adding 8 mm expansion on the new seed contours combined with residual disease contours. PTV volume, PTV volume receiving 100% and 150% of prescription dose (V100, V150), and percentage of the prescription dose received by 90% of the PTV (D90) are calculated to evaluate dose reporting variability due to seed segmentation uncertainty.

RESULTS

Mean PTV volume decreases by 8.4 cc & 10.2 cc for PTVs generated from seed contours with -0.5 mm & -1.0 mm margin, respectively, and increases by 5.8 cc & 8.2 cc, respectively, when +0.5 mm & +1 mm margins are applied to the original seed contours. We observe up to 10% change in V100 due to seed segmentation uncertainty. Mean V100 increases by 4.0% (range: 0.2% - 8.9%) & 4.9% (range: 0.5% - 11.0%) for cases with -0.5 m & -1.0 mm seed margin, respectively, and reduces by 4.2% (range: 0.5% - 6.7%) & 5.9% (range: 0.6% - 10.4%) for cases with +0.5 mm and +1.0 mm seed margin, respectively. Mean D90 increases by 7.7% (range: 4.0% - 12.6%) & 9.9% (range: 4.0% - 17.4%) for cases with -0.5 m & -1.0 mm seed margin, respectively, and reduces by 5.5% (range: 3.6% - 7.8%) & 7.4% (range: 5.2% - 9.6%) for cases with +0.5 mm and +1.0 mm seed margin, respectively. We also observe up to 8.0% changes in mean V150 when margins are applied to the seed contours.

CONCLUSION

Our results show significant impact of seed segmentation uncertainty on dose reporting in GammaTile therapy. Variability in dose reporting parameters highlight the need for a more standardized and automated approach to seed segmentation to ensure consistent and accurate dose reporting. The current manual threshold adjustment method is subject to user dependence and therefore unreliable. Development of a more robust tool could help to minimize variability and improve reliability of dose reporting.

Transcriptome-metabolome analysis reveals how sires affect meat quality in hybrid sheep populations

Crossbreeding improves and enhances meat quality and is widely used in sheep production; however, the molecular mechanisms underlying the meat quality of various crossbred sheep remain unknown. In this study, male Southdown, Suffolk and Hu sheep were crossbred with female Hu sheep, and the transcriptomes and metabolomes of the longissimus dorsi muscle of the F1 generation were sequenced to explore how different sire breeds affect meat quality. The results showed that 631 differentially expressed genes and 119 significantly altered metabolites contributed to muscle development characteristics and meat quality-related diversity (P < 0.05). These genes and metabolites were significantly enriched in lipid metabolism pathways, including arachidonic acid metabolism and PPAR signaling. Several candidate genes were associated with muscle growth, such as MYLK3, MYL10, FIGN, MYH8, MYOM3, LMCD1, and FLRT1. Among these, MYH8 and MYL10 participated in regulating muscle growth and development and were correlated with meat quality-related fatty acid levels (|r| > 0.5 and p < 0.05). We selected mRNA from four of these genes to verify the accuracy of the sequencing data via qRT-PCR. Our findings provide further insight into the key genes and metabolites involved in muscle growth and meat quality in hybrid sheep populations.

Evaluation of Mutton Quality Characteristics of Dongxiang Tribute Sheep Based on Membership Function and Gas Chromatography and Ion Mobility Spectrometry

Dongxiang tribute sheep have a history of use in food dishes such as "Dongxiang Handgrip," which dates back hundreds of years and is a favorite halal food in northwestern China. However, little is known about the mutton quality characteristics of Dongxiang tribute sheep. Here, we measured the sensory characteristics, nutritional quality, and flavor substances to comprehensively evaluate the mutton quality characteristics of these sheep. The mutton qualities of Dongxiang tribute, Tibetan, Ujumqin, and Hu sheep were comprehensively evaluated by membership function. Subsequently, the volatile components in mutton samples from 30 Dongxiang tribute sheep were detected via gas chromatography and ion mobility spectrometry (GC-IMS), and their fingerprints were established. The result of meat quality revealed that the shear force, the contents of protein, essential amino acid (EAA), non-essential amino acid (NEAA), and n-6/n-3 ratio of Dongxiang tribute mutton were better than the other three breeds. Membership functions were calculated for 10 physical and chemical indexes of mutton quality, and the comprehensive membership function values of the four breeds in order of highest to lowest mutton quality were Tibetan sheep (0.76) > Dongxiang tribute sheep (0.49) > Hu sheep (0.46) > Ujumqin sheep (0.33). Thirty volatile compounds were identified via GC-IMS: seven alcohols, eight aldehydes, five ketones, two esters, two phenols, one ether, one furan, one acid, two hydrocarbons, and one pyrazine. Ketones, aldehydes, and alcohols were the main volatile compounds forming the flavor of Dongxiang tribute sheep mutton. The reliability of the results was validated by PCA (principal component analysis) and similarity analyses. Our results provide reference value for consumers of mutton in China.

Genetic Basis of Dorper Sheep (Ovis aries) Revealed by Long-Read De Novo Genome Assembly

Dorper sheep (Ovis aries) (DPS), developed in the 1930s by crossing Dorset Horn and Blackhead Persian sheep in South Africa, is a world-famous composite breed for mutton production. The genetic basis underlying this breed is yet to be elucidated. Here, we report the sequencing and assembly of a highly contiguous Dorper sheep genome via integration of Oxford Nanopore Technology (ONT) sequencing and Hi-C (chromatin conformation capture) approaches. The assembled genome was around 2.64 Gb with a contig N50 of 73.33 Mb and 140 contigs in total. More than 99.5% of the assembled sequences could be anchored to 27 chromosomes and they were annotated with 20,450 protein-coding genes. Allele-specific expression (ASE) genes of Dorper sheep were revealed through ASE analysis and they were involved in the immune system, lipid metabolism, and environmental adaptation. A total of 5,701 and 456 allelic sites were observed in the SNP and indels loci identified from relevant whole-genome resequencing data. These allelic SNP and INDEL sites were annotated in 1,002 and 294 genes, respectively. Moreover, we calculated the number of variant sites and related genes derived from the maternal and paternal ancestors, revealing the genetic basis of outstanding phenotypic performance of Dorper sheep. In conclusion, this study reports the first reference genome of Dorper sheep and reveals its genetic basis through ASE. This study also provides a pipeline for mining genetic information of composite breeds, which has an implication for future hybrid-breeding practices.

Comprehensive Analysis of Long Non-coding RNA and mRNA Transcriptomes Related to Hypoxia Adaptation in Tibetan Sheep

Tibetan sheep have lived on the Qinghai-Tibet Plateau for a long time, and after long-term natural selection, they have shown stable genetic adaptability to high-altitude environments. However, little is known about the molecular mechanisms of the long non-coding (lnc)RNAs involved in the adaptation of Tibetan sheep to hypoxia. Here, we collected lung tissues from high-altitude Tibetan sheep and low-altitude Hu sheep for RNA sequencing to study the regulatory mechanisms of the lncRNAs and mRNAs in the adaptation of Tibetan sheep to hypoxia. We identified 254 differentially expressed lncRNAs and 1,502 differentially expressed mRNAs. We found 20 pairs of cis-regulatory relationships between 15 differentially expressed lncRNAs and 14 protein-coding genes and two pairs of trans-regulatory relationships between two differentially expressed lncRNAs and two protein-coding genes. These differentially expressed mRNAs and lncRNA target genes were mainly enriched in pathways related to lipid metabolism and immune function. Interaction network analysis showed that 17 differentially expressed lncRNAs and 15 differentially expressed mRNAs had an interactive relationship. Additionally, we used six differentially expressed lncRNAs and mRNAs to verify the accuracy of the sequencing data via qRT-PCR. Our results provide a comprehensive overview of the expression patterns of the lncRNAs and mRNAs involved in the adaptation of Tibetan sheep to hypoxia, laying a foundation for further analysis of the adaptations of plateau animals.

The Transcriptional Cell Atlas of Testis Development in Sheep at Pre-Sexual Maturity

Sheep testes undergo a dramatic rate of development with structural changes during pre-sexual maturity, including the proliferation and maturation of somatic niche cells and the initiation of spermatogenesis. To explore this complex process, 12,843 testicular cells from three males at pre-sexual maturity (three-month-old) were sequenced using the 10× Genomics Chromium^TM single-cell RNA-seq (scRNA-seq) technology. Nine testicular somatic cell types (Sertoli cells, myoid cells, monocytes, macrophages, Leydig cells, dendritic cells, endothelial cells, smooth muscle cells, and leukocytes) and an unknown cell cluster were observed. In particular, five male germ cell types (including two types of undifferentiated spermatogonia (A_pale and A_dark), primary spermatocytes, secondary spermatocytes, and sperm cells) were identified. Interestingly, A_pale and A_dark were found to be two distinct states of undifferentiated spermatogonia. Further analysis identified specific marker genes, including UCHL1, DDX4, SOHLH1, KITLG, and PCNA, in the germ cells at different states of differentiation. The study revealed significant changes in germline stem cells at pre-sexual maturation, paving the way to explore the candidate factors and pathways for the regulation of germ and somatic cells, and to provide us with opportunities for the establishment of livestock stem cell breeding programs.

Genome-Wide Selective Signatures Reveal Candidate Genes Associated with Hair Follicle Development and Wool Shedding in Sheep

Hair follicle development and wool shedding in sheep are poorly understood. This study investigated the population structures and genetic differences between sheep with different wool types to identify candidate genes related to these traits. We used Illumina ovine SNP 50K chip genotyping data of 795 sheep populations comprising 27 breeds with two wool types, measuring the population differentiation index (Fst), nucleotide diversity (θπ ratio), and extended haplotype homozygosity among populations (XP-EHH) to detect the selective signatures of hair sheep and fine-wool sheep. The top 5% of the Fst and θπ ratio values, and values of XP-EHH < −2 were considered strongly selected SNP sites. Annotation showed that the PRX, SOX18, TGM3, and TCF3 genes related to hair follicle development and wool shedding were strongly selected. Our results indicated that these methods identified important genes related to hair follicle formation, epidermal differentiation, and hair follicle stem cell development, and provide a meaningful reference for further study on the molecular mechanisms of economically important traits in sheep.

Evaluation of Bayesian alphabet and GBLUP based on different marker density for genomic prediction in Alpine Merino sheep

The marker density, the heritability level of trait and the statistical models adopted are critical to the accuracy of genomic prediction (GP) or selection (GS). If the potential of GP is to be fully utilized to optimize the effect of breeding and selection, in addition to incorporating the above factors into simulated data for analysis, it is essential to incorporate these factors into real data for understanding their impact on GP accuracy, more clearly and intuitively. Herein, we studied the GP of six wool traits of sheep by two different models, including Bayesian Alphabet (BayesA, BayesB, BayesCπ, and Bayesian LASSO) and genomic best linear unbiased prediction (GBLUP). We adopted fivefold cross-validation to perform the accuracy evaluation based on the genotyping data of Alpine Merino sheep (n = 821). The main aim was to study the influence and interaction of different models and marker densities on GP accuracy. The GP accuracy of the six traits was found to be between 0.28 and 0.60, as demonstrated by the cross-validation results. We showed that the accuracy of GP could be improved by increasing the marker density, which is closely related to the model adopted and the heritability level of the trait. Moreover, based on two different marker densities, it was derived that the prediction effect of GBLUP model for traits with low heritability was better; while with the increase of heritability level, the advantage of Bayesian Alphabet would be more obvious, therefore, different models of GP are appropriate in different traits. These findings indicated the significance of applying appropriate models for GP which would assist in further exploring the optimization of GP.

[Analysis on the imported Coronavirus Disease 2019 related cluster epidemic in rural areas of Chengdu]

An epidemiological investigation was carried out on a local cluster of outbreak caused by imported cases of Coronavirus Disease 2019 (COVID-19) in rural areas of Chengdu in December 2020, to find out the source of infection and the chain of transmission. According to Prevention and Control Protocol for COVID-19 (Version 7), field epidemiological investigation was adopted, combined with big data technology, video image investigation, gene sequencing and other methods to carry out investigation into COVID-19 cases and infections source tracing, analyze the epidemiological association, and map the chain of transmission. From December 7 to 17, 2020, 13 local COVID-19 confirmed cases and 1 asymptomatic case were diagnosed in Chengdu, of which 12 cases (85.71%) had a history of residence and activity in the village courtyard of Taiping (TP), Pidu (P) District, Chengdu. From November 8, 2020 to November 28, 2020, a group of inbound people form Nepal were transferred to the designated entry personnel quarantine hotel of P District which was adjacent to the TP village. During quarantine, there were 5 cases who tested positive for COVID-19. Through gene sequencing alignment, genes of local cases and Nepalese imported cases from the same period are homologous, all belong to the lineage of L2.2.3 (B.1.36 according to Pangolin lineage typing method). According to the results of field epidemiological investigation and gene sequencing analysis, the index case was most likely infected by contact with household waste of quarantine site. Under the situation of normalization prevention and control of COVID-19, sentinel monitoring of fever clinics in primary medical institutions is the key to early detection of the epidemic. The multi-department joint epidemiological investigation and the application of gene technology are the core links of the investigation and traceability of modern infectious diseases. The allocation of public health resources in rural areas needs to be strengthened. We need to improve the capacity for early surveillance and early warning of the epidemic in rural areas.

A survey of a COVID-19 cluster of charter flight importation

OBJECTIVES

Although a number of cases of importation with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection have been reported, there are still no data available concerning the characteristics in the coronavirus disease 2019 (COVID-19) cluster of charter flight importation. Here, we provide an analysis of COVID-19 cases and their close contacts who worked for the same company on a project in Karbala, Iraq, and returned back to Chengdu, China, by a charter flight.

METHODS

The data of imported COVID-19 cases and their close contacts were obtained from National Notifiable Disease Report System of Chinese Center for Disease Control and Prevention and field epidemiological investigation reports by Centers for Disease Control and Prevention (CDCs) in Chengdu. The information of general characteristics and laboratory findings of this cluster were collected and summarized.

RESULTS

One hundred and six (66.67%) of 159 charter flight passengers tested positive for COVID-19 before entry. Through treatment, all 159 people tested negative and meet the requirements of taking flights bound for China before boarding. However, there has been still 36 (22.64%) of them tested positive after entry. The median time from entry to confirmation was 1.0 day (Interquartile Range (IQR): 0-4.3). The Cycle threshold value (Ct value) of 36 patients' positive samples are all above 30 and most values are above 35.

CONCLUSIONS

In conclusion, there is still a risk that a number of COVID-19 cases can be imported through charter flight. However, the infectivity of confirmed patients of the charter flight was considered to be low.

Whole-genome re-sequencing association study on yearling wool traits in Chinese fine-wool sheep

To investigate single nucleotide polymorphism (SNP) loci associated with yearling wool traits of fine-wool sheep for optimizing marker-assisted selection and dissection of the genetic architecture of wool traits, we conducted a genome-wide association study (GWAS) based on the fixed and random model circulating probability unification (FarmCPU) for yearling staple length (YSL), yearling mean fiber diameter (YFD), yearling greasy fleece weight (YGFW), and yearling clean fleece rate (YCFR) by using the whole-genome re-sequenced data (totaling 577 sheep) from the following four fine-wool sheep breeds in China: Alpine Merino sheep (AMS), Chinese Merino sheep (CMS), Qinghai fine-wool sheep (QHS), and Aohan fine-wool sheep (AHS). A total of 16 SNPs were detected above the genome-wise significant threshold (P = 5.45E-09), and 79 SNPs were located above the suggestive significance threshold (P = 5.00E-07) from the GWAS results. For YFD and YGFW traits, 7 and 9 SNPs reached the genome-wise significance thresholds, whereas 10 and 12 SNPs reached the suggestive significance threshold, respectively. For YSL and YCFR traits, none of the SNPs reached the genome-wise significance thresholds, whereas 57 SNPs exceeded the suggestive significance threshold. We recorded 14 genes located at the region of ±50-kb near the genome-wise significant SNPs and 59 genes located at the region of ±50-kb near the suggestive significant SNPs. Meanwhile, we used the Average Information Restricted Maximum likelihood algorithm (AI-REML) in the “HIBLUP” package to estimate the heritability and variance components of the four desired yearling wool traits. The estimated heritability values (h²) of YSL, YFD, YGFW, and YCFR were 0.6208, 0.7460, 0.6758, and 0.5559, respectively. We noted that the genetic parameters in this study can be used for fine-wool sheep breeding. The newly detected significant SNPs and the newly identified candidate genes in this study would enhance our understanding of yearling wool formation, and significant SNPs can be applied to genome selection in fine-wool sheep breeding.

[Influence of SARS-CoV-2 vaccination on the epidemiological and clinical characteristics of imported COVID-19 cases in Chengdu]

Objective: To analyze the epidemiological and clinical characteristics of imported COVID-19 cases after SARS-CoV-2 vaccination and to provide evidence for the prevention and control of COVID-19. Methods: The imported COVID-19 cases in Chengdu as of April 15, 2021 were divided into the vaccinated group and unvaccinated group according to the history of SARS-CoV-2 vaccination. The epidemiological and clinical data of the cases were collected retrospectively, and the differences in epidemiological and clinical characteristics of the two groups were compared. Laboratory tests consisted of nucleic acid test, clinical index test, serum antibody test and lymphocyte test. Software WPS2019 was used for data management and software R 4.0.3 was used for statistical analysis. Results: A total of 75 COVID-19 cases were included in the analysis, in which 20 had received SARS-CoV-2 vaccination and only 4 with clinical symptoms, 55 patients did not receive SARS-CoV-2 vaccination, and 16 had clinical symptoms. In vaccinated group, the first injection time of vaccination ranged from July to November 2020, and 10 cases received two doses of vaccine simultaneously and 10 cases received two doses of vaccine at intervals of 14-57 days. The intervals between the completion of vaccination and the onset ranged from 87 days to 224 days. The differences in classification and clinical type between the two groups were significant. Significant differences were observed in case classification and clinical type between vaccinated group and unvaccinated group (P<0.05). The vaccinated group had a relatively high proportion of asymptomatic infections (40.00%, 8/20), while mild infections were mainly observed in the unvaccinated group(76.36%,42/55). The differences in Ct values (ORF1ab gene and N gene) at the diagnosis were not significant between vaccinated group and unvaccinated group (P>0.05), similar results were also observed in lymphocyte subtypes, procalcitonin and C-reactive protein level comparisons. Serum amyloid A level was higher in unvaccinated group than in vaccinated group (P<0.05). However, the SARS-CoV-2 related serum antibody of IgM, IgG and total antibody levels were significantly higher in vaccinated group (P<0.05). Conclusions: Risk of infection still exists with SARS-CoV-2 after vaccination, which can facilitate the production of specific serum antibody of IgM and IgG when people are exposed to the virus. It has a certain protective effect on SARS-CoV-2 infected persons. Vaccination can reduce the clinical symptoms and mitigate disease severity.

Multiple-camera defocus imaging of ultracold atomic gases

In cold atom experiments, each image of light refracted and absorbed by an atomic ensemble carries a remarkable amount of information. Numerous imaging techniques including absorption, fluorescence, and phase-contrast are commonly used. Other techniques such as off-resonance defocused imaging (ORDI, [1-4]), where an in-focus image is deconvolved from a defocused image, have been demonstrated but find only niche applications. The ORDI inversion process introduces systematic artifacts because it relies on regularization to account for missing information at some spatial frequencies. In the present work, we extend ORDI to use multiple cameras simultaneously at degrees of defocus, eliminating the need for regularization and its attendant artifacts. We demonstrate this technique by imaging Bose-Einstein condensates, and show that the statistical uncertainties in the measured column density using the multiple-camera off-resonance defocused (McORD) imaging method are competitive with absorption imaging near resonance and phase contrast imaging far from resonance. Experimentally, the McORD method may be incorporated into existing set-ups with minimal additional equipment.

Prospective evaluation of dietary and lifestyle pattern indices with risk of colorectal cancer in a cohort of younger women

BACKGROUND

Although colorectal cancer (CRC) incidence in the USA is declining overall, its incidence is increasing among those younger than 50 years of age. The reasons underlying the increasing trend are largely unknown, although behavioral changes, such as unhealthy diet and lifestyle factors, may be partially responsible.

DESIGN

A prospective cohort study included 94 217 women aged 26-45 years at baseline. Validated anthropometric measures and lifestyle information were self-reported biennially. Exposures were four recommendation-based dietary indices-the prime diet quality score and three plant-based dietary indices; and two mechanism-based indices-the empirical dietary and lifestyle index for hyperinsulinemia (EDIH and ELIH). We calculated hazard ratios (HRs) and 95% confidence intervals (CIs) for overall CRC and for early-onset (before age 50) and after age 50 CRC separately.

RESULTS

We documented 332 cases of CRC during 24 years of follow-up (2 113 655 person-years), with an average age of 52 ± 7 years at diagnosis. Hyperinsulinemic dietary and lifestyle patterns were associated with a higher risk of CRC. Multivariable-adjusted HRs (95% CIs) comparing participants in the highest versus lowest quartile were: 1.67 for EDIH (95% CI: 1.15-2.44, P-trend = 0.01) and 1.51 for ELIH (95% CI: 1.10-2.08, P-trend = 0.01). Moreover, per 75% increment in rank, ELIH appeared to be a stronger risk factor for early-onset CRC (HR = 1.86, 95% CI: 1.12-3.07) than after age 50 CRC (HR = 1.20, 95% CI: 0.83-1.73, P-heterogeneity = 0.16). The four recommendation-based indices were not significantly associated with overall, early-onset, or after age 50 CRC risk (per 75% increment in rank, HRs ranged from 0.75 to 1.28).

CONCLUSION

Dietary and lifestyle patterns contributing to hyperinsulinemia were associated with greater CRC risk in younger women. Moreover, the hyperinsulinemic lifestyle showed a suggestively stronger positive association with early-onset CRC risk, compared with after age 50 CRC. Our findings suggest that dietary and lifestyle interventions to reduce insulinemic potential may be effective for CRC prevention among younger women.

Superimposing hip extension on knee flexion evokes higher activation in biceps femoris than knee flexion alone

Hamstring muscle function during knee flexion has been linked to hamstring injury and performance. However, it is unclear whether knee flexion alone (KF) requires similar hamstring electromyography (EMG) activity pattern to simultaneous hip extension and knee flexion (HE-KF), a combination that occurs in the late swing phase of sprinting. This study examined whether HE-KF maximal voluntary isometric contraction (MVIC) evokes higher (EMG) activity in biceps femoris long head (BFlh) and semitendinosus (ST) than KF alone. Effects of shank rotation angles were also tested. Twenty-one males performed the above-mentioned MVICs while EMG activity was measured along ST and BFlh. Conditions were compared using a one-way mixed functional ANOVA model under a fully Bayesian framework. Higher EMG activity was found in HE-KF in all shank rotation positions than in KF in the middle region of BFlh (highest in the 9th channel, by 0.022 mV [95%CrI 0.014 to 0.030] in neutral shank position). For ST, this was only observed in the neutral shank position and in the most proximal channel (by 0.013 mV [95%CrI 0.001 to 0.025]). We observed muscle- and region-specific responses to HE-KF. Future studies should examine whether hamstring activation in this task is related to injury risk and sprint performance.

Selective Sweeps Uncovering the Genetic Basis of Horn and Adaptability Traits on Fine-Wool Sheep in China

Long-term natural and artificial selection leads to change in certain regions of the genome, resulting in selection signatures that can reveal genes associated with selected traits, such as horns (i.e., polled/horned), high-quality wool traits, and high-altitude hypoxia adaptability. These are complex traits determined by multiple genes, regulatory pathways, and environmental factors. A list of genes with considerable effects on horn and adaptability traits has not been found, although multiple quantitative trait loci (QTL) have been identified. Selection signatures could be identified using genetic differentiation (F_ST), polymorphism levels θπ, and Tajima’s D. This study aimed to identify selection signatures in fine-wool sheep and to investigate the genes annotated in these regions, as well as the biological pathways involved in horn and adaptability traits. For this purpose, the whole-genome sequence of 120 individuals from four breeds, which come from different elevations and habitats in China, was used to analyze selection signatures for horn and adaptability traits. Annotation of the consensus regions of F_ST and θπ ratios revealed a list of identified genes associated with polled/horned and high-altitude hypoxia adaptability traits, such as RXPF2, EERFC4, MSH6, PP1R12A, THBS1, ATP1B2, RYR2, and PLA2G2E. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis identified genes related primarily to mismatch repair, metabolism, vascular smooth muscle contraction, and cardiac muscle contraction. This is the first study to demonstrate that selection signatures play an important role in the polled/horned and high-altitude hypoxia adaptability traits of fine-wool sheep breeds that have undergone high-intensity selection and adapted to different ecological environments in China. Changes observed in the genome of fine-wool sheep may have acted on genomic regions that affect performance traits and provide a reference for genome design and breeding.

Genome-wide association studies detects candidate genes for wool traits by re-sequencing in Chinese fine-wool sheep

BACKGROUND

The quality and yield of wool determine the economic value of the fine-wool sheep. Therefore, discovering markers or genes relevant to wool traits is the cornerstone for the breeding of fine-wool sheep. In this study, we used the Illumina HiSeq X Ten platform to re-sequence 460 sheep belonging to four different fine-wool sheep breeds, namely, Alpine Merino sheep (AMS), Chinese Merino sheep (CMS), Aohan fine-wool sheep (AHS) and Qinghai fine-wool sheep (QHS). Eight wool traits, including fiber diameter (FD), fiber diameter coefficient of variance (FDCV), fiber diameter standard deviation (FDSD), staple length (SL), greasy fleece weight (GFW), clean wool rate (CWR), staple strength (SS) and staple elongation (SE) were examined. A genome-wide association study (GWAS) was performed to detect the candidate genes for the eight wool traits.

RESULTS

A total of 8.222 Tb of raw data was generated, with an average of approximately 8.59X sequencing depth. After quality control, 12,561,225 SNPs were available for analysis. And a total of 57 genome-wide significant SNPs and 30 candidate genes were detected for the desired wool traits. Among them, 7 SNPs and 6 genes are related to wool fineness indicators (FD, FDCV and FDSD), 10 SNPs and 7 genes are related to staple length, 13 SNPs and 7 genes are related to wool production indicators (GFW and CWR), 27 SNPs and 10 genes associated with staple elongation. Among these candidate genes, UBE2E3 and RHPN2 associated with fiber diameter, were found to play an important role in keratinocyte differentiation and cell proliferation. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results, revealed that multitude significant pathways are related to keratin and cell proliferation and differentiation, such as positive regulation of canonical Wnt signaling pathway (GO:0090263).

CONCLUSION

This is the first GWAS on the wool traits by using re-sequencing data in Chinese fine-wool sheep. The newly detected significant SNPs in this study can be used in genome-selective breeding for the fine-wool sheep. And the new candidate genes would provide a good theoretical basis for the fine-wool sheep breeding.

Quantitative proteomic analysis identified differentially expressed proteins with tail/rump fat deposition in Chinese thin- and fat-tailed lambs

Tail adipose as one of the important functional tissues can enhance hazardous environments tolerance for sheep. The objective of this study was to gain insight into the underlying development mechanisms of this trait. A quantitative analysis of protein abundance in ovine tail/rump adipose tissue was performed between Chinese local fat- (Kazakh, Hu and Lanzhou) and thin-tailed (Alpine Merino, Tibetan) sheep in the present study by using lable-free approach. Results showed that 3400 proteins were identified in the five breeds, and 804 were differentially expressed proteins, including 638 up regulated proteins and 83 down regulated proteins in the tail adipose tissues between fat- and thin-tailed sheep, and 8 clusters were distinguished for all the DEPs’ expression patterns. The differentially expressed proteins are mainly associated with metabolism pathways and peroxisome proliferator activated receptor signaling pathway. Furthermore, the proteomics results were validated by quantitative real-time PCR and Western Blot. Our research has also suggested that the up-regulated proteins ACSL1, HSD17β4, FABP4 in the tail adipose tissue might contribute to tail fat deposition by facilitating the proliferation of adipocytes and fat accumulation in tail/rump of sheep. Particularly, FABP4 highly expressed in the fat-tail will play an important role for tail fat deposition. Our study might provide a novel view to understanding fat accumulation in special parts of the body in sheep and other animals.

A global analysis of CNVs in Chinese indigenous fine-wool sheep populations using whole-genome resequencing

Background

Copy number variation (CNV) is an important source of genetic variation that has a significant influence on phenotypic diversity, economically important traits and the evolution of livestock species. In this study, the genome-wide CNV distribution characteristics of 32 fine-wool sheep from three breeds were analyzed using resequencing.

Results

A total of 1,747,604 CNVs were detected in this study, and 7228 CNV regions (CNVR) were obtained after merging overlapping CNVs; these regions accounted for 2.17% of the sheep reference genome. The average length of the CNVRs was 4307.17 bp. “Deletion” events took place more frequently than “duplication” or “both” events. The CNVRs obtained overlapped with previously reported sheep CNVRs to variable extents (4.39–55.46%). Functional enrichment analysis showed that the CNVR-harboring genes were mainly involved in sensory perception systems, nutrient metabolism processes, and growth and development processes. Furthermore, 1855 of the CNVRs were associated with 166 quantitative trait loci (QTL), including milk QTLs, carcass QTLs, and health-related QTLs, among others. In addition, the 32 fine-wool sheep were divided into horned and polled groups to analyze for the selective sweep of CNVRs, and it was found that the relaxin family peptide receptor 2 (RXFP2) gene was strongly influenced by selection.

Conclusions

In summary, we constructed a genomic CNV map for Chinese indigenous fine-wool sheep using resequencing, thereby providing a valuable genetic variation resource for sheep genome research, which will contribute to the study of complex traits in sheep.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07387-7.

Genomic Prediction of Additive and Dominant Effects on Wool and Blood Traits in Alpine Merino Sheep

Dominant genetic effects may provide a critical contribution to the total genetic variation of quantitative and complex traits. However, investigations of genome-wide markers to study the genomic prediction (GP) and genetic mechanisms of complex traits generally ignore dominant genetic effects. The increasing availability of genomic datasets and the potential benefits of the inclusion of non-additive genetic effects in GP have recently renewed attention to incorporation of these effects in genomic prediction models. In the present study, data from 498 genotyped Alpine Merino sheep were adopted to estimate the additive and dominant genetic effects of 9 wool and blood traits via two linear models: (1) an additive effect model (MAG) and (2) a model that included both additive and dominant genetic effects (MADG). Moreover, a method of 5-fold cross validation was used to evaluate the capability of GP in the two different models. The results of variance component estimates for each trait suggested that for fleece extension rate (73%), red blood cell count (28%), and hematocrit (25%), a large component of phenotypic variation was explained by dominant genetic effects. The results of cross validation demonstrated that the MADG model, comprising additive and dominant genetic effects, did not display an apparent advantage over the MAG model that included only additive genetic effects, i.e., the model that included dominant genetic effects did not improve the capability for prediction of the genomic model. Consequently, inclusion of dominant effects in the GP model may not be beneficial for wool and blood traits in the population of Alpine Merino sheep.

Genetic signatures of high-altitude adaptation and geographic distribution in Tibetan sheep

Most sheep breeding programs designed for the tropics and sub-tropics have to take into account the impacts of environmental adaptive traits. However, the genetic mechanism regulating the multiple biological processes driving adaptive responses remains unclear. In this study, we applied a selective sweep analysis by combing 1% top values of F_st and ZHp on both altitude and geographic subpopulations (APS) in 636 indigenous Tibetan sheep breeds. Results show that 37 genes were identified within overlapped genomic regions regarding F_st significantly associated with APS. Out of the 37 genes, we found that 8, 3 and 6 genes at chromosomes (chr.) 13, 23 and 27, respectively, were identified in the genomic regions with 1% top values of ZHp. We further analyzed the INDEL variation of 6 genes at chr.27 (X chromosome) in APS together with corresponding orthologs of 6 genes in Capra, Pantholops, and Bos Taurus. We found that an INDEL was located within 5′UTR region of HAG1 gene. This INDEL of HAG1 was strongly associated with the variation of APS, which was further confirmed by qPCR. Sheep breeds carrying “C-INDEL” of HAG1 have significantly greater body weight, shear amount, corpuscular hemoglobin and globulin levels, but lower body height, than those carrying “CA-INDEL” of HAG1. We concluded that “C-INDEL” variation of HAG1 gene confers better hypoxia tolerance in the highlands of Tibetan and explains well geographic distributions in this population. These results contribute to our understanding of adaptive responses to altitude and geographic adaptation in Tibetan sheep populations and will help to guide future conservation programs for Tibetan sheep native to Qinghai-Tibetan Plateau.

Analysis of dynamic and widespread lncRNA and miRNA expression in fetal sheep skeletal muscle

The sheep is an economically important animal, and there is currently a major focus on improving its meat quality through breeding. There are variations in the growth regulation mechanisms of different sheep breeds, making fundamental research on skeletal muscle growth essential in understanding the regulation of (thus far) unknown genes. Skeletal muscle development is a complex biological process regulated by numerous genes and non-coding RNAs, including microRNAs (miRNAs) and long non-coding RNAs (lncRNAs). In this study, we used deep sequencing data from sheep longissimus dorsi (LD) muscles sampled at day 60, 90, and 120 of gestation, as well as at day 0 and 360 following birth, to identify and examine the lncRNA and miRNA temporal expression profiles that regulate sheep skeletal myogenesis. We stained LD muscles using histological sections to analyse the area and circumference of muscle fibers from the embryonic to postnatal development stages. Our results showed that embryonic skeletal muscle growth can be characterized by time. We obtained a total of 694 different lncRNAs and compared the differential expression between the E60 vs. E90, E90 vs. E120, E120 vs. D0, and D0 vs. D360 lncRNA and gene samples. Of the total 701 known sheep miRNAs we detected, the following showed a wide range of expression during the embryonic stage: miR-2387, miR-105, miR-767, miR-432, and miR-433. We propose that the detected lncRNA expression was time-specific during the gestational and postnatal stages. GO and KEGG analyses of the genes targeted by different miRNAs and lncRNAs revealed that these significantly enriched processes and pathways were consistent with skeletal muscle development over time across all sampled stages. We found four visual lncRNA–gene regulatory networks that can be used to explore the function of lncRNAs in sheep and may be valuable in helping improve muscle growth. This study also describes the function of several lncRNAs that interact with miRNAs to regulate myogenic differentiation.

Enhanced transport of spin-orbit-coupled Bose gases in disordered potentials

Anderson localization is a single-particle localization phenomena in disordered media that is accompanied by an absence of diffusion. Spin-orbit coupling (SOC) describes an interaction between a particle's spin and its momentum that directly affects its energy dispersion, for example, creating dispersion relations with gaps and multiple local minima. We show theoretically that combining one-dimensional spin-orbit coupling with a transverse Zeeman field suppresses the effects of disorder, thereby increasing the localization length and conductivity. This increase results from a suppression of backscattering between states in the gap of the SOC dispersion relation. Here, we focus specifically on the interplay of disorder from an optical speckle potential and SOC generated by two-photon Raman processes in quasi-one-dimensional Bose-Einstein condensates. We first describe backscattering by using a Fermi golden rule approach, and then numerically confirm this picture by solving the time-dependent one-dimensional Gross-Pitaevskii equation for a weakly interacting Bose-Einstein condensate with SOC and disorder. We find that on the tens of millisecond timescale of typical cold atom experiments moving in harmonic traps, initial states with momentum in the zero-momentum SOC gap evolve with negligible backscattering, while without SOC these same states rapidly localize.

Genome-Wide Association Study Using Individual Single-Nucleotide Polymorphisms and Haplotypes for Erythrocyte Traits in Alpine Merino Sheep

Adaptation to high-altitude hypoxia is essential for domestic animals, such as yak, Tibetan chicken, and Tibetan sheep, living on high plateaus, as it ensures efficient oxygen absorption and utilization. Red blood cells are the primary medium for transporting oxygen in the blood. However, little is known about the genetic mechanism of erythrocyte traits. Genome-wide association studies (GWASs) based on single markers or haplotypes have identified potential mechanisms for genetic variation and quantitative traits. To identify loci associated with erythrocyte traits, we performed a GWAS based on the method of the single marker and haplotype in 498 Alpine Merino sheep for six erythrocyte traits: red blood cell count (RBC), hemoglobin (HGB), hematocrit (HCT), mean corpuscular hemoglobin (MCH), mean corpuscular hemoglobin concentration (MCHC), and RBC volume distribution width coefficient of variation (RWD_CV). Forty-two significant single-nucleotide polymorphisms (SNPs) associated with the six erythrocyte traits were detected by means of a single-marker GWAS, and 34 significant haplotypes associated with five erythrocyte traits were detected by means of haplotype analysis. We identified six genes (DHCR24, SPATA9, FLI1, PLCB1, EFNB2, and SH2B3) as potential genes of interest via gene function annotations, location, and expression variation. In particular, FLI1 and PLCB1 were associated with hematopoiesis and erythropoiesis, respectively. These results provide a theoretical basis for analyzing erythrocyte traits and high-altitude hypoxia adaptation in Alpine Merino sheep and will be a useful reference for future studies of plateau-dwelling livestock.

Activation of the STAT1 Pathway Accelerates Periodontitis in Nos3-/- Mice

Early studies on the etiology and pathogenesis of hypertension have shown that it has a considerable association with inflammation and the immune response as well as periodontitis. Clinical studies have also shown that hypertension can promote the periodontal tissue destruction caused by periodontitis. However, the underlying mechanisms remain unclear. This study aimed to explore the possible mechanisms of how hypertension aggravates periodontitis. Treatment with or without the signal transducer and activator of transcription 1 (STAT1) inhibitor fludarabine was performed in an endothelial nitric oxide synthase gene knockout-related (Nos3^-/-) mouse model with the hypertension phenotype of periodontitis induced by bacteria. Micro-computed tomography, immunohistochemistry, Western blot, quantitative reverse transcription polymerase chain reaction, immunofluorescence, and ELISA were performed. We demonstrated that Nos3^-/--related hypertension increases bone resorption and periodontal destruction in periodontitis lesion areas, which can be inhibited by the STAT1 inhibitor. Experimental data also showed that Nos3^-/- significantly increased macrophage infiltration and proinflammatory cytokine expression in the periodontitis lesion area, which is dependent on the angiotensin II-induced STAT1 pathway. Inhibition of STAT1 in vivo can decrease the expression of proinflammatory cytokines and macrophage infiltration. Furthermore, data in this study showed that Nos3^-/--related hypertension further downregulated the STAT3 anti-inflammatory function and its downstream chemokine expression in a STAT1-dependent manner. By applying RAW 264.7 and L929 cell lines and monocytes isolated from Nos3^-/- mice, we confirmed that activation of the STAT1 pathway inhibits STAT3 and its downstream pathway and promotes inflammatory cytokine expression in vitro. Collectively, our current study demonstrated that STAT1 plays an indispensable role in the Nos3^-/--related hypertension with aggravation of periodontitis, suggesting that STAT1 may be a key target for the treatment of periodontitis with hypertension.

Follicle stimulating hormone and estradiol trajectories from menopausal transition to late postmenopause in indigenous Chinese women

OBJECTIVE

The aim of this study was to examine follicle stimulating hormone (FSH) and estradiol (E2) trajectories in indigenous Chinese women through the ovarian aging process from 10 years before the final menstrual period (FMP) to 20 years after.

METHODS

A prospective community-based longitudinal cohort study of 362 Chinese middle-aged women with a clearly determined FMP was performed.

RESULTS

FSH and E2 trajectories were identified according to years from FMP and the Stages of Reproductive Aging Workshop + 10 (STRAW + 10), and further classified by body mass index. E2 decreases relatively steadily from Stage -2 to +1c, while FSH stabilizes quickly from Stage -2 to +1a. Obese women have higher E2 levels and lower FSH levels during menopausal transition, which last for decades after the FMP. There was no obvious decline in FSH levels for at least 20 years after the FMP.

CONCLUSIONS

The FSH and E2 trajectories in this indigenous Chinese cohort appear similar to those discussed in the Study of Women's Health Across the Nation, with ethnic differences. STRAW + 10 criteria may be used to ascertain the ovarian aging process in Chinese women. Obese women may experience different levels of hormonal fluctuations during the menopause transition, while FSH levels in all women appear to remain high even at late postmenopause.

Comparative proteomics reveals genetic mechanisms underlying secondary hair follicle development in fine wool sheep during the fetal stage

The aim of this study was to investigate the genetic mechanisms underlying wool production by characterizing the skin protein profile and determining the proteomic changes that occur as a consequence of development in wool-producing sheep using a label-free proteomics approach. Samples were collected at four stages during gestation (87, 96, 102, and 138 days), and every two consecutive stages were statistically compared (87 versus 96, 96 versus 102, and 102 versus 138 days). We identified 227 specific proteins in the sheep proteome that were present in all four stages, and 123 differentially abundant proteins (DAPs). We also observed that the microstructure of the secondary follicles changed significantly during the development of the fetal skin hair follicle. The screened DAPs were strictly related to metabolic and skin development pathways, and were associated with pathways such as the glycolysis/gluconeogenesis. These analyses indicated that the wool production of fine wool sheep is regulated via a variety of pathways. These findings provide an important resource that can be used in future studies of the genetic mechanisms underlying wool traits in fine wool sheep, and the identified DAPs should be further investigated as candidate markers for predicting wool traits in sheep. SIGNIFICANCE: Wool quality (fiber diameter, length, etc.) is an important economic trait of fine wool sheep that is determined by secondary follicle differentiation and re-differentiation. Secondary follicles of fine wool sheep developed from a bud (87 days), and underwent differentiation (96 days) and rapid growth (102 days) until maturity (138 days) during gestation. Comparative analysis based on differential proteomics of these four periods could provide a better understanding of the wool growth mechanism of fine wool sheep and offer novel strategies for improving fine wool quality by breeding.

The complete mitochondrial genome sequence and phylogenetic analysis of Alpine Merino sheep (Ovis aries)

Alpine Merino sheep is one of the most important fine-wool sheep breeds in China. In this study, we present the complete mitogenome of Alpine Merino sheep for the first time. The genome has a length of 16,619bp, containing 13 protein-coding genes, 2 ribosomal RNA genes, 22 transfer RNA genes, and a control region (D-loop). Phylogenetically, the Alpine Merino sheep is closer to Oula Tibetan sheep and Tashkurgan sheep. This report provided new data for the phylogeny of Alpine Merino sheep.

Genome-Wide Association Study of Body Weight Traits in Chinese Fine-Wool Sheep

Body weight is an important economic trait for sheep and it is vital for their successful production and breeding. Therefore, identifying the genomic regions and biological pathways that contribute to understanding variability in body weight traits is significant for selection purposes. In this study, the genome-wide associations of birth, weaning, yearling, and adult weights of 460 fine-wool sheep were determined using resequencing technology. The results showed that 113 single nucleotide polymorphisms (SNPs) reached the genome-wide significance levels for the four body weight traits and 30 genes were annotated effectively, including AADACL3, VGF, NPC1, and SERPINA12. The genes annotated by these SNPs significantly enriched 78 gene ontology terms and 25 signaling pathways, and were found to mainly participate in skeletal muscle development and lipid metabolism. These genes can be used as candidate genes for body weight in sheep, and provide useful information for the production and genomic selection of Chinese fine-wool sheep.

Long non-coding RNA expression profile in permanent atrial fibrillation patients with rheumatic heart disease

OBJECTIVE

Atrial fibrillation (AF) is the most common type of arrhythmia, especially in rheumatic heart disease (RHD) patients. The differences in structural remodeling and electrical remodeling between the left and right atrium associated with AF in RHD patients are well known, and alterations in the expression profiles of long noncoding RNAs (lncRNAs) in the left atrium have also been investigated. However, the role of lncRNAs in the right atrium (RA) remains largely unknown.

PATIENTS AND METHODS

We identified differentially expressed lncRNAs in RA tissues of RHD patients with AF or a normal sinus rhythm (NSR) using microarray analysis. Then, we performed gene ontology (GO) and KEGG pathway analyses for functional annotation of the deregulated lncRNAs. Finally, we constructed a lncRNA-mRNA co-expression network.

RESULTS

Of the 22,829 human non-coding RNAs analyzed, a total of 1,909 long non-coding RNAs were detected. A total of 182 lncRNAs (117 downregulated and 65 upregulated) were shown to be differentially expressed (fold-change > 1.5) in AF patients compared with NSR patients. Many lncRNAs might be partially involved in an AF-related pathway.

CONCLUSIONS

AF dysregulates the expression of lncRNAs in the RA of RHD patients. These findings may be useful for exploring potential therapeutic treatments for AF in RHD patients.

Two linked TBXT (brachyury) gene polymorphisms are associated with the tailless phenotype in fat-rumped sheep

T‐box transcription factor T (TBXT), encoding the brachyury protein, is an embryonic nuclear transcription factor involved in mesoderm formation and differentiation. Previous studies indicate that TBXT mutations are responsible for the tailless or short‐tailed phenotype of many vertebrates. To verify whether the tailless phenotype in fat‐rumped sheep is associated with TBXT mutations, exon 2 of the TBXT gene for 301 individuals belonging to 13 Chinese and Iranian sheep breeds was directly sequenced. Meanwhile, 380 samples were used to detect the genotypes of the candidate variations by mapping to their reads databases in the Sequence Read Archive repository of GenBank. The results showed that one missense mutation, c.334G>T (GGG>TGG) with a completely linked synonymous variant c.333G>C (CCG>CCC) was found to be associated with the ‘tailless’ characteristic in typical fat‐rumped sheep breeds. The c.334G>T transversion led to the conversion of glycine to tryptophan at the 112th amino acid in the T‐box domain of the brachyury protein. In addition, crossbreeding experiments for long‐tailed and tailless sheep showed that CT/CT allele of nucleotides (nt) 333 and 334, a recessive mutation, would cause sheep tails to be shorter, suggesting that these two linked variants at nucleotides 333 and 334 in TBXT are probably causative mutations responsible for the tailless phenotype in sheep.

[The comparison of preemptive analgesic efficacy between short-acting ketamine and long-acting parecoxib]

Objective: To compare the analgesic efficacy between preoperative single-dose ketamine, a short-acting medicine and parecoxib, a long-acting medicine for reducing analgesic consumption in the first 24 h after-operation. Methods: Eighty-one patients from Beijing Chaoyang Hospital undergoing laparoscopic uterus surgery between April and December 2015 were randomly divided into three groups: control group (group C), ketamine group (group K) and parecoxib sodium group (group P). All patients were anesthetized with general anesthesia and received sufentanil-based patient-controlled intravenous analgesia (PCIA). After induction and 10 min before incision, patients in group K, P and C were injected with intravenous 0.5 mg/kg ketamine, 40 mg parecoxib and 2 ml normal saline, respectively. The primary outcome was sufentanil consumption within 1 h and 24 h after surgery. Other outcomes included the visual analog scale (VAS) pain score at 0, 15, 30, 45min and 1, 2, 4, 8, 24 h after surgery, PCIA effective trigger times, and adverse reactions. Results: The postoperative sufentanil consumptions within 1 h in group K and P were(4.420±1.836)μg and (2.878±1.984)μg, respectively, and consumptions within 24 h were(28.200±3.712)μg and (25.511±4.037)μg, respectively, which were significantly less than that in group C with (6.144±2.346)μg within 1 h and (31.505±7.042)μg within 24 h (F=15.360, 8.406, all P<0.05). Patients in group P needed less sufentanil than group K in 1 h after surgery (P<0.05), however, the difference was not statistically significant in 24 h(P>0.05). The PCIA trigger times were 2(3.75) in group C, 0(1.50) in group K, and 0(1.00) in group P. Group K and P had less PCIA trigger times compared to group C (all P<0.05). Compared to group C, group K and P had lower VAS scores at 0, 15, 30 min after surgery and group P had lower VAS scores at 1, 15, 30, 45 min, 2 h after surgery, respectively (all P<0.05). There were no differences between groups in the incidence of any adverse effects(all P>0.05). Conclusion: A single injection of short-acting ketamine before laparoscopic uterus surgery, has the same efficacy as long-acting parecoxib for opioid-sparing effect in the first 24 h after-operation. However, parecoxib has better analgesic effect in the early postoperative period.

Accuracy of interferon-γ-induced protein 10 for diagnosing latent tuberculosis infection: a systematic review and meta-analysis

BACKGROUND

Effective diagnostic methods for detecting latent tuberculosis infection (LTBI) are important for its eradication. A number of studies have evaluated the use of interferon-γ-induced protein 10 (IP-10), which is elevated after tuberculosis infection, as a biomarker for LTBI, but conclusive results regarding its effectiveness have not been reported.

OBJECTIVES

Our objective was to assess the diagnostic value of IP-10 for LTBI.

DATA SOURCES

We searched the PubMed, Embase, the Cochrane Library and Web of Science databases to find eligible studies.

STUDY ELIGIBILITY CRITERIA

We included cohort, case-control and cross-sectional studies that evaluated IP-10 in LTBI participants in comparison with tuberculin skin tests (TST) and interferon-γ release assays (IGRA).

PARTICIPANTS

Individuals with LTBI and uninfected participants.

INTERVENTIONS

IP-10 (index test) compared with TST and IGRA (reference standard) for diagnosing LTBI.

METHODS

PubMed, Embase, the Cochrane Library, and Web of Science databases were searched up to June 2018. A hierarchical summary receiver operating characteristic (HSROC) model was used to evaluate the pooled sensitivity, specificity, positive likelihood ratio (PLR), negative likelihood ratio (NLR), diagnostic odds ratio (DOR), and HSROC curve for the diagnostic efficiency of IP-10.

RESULTS

Twelve studies including 1023 participants and 1122 samples were included. The overall pooled sensitivity was 0.85 (95% CI 0.80-0.88), specificity was 0.89 (95% CI 0.84-0.92), PLR was 7.55 (95% CI 5.20-10.97), NLR was 0.17 (95% CI 0.13-0.22) and DOR was 44.23 (95% CI 28.86-67.79), indicating a high accuracy for diagnosing LTBI. Based on a meta-regression analysis, high-burden countries, study design, IP-10 method, reference standard and the IP-10 cut-off could not explain the heterogeneity (p >0.05).

CONCLUSIONS

Our results suggested that IP-10 is a promising biomarker for the diagnosis of LTBI.