Characterization of the genome and silk-gland transcriptomes of Darwin's bark spider (Caerostris darwini).

Natural silks crafted by spiders comprise some of the most versatile materials known. Artificial silks-based on the sequences of their natural brethren-replicate some desirable biophysical properties and are increasingly utilized in commercial and medical applications today. To characterize the repertoire of protein sequences giving silks their biophysical properties and to determine the set of expressed genes across each unique silk gland contributing to the formation of natural silks, we report here draft genomic and transcriptomic assemblies of Darwin's bark spider, Caerostris darwini, an orb-weaving spider whose dragline is one of the toughest known biomaterials on Earth. We identify at least 31 putative spidroin genes, with expansion of multiple spidroin gene classes relative to the golden orb-weaver, Trichonephila clavipes. We observed substantial sharing of spidroin repetitive sequence motifs between species as well as new motifs unique to C. darwini. Comparative gene expression analyses across six silk gland isolates in females plus a composite isolate of all silk glands in males demonstrated gland and sex-specific expression of spidroins, facilitating putative assignment of novel spidroin genes to classes. Broad expression of spidroins across silk gland types suggests that silks emanating from a given gland represent composite materials to a greater extent than previously appreciated. We hypothesize that the extraordinary toughness of C. darwini major ampullate dragline silk may relate to the unique protein composition of major ampullate spidroins, combined with the relatively high expression of stretchy flagelliform spidroins whose union into a single fiber may be aided by novel motifs and cassettes that act as molecule-binding helices. Our assemblies extend the catalog of sequences and sets of expressed genes that confer the unique biophysical properties observed in natural silks.

Golden orb-weaving spider (Trichonephila clavipes) silk genes with sex-biased expression and atypical architectures.

Spider silks are renowned for their high-performance mechanical properties. Contributing to these properties are proteins encoded by the spidroin (spider fibroin) gene family. Spidroins have been discovered mostly through cDNA studies of females based on the presence of conserved terminal regions and a repetitive central region. Recently, genome sequencing of the golden orb-web weaver, Trichonephila clavipes, provided a complete picture of spidroin diversity. Here, we refine the annotation of T. clavipes spidroin genes including the reclassification of some as non-spidroins. We rename these non-spidroins as spidroin-like (SpL) genes because they have repetitive sequences and amino acid compositions like spidroins, but entirely lack the archetypal terminal domains of spidroins. Insight into the function of these spidroin and SpL genes was then examined through tissue- and sex-specific gene expression studies. Using qPCR, we show that some silk genes are upregulated in male silk glands compared to females, despite males producing less silk in general. We also find that an enigmatic spidroin that lacks a spidroin C-terminal domain is highly expressed in silk glands, suggesting that spidroins could assemble into fibers without a canonical terminal region. Further, we show that two SpL genes are expressed in silk glands, with one gene highly evolutionarily conserved across species, providing evidence that particular SpL genes are important to silk production. Together, these findings challenge long-standing paradigms regarding the evolutionary and functional significance of the proteins and conserved motifs essential for producing spider silks.

Benchmarking Single-Cell mRNA-Sequencing Technologies Uncovers Differences in Sensitivity and Reproducibility in Cell Types With Low RNA Content.

Single-cell RNA sequencing (scRNA-seq) has the ability to classify each cell and determine the transcriptomic profile of specific cell types and cells of a given disease state; however, sensitivity of the gene count for each cell can be a critical component to the success of a single-cell study. The recently introduced SMART-Seq Single Cell PLUS Kit (SSsc PLUS) claims to provide higher sensitivity and reproducibility versus popular methods for the sequencing analysis of single cells. Here, the cDNA-generation component of the kit, SMART-Seq Single Cell Kit (SSsc), was compared with the popular homebrew protocol, Smart-seq2, and its update, Smart-seq3. The SMART-Seq Library Prep Kit from SSsc PLUS was benchmarked against a commonly used scRNA-seq library preparation method, Illumina Nextera XT. Finally, the SSsc chemistry was tested in both full and fractional volumes on 2 popular liquid-handler devices to investigate whether the high sensitivity was maintained in miniaturization. We demonstrate that SSsc PLUS outperforms these other full-length methods in convenience, sensitivity, gene identification, and reproducibility while also offering full compatibility with automation platforms.

Multiplexed Targeted Resequencing Identifies Coding and Regulatory Variation Underlying Phenotypic Extremes of High-Density Lipoprotein Cholesterol in Humans.

BACKGROUND: Genome-wide association studies have uncovered common variants at many loci influencing human complex traits, such as high-density lipoprotein cholesterol (HDL-C). However, the contribution of the identified genes is difficult to ascertain from current efforts interrogating common variants with small effects. Thus, there is a pressing need for scalable, cost-effective strategies for uncovering causal variants, many of which may be rare and noncoding. METHODS: Here, we used a molecular inversion probe target capture approach to resequence both coding and regulatory regions at 7 HDL-C-associated loci in 797 individuals with extremely high HDL-C versus 735 low-to-normal HDL-C controls. Our targets included protein-coding regions of GALNT2, APOA5, APOC3, SCARB1, CCDC92, ZNF664, CETP, and LIPG (>9 kb) and proximate noncoding regulatory features (>42 kb). RESULTS: Exome-wide genotyping in 1114 of the 1532 participants yielded a >90% genotyping concordance rate with molecular inversion probe-identified variants in ≈90% of participants. This approach rediscovered nearly all established genome-wide association studies associations in GALNT2, CETP, and LIPG loci with significant and concordant associations with HDL-C from our phenotypic extremes design at 0.1% of the sample size of lipid genome-wide association studies. In addition, we identified a novel, rare, CETP noncoding variant enriched in the extreme high HDL-C group (P<0.01, score test). CONCLUSIONS: Our targeted resequencing of individuals at the HDL-C phenotypic extremes offers a novel, efficient, and cost-effective approach for identifying rare coding and noncoding variation differences in extreme phenotypes and supports the rationale for applying this methodology to uncover rare variation-particularly noncoding variation-underlying myriad complex traits.

The Nephila clavipes genome highlights the diversity of spider silk genes and their complex expression.

Spider silks are the toughest known biological materials, yet are lightweight and virtually invisible to the human immune system, and they thus have revolutionary potential for medicine and industry. Spider silks are largely composed of spidroins, a unique family of structural proteins. To investigate spidroin genes systematically, we constructed the first genome of an orb-weaving spider: the golden orb-weaver (Nephila clavipes), which builds large webs using an extensive repertoire of silks with diverse physical properties. We cataloged 28 Nephila spidroins, representing all known orb-weaver spidroin types, and identified 394 repeated coding motif variants and higher-order repetitive cassette structures unique to specific spidroins. Characterization of spidroin expression in distinct silk gland types indicates that glands can express multiple spidroin types. We find evidence of an alternatively spliced spidroin, a spidroin expressed only in venom glands, evolutionary mechanisms for spidroin diversification, and non-spidroin genes with expression patterns that suggest roles in silk production.

Progesterone Response Element Variation in the OXTR Promoter Region and Paternal Care in New World Monkeys.

Paternal care is a complex social behavior common in primate species with socially monogamous mating systems and twin births. Evolutionary causes and consequences of such behavior are not well understood, nor are their neuroendocrine and genetic bases. However, the neuropeptide oxytocin (OXT) and its receptor (OXTR) are associated with parental care in mammalian lineages. Here we investigated the interspecific variation in the number of progesterone response elements (PREs) in the OXTR promoter region of 32 primate species, correlating genetic data with behavior, social systems, and ecological/life-history parameters, while controlling for phylogeny. We verified that PREs are only present in New World monkeys and that PRE number is significantly correlated with the presence of paternal care in this branch. We suggest that PRE number could be an essential part of the genetic repertoire that allowed the emergence of taxon-specific complex social behaviors, such as paternal care in marmosets and tamarins.

Oxytocin receptor gene sequences in owl monkeys and other primates show remarkable interspecific regulatory and protein coding variation.

The oxytocin (OT) hormone pathway is involved in numerous physiological processes, and one of its receptor genes (OXTR) has been implicated in pair bonding behavior in mammalian lineages. This observation is important for understanding social monogamy in primates, which occurs in only a small subset of taxa, including Azara's owl monkey (Aotus azarae). To examine the potential relationship between social monogamy and OXTR variation, we sequenced its 5' regulatory (4936bp) and coding (1167bp) regions in 25 owl monkeys from the Argentinean Gran Chaco, and examined OXTR sequences from 1092 humans from the 1000 Genomes Project. We also assessed interspecific variation of OXTR in 25 primate and rodent species that represent a set of phylogenetically and behaviorally disparate taxa. Our analysis revealed substantial variation in the putative 5' regulatory region of OXTR, with marked structural differences across primate taxa, particularly for humans and chimpanzees, which exhibited unique patterns of large motifs of dinucleotide A+T repeats upstream of the OXTR 5' UTR. In addition, we observed a large number of amino acid substitutions in the OXTR CDS region among New World primate taxa that distinguish them from Old World primates. Furthermore, primate taxa traditionally defined as socially monogamous (e.g., gibbons, owl monkeys, titi monkeys, and saki monkeys) all exhibited different amino acid motifs for their respective OXTR protein coding sequences. These findings support the notion that monogamy has evolved independently in Old World and New World primates, and that it has done so through different molecular mechanisms, not exclusively through the oxytocin pathway.

Correlates of genetic monogamy in socially monogamous mammals: insights from Azara's owl monkeys.

Understanding the evolution of mating systems, a central topic in evolutionary biology for more than 50 years, requires examining the genetic consequences of mating and the relationships between social systems and mating systems. Among pair-living mammals, where genetic monogamy is extremely rare, the extent of extra-group paternity rates has been associated with male participation in infant care, strength of the pair bond and length of the breeding season. This study evaluated the relationship between two of those factors and the genetic mating system of socially monogamous mammals, testing predictions that male care and strength of pair bond would be negatively correlated with rates of extra-pair paternity (EPP). Autosomal microsatellite analyses provide evidence for genetic monogamy in a pair-living primate with bi-parental care, the Azara's owl monkey (Aotus azarae). A phylogenetically corrected generalized least square analysis was used to relate male care and strength of the pair bond to their genetic mating system (i.e. proportions of EPP) in 15 socially monogamous mammalian species. The intensity of male care was correlated with EPP rates in mammals, while strength of pair bond failed to reach statistical significance. Our analyses show that, once social monogamy has evolved, paternal care, and potentially also close bonds, may facilitate the evolution of genetic monogamy.

An optimized microsatellite genotyping strategy for assessing genetic identity and kinship in Azara's owl monkeys (Aotus azarai).

In this study, we characterize a panel of 20 microsatellite markers that reproducibly amplify in Azara's owl monkeys (Aotus azarai) for use in genetic profiling analyses. A total of 128 individuals from our study site in Formosa, Argentina, were genotyped for 20 markers, 13 of which were found to be polymorphic. The levels of allelic variation at these loci provided paternity exclusion probabilities of 0.852 when neither parent was known, and 0.981 when one parent was known. In addition, our analysis revealed that, although genotypes can be rapidly scored using fluorescence-based fragment analysis, the presence of complex or multiple short tandem repeat (STR) motifs at a microsatellite locus could generate similar fragment patterns from alleles that have different nucleotide sequences and perhaps different evolutionary origins. Even so, this collection of microsatellite loci is suitable for parentage analyses and will allow us to test various hypotheses about the relationship between social behavior and kinship in wild owl monkey populations. Furthermore, given the limited number of platyrrhine-specific microsatellite loci available in the literature, this STR panel represents a valuable tool for population studies of other cebines and callitrichines.

mtDNA diversity in Azara's owl monkeys (Aotus azarai azarai) of the Argentinean Chaco.

Owl monkeys (Aotus spp.) inhabit much of South America yet represent an enigmatic evolutionary branch among primates. While morphological, cytogenetic, and immunological evidence suggest that owl monkey populations have undergone isolation and diversification since their emergence in the New World, problems with adjacent species ranges, and sample provenance have complicated efforts to characterize genetic variation within the genus. As a result, the phylogeographic history of owl monkey species and subspecies remains unclear, and the extent of genetic diversity at the population level is unknown. To explore these issues, we analyzed mitochondrial DNA (mt DNA) variation in a population of wild Azara's owl monkeys (Aotus azarai azarai) living in the Gran Chaco region of Argentina. We sequenced the complete mitochondrial genome from one individual (16,585 base pairs (bp)) and analyzed 1,099 bp of the hypervariable control region (CR) and 696 bp of the cytochrome oxidase II (COII) gene in 117 others. In addition, we sequenced the mitochondrial genome (16,472 bp) of one Nancy Ma's owl monkey (A. nancymaae). Based on the whole mtDNA and COII data, we observed an ancient phylogeographic discontinuity among Aotus species living north, south, and west of the Amazon River that began more than eight million years ago. Our population analyses identified three major CR lineages and detected a high level of haplotypic diversity within A. a. azarai. These data point to a recent expansion of Azara's owl monkeys into the Argentinean Chaco. Overall, we provide a detailed view of owl monkey mtDNA variation at genus, species, and population levels.

Refinement of primate copy number variation hotspots identifies candidate genomic regions evolving under positive selection.

BACKGROUND: Copy number variants (CNVs), defined as losses and gains of segments of genomic DNA, are a major source of genomic variation. RESULTS: In this study, we identified over 2,000 human CNVs that overlap with orthologous chimpanzee or orthologous macaque CNVs. Of these, 170 CNVs overlap with both chimpanzee and macaque CNVs, and these were collapsed into 34 hotspot regions of CNV formation. Many of these hotspot regions of CNV formation are functionally relevant, with a bias toward genes involved in immune function, some of which were previously shown to evolve under balancing selection in humans. The genes in these primate CNV formation hotspots have significant differential expression levels between species and show evidence for positive selection, indicating that they have evolved under species-specific, directional selection. CONCLUSIONS: These hotspots of primate CNV formation provide a novel perspective on divergence and selective pressures acting on these genomic regions.

AVPR1A sequence variation in monogamous owl monkeys (Aotus azarai) and its implications for the evolution of platyrrhine social behavior.

The arginine vasopressin V1a receptor gene (AVPR1A) has been implicated in increased partner preference and pair bonding behavior in mammalian lineages. This observation is of considerable importance for studies of social monogamy, which only appears in a small subset of primate taxa, including the Argentinean owl monkey (Aotus azarai). Thus, to investigate the possible influence of AVPR1A on the evolution of social behavior in owl monkeys, we sequenced this locus in a wild population from the Gran Chaco. We also assessed the interspecific variation of AVPR1A in platyrrhine species that represent a set of phylogenetically and behaviorally disparate taxa. The resulting data revealed A. azarai to have a unique genic structure for AVPR1A that varies in coding sequence and microsatellite repeat content relative to other primate and mammalian species. Specifically, one repetitive region that has been the focus in studies of human AVPR1A diversity, "RS3," is completely absent in A. azarai and all other platyrrhines examined. This finding suggests that, if AVPR1A modulates behavior in owl monkeys and other neotropical primates, it does so independent of this region. These observations have also provided clues about the process by which the range of social behavior in the Order Primates evolved through lineage-specific neurogenetic variation.

Expanding Southwest Pacific mitochondrial haplogroups P and Q.

Modern humans have occupied New Guinea and the nearby Bismarck and Solomon archipelagos of Island Melanesia for at least 40,000 years. Previous mitochondrial DNA (mtDNA) studies indicated that two common lineages in this region, haplogroups P and Q, were particularly diverse, with the coalescence for P considered significantly older than that for Q. In this study, we expand the definition of haplogroup Q so that it includes three major branches, each separated by multiple mutational distinctions (Q1, equivalent to the earlier definition of Q, plus Q2 and Q3). We report three whole-mtDNA genomes that establish Q2 as a major Q branch. In addition, we describe 314 control region sequences that belong to the expanded haplogroups P and Q from our Southwest Pacific collection. The coalescence dates for the largest P and Q branches (P1 and Q1) are similar to each other (approximately 50,000 years old) and considerably older than prior estimates. Newly identified Q2, which was found in Island Melanesian samples just to the east, is somewhat younger by more than 10,000 years. Our coalescence estimates should be more reliable than prior ones because they were based on significantly larger samples as well as complete mtDNA-coding region sequencing. Our estimates are roughly in accord with the current suggested dates for the first settlement of New Guinea-Sahul. The phylogeography of P and Q indicates almost total (female) isolation of ancient New Guinea-Island Melanesia from Australia that may have existed from the time of the first settlement. While Q subsequently diversified extensively in New Guinea-Island Melanesia, it has not been found in Australia. The only shared mtDNA haplogroup between Australia and New Guinea identified to date remains one minor branch of P.

Mitochondrial DNA variation and the origins of the Aleuts.

The mitochondrial DNA (mtDNA) variation in 179 Aleuts from five different islands (Atka, Unalaska, Umnak, St. Paul, and St. George) and Anchorage was analyzed to better understand the origins of Aleuts and their role in the peopling of the Americas. Mitochondrial DNA samples were characterized using polymerase chain reaction amplification, restriction fragment length polymorphism analysis, and direct sequencing of the first hypervariable segment (HVS-I) of the control region. This study showed that Aleut mtDNAs belonged to two of the four haplogroups (A and D) common among Native Americans. Haplogroup D occurred at a very high frequency in Aleuts, and this, along with their unique HVS-I sequences, distinguished them from Eskimos, Athapaskan Indians, and other northern Amerindian populations. While sharing several control region sequences (CIR11, CHU14, CIR60, and CIR61) with other circumarctic populations, Aleuts lacked haplogroup A mtDNAs having the 16265G mutation that are specific to Eskimo populations. R-matrix and median network analyses indicated that Aleuts were closest genetically to Chukotkan (Chukchi and Siberian Eskimos) rather than to Native American or Kamchatkan populations (Koryaks and Itel'men). Dating of the Beringian branch of haplogroup A (16192T) suggested that populations ancestral to the Aleuts, Eskimos, and Athapaskan Indians emerged approximately 13,120 years ago, while Aleut-specific A and D sublineages were dated at 6539 +/- 3511 and 6035 +/- 2885 years, respectively. Our findings support the archaeologically based hypothesis that ancestral Aleuts crossed the Bering Land Bridge or Beringian platform and entered the Aleutian Islands from the east, rather than island hopping from Kamchatka into the western Aleutians. Furthermore, the Aleut migration most likely represents a separate event from those responsible for peopling the remainder of the Americas, meaning that the New World was colonized through multiple migrations.