Late-life dietary folate restriction reduces biosynthetic processes without compromising healthspan in mice.

Folate is a vitamin required for cell growth and is present in fortified foods in the form of folic acid to prevent congenital abnormalities. The impact of low folate status on life-long health is poorly understood. We found that limiting folate levels with the folate antagonist methotrexate increased the lifespan of yeast and worms. We then restricted folate intake in aged mice and measured various health metrics, metabolites, and gene expression signatures. Limiting folate intake decreased anabolic biosynthetic processes in mice and enhanced metabolic plasticity. Despite reduced serum folate levels in mice with limited folic acid intake, these animals maintained their weight and adiposity late in life, and we did not observe adverse health outcomes. These results argue that the effectiveness of folate dietary interventions may vary depending on an individual's age and sex. A higher folate intake is advantageous during the early stages of life to support cell divisions needed for proper development. However, a lower folate intake later in life may result in healthier aging.

Genetic mapping and prediction for novel lesion mimic in maize demonstrates quantitative effects from genetic background, environment and epistasis.

KEY MESSAGE: A novel locus was discovered on chromosome 7 associated with a lesion mimic in maize; this lesion mimic had a quantitative and heritable phenotype and was predicted better via subset genomic markers than whole genome markers across diverse environments. Lesion mimics are a phenotype of leaf micro-spotting in maize (Zea mays L.), which can be early signs of biotic or abiotic stresses. Dissecting its inheritance is helpful to understand how these loci behave across different genetic backgrounds. Here, 538 maize recombinant inbred lines (RILs) segregating for a novel lesion mimic were quantitatively phenotyped in Georgia, Texas, and Wisconsin. These RILs were derived from three bi-parental crosses using a tropical pollinator (Tx773) as the common parent crossed with three inbreds (LH195, LH82, and PB80). While this lesion mimic was heritable across three environments based on phenotypic ([Formula: see text] = 0.68) and genomic ([Formula: see text] = 0.91) data, transgressive segregation was observed. A genome-wide association study identified a single novel locus on chromosome 7 (at 70.6 Mb) also covered by a quantitative trait locus interval (69.3-71.0 Mb), explaining 11-15% of the variation, depending on the environment. One candidate gene identified in this region, Zm00001eb308070, is related to the abscisic acid pathway involving in cell death. Genomic predictions were applied to genome-wide markers (39,611 markers) contrasted with a marker subset (51 markers). Population structure explained more variation than environment in genomic prediction, but other substantial genetic background effects were additionally detected. Subset markers explained substantially less genetic variation (24.9%) for the lesion mimic than whole genome markers (55.4%) in the model, yet predicted the lesion mimic better (0.56-0.66 vs. 0.26-0.29). These results indicate this lesion mimic phenotype was less affected by environment than by epistasis and genetic background effects, which explain its transgressive segregation.

Horn fly transcriptome data of ten populations from the southern United States with varying degrees and molecular mechanisms of pesticide resistance.

Haematobia irritans irritans (Linnaeus, 1758: Diptera: Muscidae), the horn fly, is an external parasite of penned and pastured livestock that causes a major economic impact on cattle production worldwide. Pesticides such as synthetic pyrethroids and organophosphates are routinely used to control horn flies; however, resistance to these chemicals has become a concern in several countries. To further elucidate the molecular mechanisms of resistance in horn fly populations, we sequenced the transcriptomes of ten populations of horn flies from the southern US possessing varying degrees of pesticide resistance levels to pyrethroids, organophosphates, and endosulfans. We employed an Illumina paired end HiSeq approach, followed by de novo assembly of the transcriptomes using CLC Genomics Workbench 8.0.1 De Novo Assembler using multiple kmers, and annotation using Blast2GO PRO version 5.2.5. The Gene Ontology biological process term Response to Insecticide was found in all the populations, but at an increased frequency in the populations with higher levels of insecticide resistance. The raw sequence reads are archived in the Sequence Read Archive (SRA) and assembled population transcriptomes in the Transcriptome Shotgun Assembly (TSA) at the National Center for Biotechnology Information (NCBI).

Comparative transcriptomic analysis of germinating rice seedlings to individual and combined anaerobic and cold stress.

BACKGROUND: Rice is one of the most important cereals consumed worldwide. Two major abiotic factors affecting rice plants in different growth stages are flooding stress and cold stress. These abiotic stresses can take place independently or simultaneously and significantly affect rice plants during germination and seedling growth. Fortunately, a wide array of phenotypic responses conferring flooding stress and chilling stress tolerance exist within the rice germplasm, indicating the presence of different molecular mechanisms underlying tolerance to these stresses. Understanding these differences may assist in developing improved rice cultivars having higher tolerance to both stresses. In this study, we conducted a comparative global gene expression analysis of two rice genotypes with contrasting phenotypes under cold stress, anaerobic stress, and combined cold and anaerobic stress during germination. RESULTS: The differential gene expression analysis revealed that 5571 differentially expressed genes (DEGs), 7206 DEGs, and 13279 DEGs were identified under anaerobic stress, cold stress, and combined stress, respectively. Genes involved in the carbohydrate metabolic process, glucosyltransferase activity, regulation of nitrogen compound metabolic process, protein metabolic process, lipid metabolic process, cellular nitrogen compound biosynthetic process, lipid biosynthetic process, and a microtubule-based process were enriched across all stresses. Notably, the common Gene Ontology (GO) analysis identified three hub genes, namely Os08g0176800 (similar to mRNA-associated protein mrnp 41), Os11g0454200 (dehydrin), and OS10g0505900 (expressed protein). CONCLUSION: A large number of differentially expressed genes were identified under anaerobic, cold conditions during germination and the combination of the two stress conditions in rice. These results will assist in the identification of promising candidate genes for possible manipulation toward rice crops that are more tolerant under flooding and cold during germination, both independently and concurrently.

Transcriptome Analysis of Developing Grains from Wheat Cultivars TAM 111 and TAM 112 Reveal Cultivar-Specific Regulatory Networks.

Wheat flour's end-use quality is tightly linked to the quantity and composition of storage proteins in the endosperm. TAM 111 and TAM 112 are two popular cultivars grown in the Southern US Great Plains with significantly different protein content. To investigate regulatory differences, transcriptome data were analyzed from developing grains at early- and mid-filling stages. At the mid-filling stage, TAM 111 preferentially upregulated starch metabolism-related pathways compared to TAM 112, whereas amino acid metabolism and transporter-related pathways were over-represented in TAM 112. Elemental analyses also indicated a higher N percentage in TAM 112 at the mid-filling stage. To explore the regulatory variation, weighted correlation gene network was constructed from publicly available RNAseq datasets to identify the modules differentially regulated in TAM 111 and TAM 112. Further, the potential transcription factors (TFs) regulating those modules were identified using graphical least absolute shrinkage and selection operator (GLASSO). Homologs of the OsNF-Y family members with known starch metabolism-related functions showed higher connectivities in TAM 111. Multiple TFs with high connectivity in TAM 112 had predicted functions associated with ABA response in grain. These results will provide novel targets for breeders to explore and further our understanding in mechanisms regulating grain development.

Characterization of triatomine bloodmeal sources using direct Sanger sequencing and amplicon deep sequencing methods.

Knowledge of host associations of blood-feeding vectors may afford insights into managing disease systems and protecting public health. However, the ability of methods to distinguish bloodmeal sources varies widely. We used two methods-Sanger sequencing and amplicon deep sequencing-to target a 228 bp region of the vertebrate Cytochrome b gene and determine hosts fed upon by triatomines (n = 115) collected primarily in Texas, USA. Direct Sanger sequencing of PCR amplicons was successful for 36 samples (31%). Sanger sequencing revealed 15 distinct host species, which included humans, domestic animals (Canis lupus familiaris, Ovis aries, Gallus gallus, Bos taurus, Felis catus, and Capra hircus), wildlife (Rattus rattus, Incilius nebulifer, Sciurus carolinensis, Sciurus niger, and Odocoileus virginianus), and captive animals (Panthera tigris, Colobus spp., and Chelonoidis carbonaria). Samples sequenced by the Sanger method were also subjected to Illumina MiSeq amplicon deep sequencing. The amplicon deep sequencing results (average of 302,080 usable reads per sample) replicated the host community revealed using Sanger sequencing, and detected additional hosts in five triatomines (13.9%), including two additional blood sources (Procyon lotor and Bassariscus astutus). Up to four bloodmeal sources were detected in a single triatomine (I. nebulifer, Homo sapiens, C. lupus familiaris, and S. carolinensis). Enhanced understanding of vector-host-parasite networks may allow for integrated vector management programs focusing on highly-utilized and highly-infected host species.

Brucella activates the host RIDD pathway to subvert BLOS1-directed immune defense.

The phagocytosis and destruction of pathogens in lysosomes constitute central elements of innate immune defense. Here, we show that Brucella, the causative agent of brucellosis, the most prevalent bacterial zoonosis globally, subverts this immune defense pathway by activating regulated IRE1α-dependent decay (RIDD) of Bloc1s1 mRNA encoding BLOS1, a protein that promotes endosome-lysosome fusion. RIDD-deficient cells and mice harboring a RIDD-incompetent variant of IRE1α were resistant to infection. Inactivation of the Bloc1s1 gene impaired the ability to assemble BLOC-1-related complex (BORC), resulting in differential recruitment of BORC-related lysosome trafficking components, perinuclear trafficking of Brucella-containing vacuoles (BCVs), and enhanced susceptibility to infection. The RIDD-resistant Bloc1s1 variant maintains the integrity of BORC and a higher-level association of BORC-related components that promote centrifugal lysosome trafficking, resulting in enhanced BCV peripheral trafficking and lysosomal destruction, and resistance to infection. These findings demonstrate that host RIDD activity on BLOS1 regulates Brucella intracellular parasitism by disrupting BORC-directed lysosomal trafficking. Notably, coronavirus murine hepatitis virus also subverted the RIDD-BLOS1 axis to promote intracellular replication. Our work establishes BLOS1 as a novel immune defense factor whose activity is hijacked by diverse pathogens.

A new strategy for using historical imbalanced yield data to conduct genome-wide association studies and develop genomic prediction models for wheat breeding.

Using imbalanced historical yield data to predict performance and select new lines is an arduous breeding task. Genome-wide association studies (GWAS) and high throughput genotyping based on sequencing techniques can increase prediction accuracy. An association mapping panel of 227 Texas elite (TXE) wheat breeding lines was used for GWAS and a training population to develop prediction models for grain yield selection. An imbalanced set of yield data collected from 102 environments (year-by-location) over 10 years, through testing yield in 40-66 lines each year at 6-14 locations with 38-41 lines repeated in the test in any two consecutive years, was used. Based on correlations among data from different environments within two adjacent years and heritability estimated in each environment, yield data from 87 environments were selected and assigned to two correlation-based groups. The yield best linear unbiased estimation (BLUE) from each group, along with reaction to greenbug and Hessian fly in each line, was used for GWAS to reveal genomic regions associated with yield and insect resistance. A total of 74 genomic regions were associated with grain yield and two of them were commonly detected in both correlation-based groups. Greenbug resistance in TXE lines was mainly controlled by Gb3 on chromosome 7DL in addition to two novel regions on 3DL and 6DS, and Hessian fly resistance was conferred by the region on 1AS. Genomic prediction models developed in two correlation-based groups were validated using a set of 105 new advanced breeding lines and the model from correlation-based group G2 was more reliable for prediction. This research not only identified genomic regions associated with yield and insect resistance but also established the method of using historical imbalanced breeding data to develop a genomic prediction model for crop improvement. Supplementary Information: The online version contains supplementary material available at 10.1007/s11032-022-01287-8.

Reference-Guided De Novo Genome Assembly to Dissect a QTL Region for Submergence Tolerance Derived from Ciherang-Sub1.

Global climate change has increased the number of severe flooding events that affect agriculture, including rice production in the U.S. and internationally. Heavy rainfall can cause rice plants to be completely submerged, which can significantly affect grain yield or completely destroy the plants. Recently, a major effect submergence tolerance QTL during the vegetative stage, qSub8.1, which originated from Ciherang-Sub1, was identified in a mapping population derived from a cross between Ciherang-Sub1 and IR10F365. Ciherang-Sub1 was, in turn, derived from a cross between Ciherang and IR64-Sub1. Here, we characterize the qSub8.1 region by analyzing the sequence information of Ciherang-Sub1 and its two parents (Ciherang and IR64-Sub1) and compare the whole genome profile of these varieties with the Nipponbare and Minghui 63 (MH63) reference genomes. The three rice varieties were sequenced with 150 bp pair-end whole-genome shotgun sequencing (Illumina HiSeq4000), followed by performing the Trimmomatic-SOAPdenovo2-MUMmer3 pipeline for genome assembly, resulting in approximate genome sizes of 354.4, 343.7, and 344.7 Mb, with N50 values of 25.1, 25.4, and 26.1 kb, respectively. The results showed that the Ciherang-Sub1 genome is composed of 59-63% Ciherang, 22-24% of IR64-Sub1, and 15-17% of unknown sources. The genome profile revealed a more detailed genomic composition than previous marker-assisted breeding and showed that the qSub8.1 region is mostly from Ciherang, with some introgressed segments from IR64-Sub1 and currently unknown source(s).

Genome-wide QTL mapping of yield and agronomic traits in two widely adapted winter wheat cultivars from multiple mega-environments.

Quantitative trait loci (QTL) analysis could help to identify suitable molecular markers for marker-assisted breeding (MAB). A mapping population of 124 F5:7recombinant inbred lines derived from the cross 'TAM 112'/'TAM 111' was grown under 28 diverse environments and evaluated for grain yield, test weight, heading date, and plant height. The objective of this study was to detect QTL conferring grain yield and agronomic traits from multiple mega-environments. Through a linkage map with 5,948 single nucleotide polymorphisms (SNPs), 51 QTL were consistently identified in two or more environments or analyses. Ten QTL linked to two or more traits were also identified on chromosomes 1A, 1D, 4B, 4D, 6A, 7B, and 7D. Those QTL explained up to 13.3% of additive phenotypic variations with the additive logarithm of odds (LOD(A)) scores up to 11.2. The additive effect increased yield up to 8.16 and 6.57 g m-2 and increased test weight by 2.14 and 3.47 kg m-3 with favorable alleles from TAM 111 and TAM 112, respectively. Seven major QTL for yield and six for TW with one in common were of our interest on MAB as they explained 5% or more phenotypic variations through additive effects. This study confirmed previously identified loci and identified new QTL and the favorable alleles for improving grain yield and agronomic traits.

Case Report: Paucisymptomatic College-Age Population as a Reservoir for Potentially Neutralization-Resistant Severe Acute Respiratory Syndrome Coronavirus 2 Variants.

To better understand the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variant lineage distribution in a college campus population, we carried out viral genome surveillance over a 7-week period from January to March 2021. Among the sequences were three novel viral variants: BV-1 with a B.1.1.7/20I genetic background and an additional spike mutation Q493R, associated with a mild but longer-than-usual COVID-19 case in a college-age person, BV-2 with a T478K mutation on a 20B genetic background, and BV-3, an apparent recombinant lineage. This work highlights the potential of an undervaccinated younger population as a reservoir for the spread and generation of novel variants. This also demonstrates the value of whole genome sequencing as a routine disease surveillance tool.

Raw pacific biosciences and illumina sequencing reads and assembled genome data for the cattle ticks Rhipicephalus microplus and Rhipicephalus annulatus.

Ticks from the genus Rhipicephalus have enormous global economic impact as ectoparasites of cattle. Rhipicephalus microplus and Rhipicephalus annulatus are known to harbor infectious pathogens such as Babesia bovis, Babesia bigemina, and Anaplasma marginale. Having reference quality genomes of these ticks would advance research to identify druggable targets for chemical entities with acaricidal activity and refine anti-tick vaccine approaches. We sequenced and assembled the genomes of R. microplus and R. annulatus, using Pacific Biosciences and HiSeq 4000 technologies on very high molecular weight genomic DNA. We used 22 and 29 SMRT cells on the Pacific Biosciences Sequel for R. microplus and R. annulatus, respectively, and 3 lanes of the Illumina HiSeq 4000 platform for each tick. The PacBio sequence yields for R. microplus and R. annulatus were 21.0 and 27.9 million subreads, respectively, which were assembled with Canu v. 1.7. The final Canu assemblies consisted of 92,167 and 57,796 contigs with an average contig length of 39,249 and 69,055 bp for R. microplus and R. annulatus, respectively. Annotated genome quality was assessed by BUSCO analysis to provide quantitative measures for each assembled genome. Over 82% and 92% of the 1066 member BUSCO gene set was found in the assembled genomes of R. microplus and R. annulatus, respectively. For R. microplus, only 189 of the 1066 BUSCO genes were missing and only 140 were present in a fragmented condition. For R. annulatus, only 75 of the BUSCO genes were missing and only 109 were present in a fragmented condition. The raw sequencing reads and the assembled contigs/scaffolds are archived at the National Center for Biotechnology Information.

RNA-seq analysis reveals different drought tolerance mechanisms in two broadly adapted wheat cultivars 'TAM 111' and 'TAM 112'.

Wheat cultivars 'TAM 111' and 'TAM 112' have been dominantly grown in the Southern U.S. Great Plains for many years due to their high yield and drought tolerance. To identify the molecular basis and genetic control of drought tolerance in these two landmark cultivars, RNA-seq analysis was conducted to compare gene expression difference in flag leaves under fully irrigated (wet) and water deficient (dry) conditions. A total of 2254 genes showed significantly altered expression patterns under dry and wet conditions in the two cultivars. TAM 111 had 593 and 1532 dry-wet differentially expressed genes (DEGs), and TAM 112 had 777 and 1670 at heading and grain-filling stages, respectively. The two cultivars have 1214 (53.9%) dry-wet DEGs in common, which agreed with their excellent adaption to drought, but 438 and 602 dry-wet DEGs were respectively shown only in TAM 111 and TAM 112 suggested that each has a specific mechanism to cope with drought. Annotations of all 2254 genes showed 1855 have functions related to biosynthesis, stress responses, defense responses, transcription factors and cellular components related to ion or protein transportation and signal transduction. Comparing hierarchical structure of biological processes, molecule functions and cellular components revealed the significant regulation differences between TAM 111 and TAM 112, particularly for genes of phosphorylation and adenyl ribonucleotide binding, and proteins located in nucleus and plasma membrane. TAM 112 showed more active than TAM 111 in response to drought and carried more specific genes with most of them were up-regulated in responses to stresses of water deprivation, heat and oxidative, ABA-induced signal pathway and transcription regulation. In addition, 258 genes encoding predicted uncharacterized proteins and 141 unannotated genes with no similar sequences identified in the databases may represent novel genes related to drought response in TAM 111 or TAM 112. This research thus revealed different drought-tolerance mechanisms in TAM 111 and TAM 112 and identified useful drought tolerance genes for wheat adaption. Data of gene sequence and expression regulation from this study also provided useful information of annotating novel genes associated with drought tolerance in the wheat genome.

Genome wide identification of QTL associated with yield and yield components in two popular wheat cultivars TAM 111 and TAM 112.

Two drought-tolerant wheat cultivars, 'TAM 111' and 'TAM 112', have been widely grown in the Southern Great Plains of the U.S. and used as parents in many wheat breeding programs worldwide. This study aimed to reveal genetic control of yield and yield components in the two cultivars under both dryland and irrigated conditions. A mapping population containing 124 F5:7 recombinant inbred lines (RILs) was developed from the cross of TAM 112/TAM 111. A set of 5,948 SNPs from the wheat 90K iSelect array and double digest restriction-site associated DNA sequencing was used to construct high-density genetic maps. Data for yield and yield components were obtained from 11 environments. QTL analyses were performed based on 11 individual environments, across all environments, within and across mega-environments. Thirty-six unique consistent QTL regions were distributed on 13 chromosomes including 1A, 1B, 1D, 2A, 2D, 3D, 4B, 4D, 6A, 6B, 6D, 7B, and 7D. Ten unique QTL with pleiotropic effects were identified on four chromosomes and eight were in common with the consistent QTL. These QTL increased dry biomass grain yield by 16.3 g m-2, plot yield by 28.1 g m-2, kernels spike-1 by 0.7, spikes m-2 by 14.8, thousand kernel weight by 0.9 g with favorable alleles from either parent. TAM 112 alleles mainly increased spikes m-2 and thousand kernel weight while TMA 111 alleles increased kernels spike-1, harvest index and grain yield. The saturated genetic map and markers linked to significant QTL from this study will be very useful in developing high throughput genotyping markers for tracking the desirable haplotypes of these important yield-related traits in popular parental cultivars.

Universal promoter scanning by Pol II during transcription initiation in Saccharomyces cerevisiae.

BACKGROUND: The majority of eukaryotic promoters utilize multiple transcription start sites (TSSs). How multiple TSSs are specified at individual promoters across eukaryotes is not understood for most species. In Saccharomyces cerevisiae, a pre-initiation complex (PIC) comprised of Pol II and conserved general transcription factors (GTFs) assembles and opens DNA upstream of TSSs. Evidence from model promoters indicates that the PIC scans from upstream to downstream to identify TSSs. Prior results suggest that TSS distributions at promoters where scanning occurs shift in a polar fashion upon alteration in Pol II catalytic activity or GTF function. RESULTS: To determine the extent of promoter scanning across promoter classes in S. cerevisiae, we perturb Pol II catalytic activity and GTF function and analyze their effects on TSS usage genome-wide. We find that alterations to Pol II, TFIIB, or TFIIF function widely alter the initiation landscape consistent with promoter scanning operating at all yeast promoters, regardless of promoter class. Promoter architecture, however, can determine the extent of promoter sensitivity to altered Pol II activity in ways that are predicted by a scanning model. CONCLUSIONS: Our observations coupled with previous data validate key predictions of the scanning model for Pol II initiation in yeast, which we term the shooting gallery. In this model, Pol II catalytic activity and the rate and processivity of Pol II scanning together with promoter sequence determine the distribution of TSSs and their usage.

Identification and quantification of coding and long non-coding RNAs in stallion spermatozoa separated by density.

BACKGROUND: It is not unusual for stallions to have fertility problems. For many, artificial insemination with more dense spermatozoa (isolated by density gradient centrifugation) results in greater pregnancy rates compared with the rates when using unfractionated spermatozoa. RNAs in spermatozoa delivered to the oocyte at conception are required for embryo development. Novel molecular assays of spermatozoa that reflect function are needed to predict the fertility of stallions. OBJECTIVES: To describe and compare the RNA populations in more dense and less dense spermatozoa from stallions. MATERIALS AND METHODS: Spermatozoa from five stallions were separated into more dense and less dense populations by density gradient centrifugation. Complementary DNA libraries were made from each of the ten total RNA samples after ribosomal RNA removal. Next-generation sequencing characterized the RNA populations in more and less dense spermatozoa. Quantitative reverse transcription-PCR was used to confirm differential expression of selected RNAs. RESULTS: Stallion spermatozoa contain 11 215 RNAs, with the most prevalent RNA being a 1492 base long non-coding RNA. The levels of 159 RNAs were greater in more dense spermatozoa, while levels of seven other RNAs were greater in less dense spermatozoa. Quantitative reverse transcription-PCR confirmed the threefold greater levels of solute carrier family 26 member 8 (SLC26A8) mRNA in less dense spermatozoa, and sixfold and threefold greater expression levels of the SCP2 sterol binding domain containing 1 (SCP2D1) and spermatogenesis-associated protein 31D1 (SPATA31D1) mRNAs in more dense spermatozoa, respectively. DISCUSSION AND CONCLUSION: We identified 11 215 RNAs in stallion spermatozoa and 166 with differential expression between more dense and less dense fractions. Many prevalent RNAs were also found in bull, boar, and human spermatozoa. Many differentially expressed RNAs are known to be testis- or spermatozoa-specific. Our results may lead to identification of an RNA population in spermatozoa that is optimal for establishing successful pregnancies.

Accurate Genomic Predictions for Chronic Wasting Disease in U.S. White-Tailed Deer.

The geographic expansion of chronic wasting disease (CWD) in U.S. white-tailed deer (Odocoileus virginianus) has been largely unabated by best management practices, diagnostic surveillance, and depopulation of positive herds. Using a custom Affymetrix Axiom single nucleotide polymorphism (SNP) array, we demonstrate that both differential susceptibility to CWD, and natural variation in disease progression, are moderately to highly heritable ([Formula: see text] among farmed U.S. white-tailed deer, and that loci other than PRNP are involved. Genome-wide association analyses using 123,987 quality filtered SNPs for a geographically diverse cohort of 807 farmed U.S. white-tailed deer (n = 284 CWD positive; n = 523 CWD non-detect) confirmed the prion gene (PRNP; G96S) as a large-effect risk locus (P-value < 6.3E-11), as evidenced by the estimated proportion of phenotypic variance explained (PVE ≥ 0.05), but also demonstrated that more phenotypic variance was collectively explained by loci other than PRNP Genomic best linear unbiased prediction (GBLUP; n = 123,987 SNPs) with k-fold cross validation (k = 3; k = 5) and random sampling (n = 50 iterations) for the same cohort of 807 farmed U.S. white-tailed deer produced mean genomic prediction accuracies ≥ 0.81; thereby providing the necessary foundation for exploring a genomically-estimated CWD eradication program.

The Pacific Biosciences de novo assembled genome dataset from a parthenogenetic New Zealand wild population of the longhorned tick, Haemaphysalis longicornis Neumann, 1901.

The longhorned tick, Haemaphysalis longicornis, feeds upon a wide range of bird and mammalian hosts. Mammalian hosts include cattle, deer, sheep, goats, humans, and horses. This tick is known to transmit a number of pathogens causing tick-borne diseases, and was the vector of a recent serious outbreak of oriental theileriosis in New Zealand. A New Zealand-USA consortium was established to sequence, assemble, and annotate the genome of this tick, using ticks obtained from New Zealand's North Island. In New Zealand, the tick is considered exclusively parthenogenetic and this trait was deemed useful for genome assembly. Very high molecular weight genomic DNA was sequenced on the Illumina HiSeq4000 and the long-read Pac Bio Sequel platforms. Twenty-eight SMRT cells produced a total of 21.3 million reads which were assembled with Canu on a reserved supercomputer node with access to 12 TB of RAM, running continuously for over 24 days. The final assembly dataset consisted of 34,211 contigs with an average contig length of 215,205 bp. The quality of the annotated genome was assessed by BUSCO analysis, an approach that provides quantitative measures for the quality of an assembled genome. Over 95% of the BUSCO gene set was found in the assembled genome. Only 48 of the 1066 BUSCO genes were missing and only 9 were present in a fragmented condition. The raw sequencing reads and the assembled contigs/scaffolds are archived at the National Center for Biotechnology Information.

Young and aged TLR4 deficient mice show sex-dependent enhancements in spatial memory and alterations in interleukin-1 related genes.

Toll-like receptor-4 (TLR4) is a transmembrane receptor that initiates an immune response following a bacterial infection or host derived molecules associated with cellular distress. Beyond triggering inflammation, TLR4 has been implicated in modulating behavioral and cognitive processes in a physiologically normal state, as young adult TLR4 deficient mice show learning enhancements in select tasks. Currently unknown is whether these benefits are present in both sexes and persist with aging. The present study evaluated spatial memory, anxiety-like behavior, and central levels of pro- and anti-inflammatory molecules in young (4-5 months) and aged (18-19 months) TLR4 deficient (TLR4-/-) and wild-type (WT) male and female mice. Results confirmed that TLR4-/- mice show enhanced spatial memory compared to WT mice. These effects were age- and sex-specific, as memory retention was superior in the TLR4-/- young males and aged females. While TLR4-/- mice showed age-related changes in behavior, these changes were attenuated relative to aged WT mice. Further, aged TLR4-/- mice showed differential expression of molecules involved in interleukin (IL)-1 signaling in the hippocampus. For instance, aged TLR4-/- females showed heightened expression of IL-1 receptor antagonist (IL-1ra) and the IL-1 accessory proteins AcP and AcPb. Collectively, these data provide the initial evidence that TLR4 deficiency enhances cognitive function and modulates the inflammatory profile of the hippocampus in a sex- and age-dependent manner.

Horse Y chromosome assembly displays unique evolutionary features and putative stallion fertility genes.

Dynamic evolutionary processes and complex structure make the Y chromosome among the most diverse and least understood regions in mammalian genomes. Here, we present an annotated assembly of the male specific region of the horse Y chromosome (eMSY), representing the first comprehensive Y assembly in odd-toed ungulates. The eMSY comprises single-copy, equine specific multi-copy, PAR transposed, and novel ampliconic sequence classes. The eMSY gene density approaches that of autosomes with the highest number of retained X-Y gametologs recorded in eutherians, in addition to novel Y-born and transposed genes. Horse, donkey and mule testis RNAseq reveals several candidate genes for stallion fertility. A novel testis-expressed XY ampliconic sequence class, ETSTY7, is shared with the parasite Parascaris genome, providing evidence for eukaryotic horizontal transfer and inter-chromosomal mobility. Our study highlights the dynamic nature of the Y and provides a reference sequence for improved understanding of equine male development and fertility.